| 查看: 251 | 回复: 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)叠加，变量不变，我试了试，不出图，也不显示哪儿出错，就出来一堆数。。。请问如何弄呢？ |
回复此楼» 猜你喜欢
材料专硕(0856) 339分求调剂
已经有7人回复
-
081700学硕，323分，一志愿中国海洋大学求调剂学校
已经有15人回复
-
295求调剂
已经有10人回复
-
一志愿9材料学硕297已过六级求调剂推荐
已经有10人回复
-
312求调剂
已经有3人回复
-
一志愿北京化工085600 310分求调剂
已经有12人回复
-
一志愿 南京航空航天大学 ，080500材料科学与工程学硕
已经有8人回复
-
工科277分求调剂材料
已经有5人回复
-
电子信息调剂交叉学科有推荐吗
已经有9人回复
-
材料工程专硕求调剂
已经有11人回复
3楼2010-05-30 09:34:07
eyouwcy
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4505.3
- 红花: 3
- 帖子: 658
- 在线: 87.4小时
- 虫号: 260996
- 注册: 2006-06-18
- 性别: GG
- 专业: 金属材料的微观结构
2楼2010-05-29 01:14:11
