版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (2663)
>虫友互识 (305)
>导师招生 (112)
>考博 (108)
>文献求助 (80)
>基金申请 (67)
>硕博家园 (56)
>仿真模拟 (44)
>博后之家 (43)
>论文投稿 (43)
>休闲灌水 (39)
>教师之家 (30)
>公派出国 (30)
>找工作 (29)
>考研 (27)
>招聘信息布告栏 (21)

返回列表

Math露珠

木虫 (知名作家)

应助: 10 (幼儿园)
金币: 7725.3
帖子: 8879
在线: 886小时
虫号: 4010226

[交流] 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金币奖励，谢谢大家

回复此楼

» 猜你喜欢

河北省自然科学基金已经有8人回复
西安交大新媒学院副院长用撤稿论文结题已经有5人回复
论文撤稿了已经有5人回复
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！已经有5人回复
26/27申博自荐已经有9人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有3人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有7人回复
揭秘青基评审内幕：几个A才能顺利中标已经有4人回复
青B发送上会通知了吗已经有7人回复
博士申请已经有3人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2017-10-24 03:03:34

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

简单回复

M塞北2楼

2017-10-24 05:09 回复

Math露珠(金币+1): 谢谢参与

发自小木虫Android客户端

dsctg3楼

2017-10-24 08:40 回复