24小时热门版块排行榜    

查看: 333  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 岸芷汀兰- 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 19:50 by xzcv0Tehd8lQ
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 17:01 by BCdvBAEpt3zm
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 16:58 by BCdvBAEpt3zm
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 16:58 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 16:45 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 15:01 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 14:50 by BCdvBAEpt3zm
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 14:37 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 12:38 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 10:43 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 10:27 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
请填处理意见