PDA

查看完整版本 : [MATLAB图像处理] 【求助】关于用MATLAB或者Origin画图的问题 膜拜【啊】


d_kelvin
2009-05-26, 14:56
N取值0到51间隔1 k取值0到1间隔0.1 画出以k为横坐标E为纵坐标的图

函数如下

E(N,k)=2.6*sqrt(3+2cos(3.14*N*(6+8)/459.9216-(sqrt(3)/2)*(6/(2.46*sqrt(76)))*k*2.46)+2cos(3.14*N*(12+4)/459.9216+(sqrt(3)/2)*(4/(2.46*sqrt(76)))*k*2.46)+2cos(3.14*N*(2/459.9216)+(sqrt(3)/2)*(10/(2.46*sqrt(76)))*k*2.46))

求助 编程~~~~膜拜高人~~~

TTT_IOU
2009-05-26, 21:23
下列程序可否:
一.
syms N k
f=2.6*sqrt(3+2*cos(3.14*N*(6+8)/459.9216-(sqrt(3)/2)*(6/(2.46*sqrt(76)))*k*2.46)+2*cos(3.14*N*(12+4)/459.9216+(sqrt(3)/2)*(4/(2.46*sqrt(76)))*k*2.46)+2*cos(3.14*N*(2/459.9216)+(sqrt(3)/2)*(10/(2.46*sqrt(76)))*k*2.46));
f1=subs(f,N,[0:51]);
f2=vpa(f1,5);
for i=1:length(f2);
figure
ezplot(f2(i))
end
%k不去范围
二.
syms N k
f=2.6*sqrt(3+2*cos(3.14*N*(6+8)/459.9216-(sqrt(3)/2)*(6/(2.46*sqrt(76)))*k*2.46)+...
2*cos(3.14*N*(12+4)/459.9216+(sqrt(3)/2)*(4/(2.46*sqrt(76)))*k*2.46)+...
2*cos(3.14*N*(2/459.9216)+(sqrt(3)/2)*(10/(2.46*sqrt(76)))*k*2.46));
f1=subs(f,N,[0:51]);
f2=vpa(f1,5);
for i=1:length(f2);
figure
ezplot(f2(i),[0 ,1])
end
希望有帮助!!!

d_kelvin
2009-05-26, 22:02
下列程序可否:
一.
syms N k
f=2.6*sqrt(3+2*cos(3.14*N*(6+8)/459.9216-(sqrt(3)/2)*(6/(2.46*sqrt(76)))*k*2.46)+2*cos(3.14*N*(12+4)/459.9216+(sqrt(3)/2)*(4/(2.46*sqrt(76)))*k*2.46)+2*cos(3.14*N*(2/459.9216)+(sqrt(3)/2)*(10/(2.46*sqrt(76)))*k*2.46));
f1=subs(f,N,[0:51]);
f2=vpa(f1,5);
for i=1:length(f2);
figure
ezplot(f2(i))
end
%k不去范围
二.
syms N k
f=2.6*sqrt(3+2*cos(3.14*N*(6+8)/459.9216-(sqrt(3)/2)*(6/(2.46*sqrt(76)))*k*2.46)+...
2*cos(3.14*N*(12+4)/459.9216+(sqrt(3)/2)*(4/(2.46*sqrt(76)))*k*2.46)+...
2*cos(3.14*N*(2/459.9216)+(sqrt(3)/2)*(10/(2.46*sqrt(76)))*k*2.46));
f1=subs(f,N,[0:51]);
f2=vpa(f1,5);
for i=1:length(f2);
figure
ezplot(f2(i),[0,1])
end
希望有帮助!!!


谢谢高人~~~~
但是 我运行了下 它出现了52个图~~~如何把这52个图画到一个坐标上显示呢~~~:)

TTT_IOU
2009-05-26, 23:32
syms N k
f=2.6*sqrt(3+2*cos(3.14*N*(6+8)/459.9216-(sqrt(3)/2)*(6/(2.46*sqrt(76)))*k*2.46)+...
2*cos(3.14*N*(12+4)/459.9216+(sqrt(3)/2)*(4/(2.46*sqrt(76)))*k*2.46)+...
2*cos(3.14*N*(2/459.9216)+(sqrt(3)/2)*(10/(2.46*sqrt(76)))*k*2.46));
f1=subs(f,N,[0:51]);
f2=vpa(f1,5);
for i=1:length(f2);
ezplot(f2(i),[0,1])
hold on
end
不知上面程序是否是你要的,望有帮助!!!