24小时热门版块排行榜    

查看: 238  |  回复: 2

岸芷汀兰-

铁虫 (小有名气)

[交流] 【求助】画图时又出错了 已有2人参与

直接看程序。。。。
function ctn
y0=[0;0];
%tspan=[0:2.8202*10^(-5)/100:2.8202*10^(-5)];
tspan=[0:2.8202*10^(-6)/100:2.8202*10^(-6)];
[t,y]=ode45('funb',tspan,y0,[]);

x=0:0.33/100:0.33;
a=length(t);
b=length(x);
for i=1:a
for j=1:b
z(j,i)=2.7786*sin(2*pi/0.33*x(j)).*y(i,1);
end
end
mesh(t,x,z);
xlabel('t'),ylabel('x'),zlabel('z')


function equ1=funb(t,y,flag)
equ1=[y(2);(-1.6561*10^5*y(2)-3.2257*10^10*y(1)+(3.5454e+008./0.1295)*2.7786*sin(pi/0.33*8.2980*10^10*t^2))]




function cttttn
w0=[0;0];
%tspan=[0:2.8202*10^(-5)/100:2.8202*10^(-5)];
tspan=[0:2.8202*10^(-6)/100:2.8202*10^(-6)];
[t,w]=ode45('funa',tspan,w0,[]);

x=0:0.33/100:0.33;
a=length(t);
b=length(x);
for i=1:a
for j=1:b
z(j,i)= 2.7786*sin(pi/0.33*x(j)).*w(i,1);
end
end
mesh(t,x,z);
xlabel('t'),ylabel('x'),zlabel('z')


function equ1=funa(t,w,flag)
equ1=[w(2);(-1.6561*10^5*w(2)-3.2255*10^10*w(1)+(3.5454e+008./0.2028)*2.7786*sin(pi/0.33*0.5*8.2980*10^10*t^2))]
上下是两个程序,但是我想把这个图画到一起,就是让两个z(j,i)叠加,变量不变,我试了试,不出图,也不显示哪儿出错,就出来一堆数。。。请问如何弄呢?
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

eyouwcy

木虫 (正式写手)

★ ★
robert2020(金币+1):多谢提供建议! 2010-05-29 08:21:38
zzuwangshilei(金币+1):多谢提供建议 2010-05-29 13:59:05
岸芷汀兰-(金币+1):谢谢。。。 2010-05-30 09:33:24
岸芷汀兰-(金币+1): 2010-05-31 10:51:21
看过,是否Z下标有误。检查一下看看
2楼2010-05-29 01:14:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

岸芷汀兰-

铁虫 (小有名气)

zzuwangshilei:积极关注自己的帖子,欢迎具体指出弄混的变量 2010-05-30 11:26:29
其中有个变量给弄混了,,,已解决。。。
3楼2010-05-30 09:34:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 岸芷汀兰- 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见