【求助】关于用MATLAB或者Origin画图的问题 膜拜【啊】
 

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))

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

回复: 【求助】关于用MATLAB或者Origin画图的问题 膜拜【啊】
 

下列程序可否：
一.
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
希望有帮助！！！

回复: 【求助】关于用MATLAB或者Origin画图的问题 膜拜【啊】
 

[QUOTE=TTT_IOU;28655]下列程序可否：
一.
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
希望有帮助！！！[/QUOTE]


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

回复: 【求助】关于用MATLAB或者Origin画图的问题 膜拜【啊】
 

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
不知上面程序是否是你要的，望有帮助！！！