24小时热门版块排行榜    

查看: 845  |  回复: 2

Math露珠

木虫 (知名作家)


[交流] Error using plot Vectors must be the same length.如下程序如何更改?

Error using plot Vectors must be the same length.如下程序如何更改?


hold on %figureaxes
>> q1=0.94;q2=0.94;q3=0.94; %x,y,z
>> h=0.05;T=250;  %hT
>> a=0.45;b=0.75; %
>> p1=16.5;p2=17.5;  %pp1p2
>> m1=1;m2=1;  %
>> x0=0.1;y0=0.001;z0=0.001;  %
>> x0=1;y0=1;z0=2; %
>> p=p1; %p%----  -----
>> x(1)=x0+h^q1*y0;
>> y(1)=y0+h^q2*(x0-b*y0-x0*z0);
>> z(1)=z0+h^q3*(-a*z0+x0*x0);
>> t=m1*h;n=1;  %
>> % ------  --------
>> while t<=T  %
    n=n+1;   %
    M1=0;M2=0;M3=0;   %
  % xyz
    for k=1:n-1
        M1=M1+(-1)^(k)*gamma(q1+1)/(gamma(k+1)*gamma(q1-k+1))*(x(n-k)-x0);
        M2=M2+(-1)^(k)*gamma(q2+1)/(gamma(k+1)*gamma(q2-k+1))*(y(n-k)-y0);
        M3=M3+(-1)^(k)*gamma(q3+1)/(gamma(k+1)*gamma(q3-k+1))*(z(n-k)-z0);
    end
   %5m1*hm2*h
    if p==p1
      p=p2;t=t+m2*h;
    else
        p=p1;t=t+m1*h;
    end
%
   x(n)=x0-M1+h^q1*y(n-1);
   y(n)=y0-M2+h^q2*(x(n)-b*y(n-1)-x(n)*z(n-1));  %pp1p2
   z(n)=z0-M3+h^q3*(-a*z(n-1)+x(n)*x(n));
end  %
>> t=0:h:T; %
>> figure(13)  %x(t)
>> plot(t,x,'K');xlabel('t'),ylabel('x');

错误提示:
Error using plot
Vectors must be the same length.

提供帮助的虫子50金币奖励,谢谢大家
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
M塞北2楼
2017-10-24 05:09   回复  
Math露珠(金币+1): 谢谢参与
发自小木虫Android客户端
dsctg3楼
2017-10-24 08:40   回复  
Math露珠(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 Math露珠 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 i3KBRD19ERkG 2026-09-16 4/200 2026-09-17 07:50 by k06BKNblrGRH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 5BDX0d0WFp7t 2026-09-15 3/150 2026-09-17 06:50 by k06BKNblrGRH
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 mibUvS8DDCwf 2026-09-15 4/200 2026-09-17 06:26 by k06BKNblrGRH
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:38 by k06BKNblrGRH
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 5/250 2026-09-17 05:14 by k06BKNblrGRH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-13 6/300 2026-09-17 02:26 by Tql5LQhh5rLK
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 21:14 by MHf40yA0Uhyb
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-16 19:15 by rRvxEj3zorCy
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 16:00 by i3KBRD19ERkG
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:48 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:36 by LH5NkK5Ud4bh
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:50 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
信息提示
请填处理意见