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

登录注册

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

24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者summer_ziye将赠送您 5 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

summer_ziye

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1404
散金: 20
红花: 1
帖子: 261
在线: 273.7小时
虫号: 1446398
注册: 2011-10-17
专业: 控制理论与方法

[求助] 求各位大神帮忙我的程序老是提示错误如何修改已有2人参与

clc
clear

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%% Example 2
A11=[0 1; -44.75 -20]; A12=[0 1; -41.19 -20];
A21=[0 1; -42.55 -20]; A22=[0 1; -38.99 -20];
Atau11=[0 0; -1 -0.4]; Atau12=Atau11;
Atau21=[0 0; -1.2 -0.5]; Atau22=Atau21;
B12=[0 0; 0.8 0]; B21=[0 0; 1 0];
C11=[0; 0.5]; C12=C11; C21=[0; 0.4]; C22=C21;
D11=[1 0]; D12=D11; D21=D11; D22=D11;
Dtau11=[0.1 0]; Dtau12=Dtau11;
Dtau21=[0.08 0]; Dtau22=Dtau21;
G11=[0 1]; G12=G11; G21=G11; G22=G11; E11=0; E12=0; E21=0; E22=0;

%gamma=1;  %%%%% for simulation
%%%%%%%%%%%%%

Af11=[-2.4409 0.7508; -193.4481  -36.0634];Af12=[ -2.4704 0.7714; -196.6598  -36.5871];
Af21=[-1.1664 0.8471; -163.0604  -34.7050];Af22=[-1.1757 0.8727; -166.4750  -35.1970];
Bf11=[-1.8255; -137.2169];Bf12=[-1.8850; -145.3440];
Bf21=[-0.6479; -113.0800];Bf22=[-0.6873; -121.2615];
Gf11=[0.0788 -0.4708];Gf12=[0.0431 -0.4466];
Gf21=[0.0825 -0.4523];Gf22=[0.0454 -0.4286];
%%%%%%%%%%%%%

A111=[A11, zeros(2); Bf11*D11, Af11]; A112=[A11, zeros(2); Bf12*D11, Af12];
A121=[A12, zeros(2); Bf11*D12, Af11]; A122=[A12, zeros(2); Bf12*D12, Af12];
Atau111=[Atau11, zeros(2); Bf11*Dtau11, zeros(2)]; Atau112=[Atau11, zeros(2); Bf12*Dtau11, zeros(2)];
Atau121=[Atau12, zeros(2); Bf11*Dtau12, zeros(2)]; Atau122=[Atau12, zeros(2); Bf12*Dtau12, zeros(2)];
C111=[C11; Bf11*E11]; C112=[C11; Bf12*E11];
C121=[C12; Bf11*E12]; C122=[C12; Bf12*E12];
G111=[G11, -Gf11]; G112=[G11, -Gf12]; G121=[G12, -Gf11]; G122=[G12, -Gf12];
M21=[B21, zeros(2); zeros(2), zeros(2)];

A211=[A21, zeros(2); Bf21*D21, Af21]; A212=[A21, zeros(2); Bf22*D21, Af22];
A221=[A22, zeros(2); Bf21*D22, Af21]; A222=[A22, zeros(2); Bf22*D22, Af22];
Atau211=[Atau21, zeros(2); Bf21*Dtau21, zeros(2)]; Atau212=[Atau21, zeros(2); Bf22*Dtau21, zeros(2)];
Atau221=[Atau22, zeros(2); Bf21*Dtau22, zeros(2)]; Atau222=[Atau22, zeros(2); Bf22*Dtau22, zeros(2)];
C211=[C21; Bf21*E21]; C212=[C21; Bf22*E21];
C221=[C22; Bf21*E22]; C222=[C22; Bf22*E22];
G211=[G21, -Gf21]; G212=[G21, -Gf22]; G221=[G22, -Gf21]; G222=[G22, -Gf22];
M12=[B12, zeros(2); zeros(2), zeros(2)];

ee1=[1, 0, 0, 0]; ee2=[0, 1, 0, 0]; ee3=[0, 0, 1, 0]; ee4=[0, 0, 0, 1];

x101=0.3; x102=0; x103=0; x104=0; x10=[x101 x102 x103 x104]';
phi1=x10;
x201=0.3; x202=0; x203=0; x204=0; x20=[x201 x202 x203 x204]';
phi2=x20;

h=0.002;
finaltime=15;

%%%%%%%%%%%%%%%%%%%%%%%% initial value of membership function h_i

h11(1:1,1)=abs(0.637*x101);
  h12(1:1,1)=1-h11(1:1,1);
  ttau1(1:1,1)=0.2*sin(0)-0.1;
  tau1(1:1,1)=(fix(ttau1(1:1,1)/h))*h;

h21(1:1,1)=abs(0.637*x201);
h22(1:1,1)=1-h21(1:1,1);
  ttau2(1:1,1)=0.2*sin(0)-0.1;
  tau2(1:1,1)=(fix(ttau2(1:1,1)/h))*h;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

dx10=(h11(1:1,1)*h11(1:1,1)*A111+h11(1:1,1)*h12(1:1,1)*A112+h12(1:1,1)*h11(1:1,1)*A121+h12(1:1,1)*h12(1:1,1)*A122)*x10+(h11(1:1,1)*h11(1:1,1)*Atau111+h11(1:1,1)*h12(1:1,1)*Atau112+h12(1:1,1)*h11(1:1,1)*Atau121+h12(1:1,1)*h12(1:1,1)*Atau122)*phi1+M21*phi2+(h11(1:1,1)*h11(1:1,1)*C111+h11(1:1,1)*h12(1:1,1)*C112+h12(1:1,1)*h11(1:1,1)*C121+h12(1:1,1)*h12(1:1,1)*C122)*0;
%dx10 is the derivative of x1 at time 0
x1(1:4,1)=dx10*h+x10;
%x1(i) is the value of x1 at step(i)

dx20=(h21(1:1,1)*h21(1:1,1)*A211+h21(1:1,1)*h22(1:1,1)*A212+h22(1:1,1)*h21(1:1,1)*A221+h22(1:1,1)*h22(1:1,1)*A222)*x20+(h21(1:1,1)*h21(1:1,1)*Atau211+h21(1:1,1)*h22(1:1,1)*Atau212+h22(1:1,1)*h21(1:1,1)*Atau221+h22(1:1,1)*h22(1:1,1)*Atau222)*phi2+M12*phi1+(h21(1:1,1)*h21(1:1,1)*C211+h21(1:1,1)*h22(1:1,1)*C212+h22(1:1,1)*h21(1:1,1)*C221+h22(1:1,1)*h22(1:1,1)*C222)*0;
%dx20 is the derivative of x2 at time 0
x2(1:4,1)=dx20*h+x20;
%x2(i) is the value of x2 at step(i)

v1(1:1,1)=0;
v2(1:1,1)=0;
e1(1:2,1)=(h11(1:1,1)*h11(1:1,1)*G111+h11(1:1,1)*h12(1:1,1)*G112+h12(1:1,1)*h11(1:1,1)*G121+h12(1:1,1)*h12(1:1,1)*G122)*x1(1:4,1);
e2(1:2,1)=(h21(1:1,1)*h21(1:1,1)*G211+h21(1:1,1)*h22(1:1,1)*G212+h22(1:1,1)*h21(1:1,1)*G221+h22(1:1,1)*h22(1:1,1)*G222)*x2(1:4,1);
vplus(1:1,1)=h*((v1(1:1,1))^2+(v2(1:1,1))^2);
eplus(1:1,1)=h*((e1(1:2,1))'*e1(1:2,1)+(e2(1:2,1))'*e2(1:2,1));
ev(1:1,1)=eplus(1:1,1)/vplus(1:1,1);

for count=1:finaltime/h
t=count*h;

h11(:,count)=abs(0.637*ee1*x1(:,count));
h12(:,count)=1-h11(:,count);

h21(:,count)=abs(0.637*ee1*x2(:,count));
h22(:,count)=1-h21(:,count);

  v1(:,count)=(sin(4*pi*t))*exp(-0.5*t);
  v2(:,count)=(sin(4*pi*t))*exp(-0.5*t);

  ttau1(:,count)=0.2*sin(t)-0.1;
%  tau1(:,count)=(fix(ttau1(:,count)/h)+1)*h;

  ttau2(:,count)=0.2*sin(t)-0.1;
%  tau2(:,count)=(fix(ttau2(:,count)/h)+1)*h;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

  if t<=ttau1(:,count)

dx1(:,count)=(h11(:,count)*h11(:,count)*A111+h11(:,count)*h12(:,count)*A112+h12(:,count)*h11(:,count)*A121+h12(:,count)*h12(:,count)*A122)*x1(:,count)+(h11(:,count)*h11(:,count)*Atau111+h11(:,count)*h12(:,count)*Atau112+h12(:,count)*h11(:,count)*Atau121+h12(:,count)*h12(:,count)*Atau122)*phi1+M21*x2(:,count)+(h11(:,count)*h11(:,count)*C111+h11(:,count)*h12(:,count)*C112+h12(:,count)*h11(:,count)*C121+h12(:,count)*h12(:,count)*C122)*v1(:,count);
x1(:,count+1)=dx1(:,count)*h+x1(:,count);
else


dx1(:,count)=(h11(:,count)*h11(:,count)*A111+h11(:,count)*h12(:,count)*A112+h12(:,count)*h11(:,count)*A121+h12(:,count)*h12(:,count)*A122)*x1(:,count)+(h11(:,count)*h11(:,count)*Atau111+h11(:,count)*h12(:,count)*Atau112+h12(:,count)*h11(:,count)*Atau121+h12(:,count)*h12(:,count)*Atau122)*x1(:,count-fix(ttau1(:,count)/h))+M21*x2(:,count)+(h11(:,count)*h11(:,count)*C111+h11(:,count)*h12(:,count)*C112+h12(:,count)*h11(:,count)*C121+h12(:,count)*h12(:,count)*C122)*v1(:,count);
x1(:,count+1)=dx1(:,count)*h+x1(:,count);
  end


  if t<=ttau2(:,count)

dx2(:,count)=(h21(:,count)*h21(:,count)*A211+h21(:,count)*h22(:,count)*A212+h22(:,count)*h21(:,count)*A221+h22(:,count)*h22(:,count)*A222)*x2(:,count)+(h21(:,count)*h21(:,count)*Atau211+h21(:,count)*h22(:,count)*Atau212+h22(:,count)*h21(:,count)*Atau221+h22(:,count)*h22(:,count)*Atau222)*phi2+M12*x1(:,count)+(h21(:,count)*h21(:,count)*C211+h21(:,count)*h22(:,count)*C212+h22(:,count)*h21(:,count)*C221+h22(:,count)*h22(:,count)*C222)*v2(:,count);
x2(:,count+1)=dx2(:,count)*h+x2(:,count);
else


dx2(:,count)=(h21(:,count)*h21(:,count)*A211+h21(:,count)*h22(:,count)*A212+h22(:,count)*h21(:,count)*A221+h22(:,count)*h22(:,count)*A222)*x2(:,count)+(h21(:,count)*h21(:,count)*Atau211+h21(:,count)*h22(:,count)*Atau212+h22(:,count)*h21(:,count)*Atau221+h22(:,count)*h22(:,count)*Atau222)*x2(:,count-fix(ttau2(:,count)/h))+M12*x1(:,count)+(h21(:,count)*h21(:,count)*C211+h21(:,count)*h22(:,count)*C212+h22(:,count)*h21(:,count)*C221+h22(:,count)*h22(:,count)*C222)*v2(:,count);
x2(:,count+1)=dx2(:,count)*h+x2(:,count);
  end

e1(:,count)=(h11(:,count)*h11(:,count)*G111+h11(:,count)*h12(:,count)*G112+h12(:,count)*h11(:,count)*G121+h12(:,count)*h12(:,count)*G122)*x1(:,count);
e2(:,count)=(h21(:,count)*h21(:,count)*G211+h21(:,count)*h22(:,count)*G212+h22(:,count)*h21(:,count)*G221+h22(:,count)*h22(:,count)*G222)*x2(:,count);
vplus(:,count+1)=h*((v1(:,count))^2+(v2(:,count))^2)+vplus(:,count);
eplus(:,count+1)=h*((e1(:,count))'*e1(:,count)+(e2(:,count))'*e2(:,count))+eplus(:,count);
ev(:,count+1)=eplus(:,count+1)/vplus(:,count+1);

end

%%%%%%%%%%%%%plot xi(1),xi(2), e1,e2, in the same fig

hold on
plot([0:h:finaltime],x1(2,

)
plot([0:h:finaltime],x1(2,

-e1(

,'--')
%plot([0:h:finaltime],x2(2,

)
%plot([0:h:finaltime],x2(2,

-e2(

,'--')
%plot([0:h:finaltime],e1(

)
%plot([0:h:finaltime],e2(

,'--')
hold off
%legend('z_1(t)','z_{f1}(t)'))
%legend('z_2(t)','z_{f2}(t)'))
%legend('e_1(t)','e_2(t)'))

%hold on
%plot([0:h:finaltime],ev(

)
%hold off
%legend('u(t)')
xlabel('Time t')

??? Index exceeds matrix dimensions.
本人是初学者，实在不会，求各位大神帮帮忙

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

matlab程序报错，请大神帮忙改正已经有3人回复
求大神帮忙解释截断能测试中的问题！已经有7人回复
求助大神！DNA-PAGE一直出现问题，不知道什么原因，请求大神回复帮忙解决。已经有15人回复
ansys有错误怎么办啊大神们帮帮忙已经有6人回复
求大神帮忙解题啊！！！已经有4人回复
求大神们帮忙看看这个程序哪儿需要改啊已经有7人回复
额，最近老是被这个问题困扰，大神们，求助啊！！！！已经有8人回复
这些错误该怎么改啊？求大神赐教已经有3人回复
求助各位大神我的情况。。。。。。。。已经有18人回复
菜鸟求助，老师给了一个分子让优化，可总是错误，求大神指点已经有22人回复
求大神帮我解答，谢谢，必定重谢！！！已经有4人回复
各位大神，请帮我看看题目解法对不对已经有6人回复
各位大神前辈为什么小弟在做模拟时总是出现这样的错误，求指导一下，不胜感激已经有4人回复
求大神看一下我这个，感觉是软件设置的问题已经有7人回复
求大神们帮我解决一个关于修改系统进化树字体的问题已经有18人回复
求大神帮忙解答！已经有8人回复
我电脑出了点小问题求大神们帮帮忙已经有8人回复
有大神帮忙看看这三个题目么？已经有3人回复
帮忙改一下这个A，B的错误已经有6人回复
求助大神帮忙解答如下题目已经有5人回复
修改图标问题，请各位大神帮忙！已经有5人回复
求助！这个程序错在哪里？已经有8人回复
求助大神帮我看看这段程序是干什么的？（就是定义这个函数算什么）已经有6人回复
求各位大神一个问题已经有4人回复
请大神帮忙看一下问题所在谢谢已经有4人回复

1楼2014-06-01 16:31:52

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

summer_ziye

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1404
散金: 20
红花: 1
帖子: 261
在线: 273.7小时
虫号: 1446398
注册: 2011-10-17
专业: 控制理论与方法

谢谢我的问题已经解决了，请斑竹封贴

赞一下

回复此楼

6楼2014-06-11 13:28:25

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

查看全部 6 个回答

yangter92

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 690.2
散金: 60
帖子: 99
在线: 34.5小时
虫号: 3182467
注册: 2014-05-05
性别: GG
专业: 计算数学与科学工程计算

【答案】应助回帖

感谢参与，应助指数 +1

你把前面有百分号的先去掉，太多了看着都晕

[ 发自小木虫客户端 ]

赞一下

回复此楼

justdomyself

2楼2014-06-01 17:28:50

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

summer_ziye

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1404
散金: 20
红花: 1
帖子: 261
在线: 273.7小时
虫号: 1446398
注册: 2011-10-17
专业: 控制理论与方法

引用回帖:

2楼: Originally posted by yangter92 at 2014-06-01 17:28:50
你把前面有百分号的先去掉，太多了看着都晕

好的

回复此楼

3楼2014-06-02 19:14:13

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

summer_ziye

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1404
散金: 20
红花: 1
帖子: 261
在线: 273.7小时
虫号: 1446398
注册: 2011-10-17
专业: 控制理论与方法

clc
clear
A11=[0 1; -44.75 -20]; A12=[0 1; -41.19 -20];
A21=[0 1; -42.55 -20]; A22=[0 1; -38.99 -20];
Atau11=[0 0; -1 -0.4]; Atau12=Atau11;
Atau21=[0 0; -1.2 -0.5]; Atau22=Atau21;
B12=[0 0; 0.8 0]; B21=[0 0; 1 0];
C11=[0; 0.5]; C12=C11; C21=[0; 0.4]; C22=C21;
D11=[1 0]; D12=D11; D21=D11; D22=D11;
Dtau11=[0.1 0]; Dtau12=Dtau11;
Dtau21=[0.08 0]; Dtau22=Dtau21;
G11=[0 1]; G12=G11; G21=G11; G22=G11; E11=0; E12=0; E21=0; E22=0;

Af11=[-2.4409 0.7508; -193.4481  -36.0634];Af12=[ -2.4704 0.7714; -196.6598  -36.5871];
Af21=[-1.1664 0.8471; -163.0604  -34.7050];Af22=[-1.1757 0.8727; -166.4750  -35.1970];
Bf11=[-1.8255; -137.2169];Bf12=[-1.8850; -145.3440];
Bf21=[-0.6479; -113.0800];Bf22=[-0.6873; -121.2615];
Gf11=[0.0788 -0.4708];Gf12=[0.0431 -0.4466];
Gf21=[0.0825 -0.4523];Gf22=[0.0454 -0.4286];

A111=[A11, zeros(2); Bf11*D11, Af11]; A112=[A11, zeros(2); Bf12*D11, Af12];
A121=[A12, zeros(2); Bf11*D12, Af11]; A122=[A12, zeros(2); Bf12*D12, Af12];
Atau111=[Atau11, zeros(2); Bf11*Dtau11, zeros(2)]; Atau112=[Atau11, zeros(2); Bf12*Dtau11, zeros(2)];
Atau121=[Atau12, zeros(2); Bf11*Dtau12, zeros(2)]; Atau122=[Atau12, zeros(2); Bf12*Dtau12, zeros(2)];
C111=[C11; Bf11*E11]; C112=[C11; Bf12*E11];
C121=[C12; Bf11*E12]; C122=[C12; Bf12*E12];
G111=[G11, -Gf11]; G112=[G11, -Gf12]; G121=[G12, -Gf11]; G122=[G12, -Gf12];
M21=[B21, zeros(2); zeros(2), zeros(2)];

A211=[A21, zeros(2); Bf21*D21, Af21]; A212=[A21, zeros(2); Bf22*D21, Af22];
A221=[A22, zeros(2); Bf21*D22, Af21]; A222=[A22, zeros(2); Bf22*D22, Af22];
Atau211=[Atau21, zeros(2); Bf21*Dtau21, zeros(2)]; Atau212=[Atau21, zeros(2); Bf22*Dtau21, zeros(2)];
Atau221=[Atau22, zeros(2); Bf21*Dtau22, zeros(2)]; Atau222=[Atau22, zeros(2); Bf22*Dtau22, zeros(2)];
C211=[C21; Bf21*E21]; C212=[C21; Bf22*E21];
C221=[C22; Bf21*E22]; C222=[C22; Bf22*E22];
G211=[G21, -Gf21]; G212=[G21, -Gf22]; G221=[G22, -Gf21]; G222=[G22, -Gf22];
M12=[B12, zeros(2); zeros(2), zeros(2)];

ee1=[1, 0, 0, 0]; ee2=[0, 1, 0, 0]; ee3=[0, 0, 1, 0]; ee4=[0, 0, 0, 1];

x101=0.3; x102=0; x103=0; x104=0; x10=[x101 x102 x103 x104]';
phi1=x10;
x201=0.3; x202=0; x203=0; x204=0; x20=[x201 x202 x203 x204]';
phi2=x20;

h=0.002;
finaltime=15;

h11(1:1,1)=abs(0.637*x101);
  h12(1:1,1)=1-h11(1:1,1);
  ttau1(1:1,1)=0.2*sin(0)-0.1;
  tau1(1:1,1)=(fix(ttau1(1:1,1)/h))*h;

h21(1:1,1)=abs(0.637*x201);
h22(1:1,1)=1-h21(1:1,1);
  ttau2(1:1,1)=0.2*sin(0)-0.1;
  tau2(1:1,1)=(fix(ttau2(1:1,1)/h))*h;

dx10=(h11(1:1,1)*h11(1:1,1)*A111+h11(1:1,1)*h12(1:1,1)*A112+h12(1:1,1)*h11(1:1,1)*A121+h12(1:1,1)*h12(1:1,1)*A122)*x10+(h11(1:1,1)*h11(1:1,1)*Atau111+h11(1:1,1)*h12(1:1,1)*Atau112+h12(1:1,1)*h11(1:1,1)*Atau121+h12(1:1,1)*h12(1:1,1)*Atau122)*phi1+M21*phi2+(h11(1:1,1)*h11(1:1,1)*C111+h11(1:1,1)*h12(1:1,1)*C112+h12(1:1,1)*h11(1:1,1)*C121+h12(1:1,1)*h12(1:1,1)*C122)*0;

x1(1:4,1)=dx10*h+x10;

dx20=(h21(1:1,1)*h21(1:1,1)*A211+h21(1:1,1)*h22(1:1,1)*A212+h22(1:1,1)*h21(1:1,1)*A221+h22(1:1,1)*h22(1:1,1)*A222)*x20+(h21(1:1,1)*h21(1:1,1)*Atau211+h21(1:1,1)*h22(1:1,1)*Atau212+h22(1:1,1)*h21(1:1,1)*Atau221+h22(1:1,1)*h22(1:1,1)*Atau222)*phi2+M12*phi1+(h21(1:1,1)*h21(1:1,1)*C211+h21(1:1,1)*h22(1:1,1)*C212+h22(1:1,1)*h21(1:1,1)*C221+h22(1:1,1)*h22(1:1,1)*C222)*0;

x2(1:4,1)=dx20*h+x20;

v1(1:1,1)=0;
v2(1:1,1)=0;
e1(1:2,1)=(h11(1:1,1)*h11(1:1,1)*G111+h11(1:1,1)*h12(1:1,1)*G112+h12(1:1,1)*h11(1:1,1)*G121+h12(1:1,1)*h12(1:1,1)*G122)*x1(1:4,1);
e2(1:2,1)=(h21(1:1,1)*h21(1:1,1)*G211+h21(1:1,1)*h22(1:1,1)*G212+h22(1:1,1)*h21(1:1,1)*G221+h22(1:1,1)*h22(1:1,1)*G222)*x2(1:4,1);
vplus(1:1,1)=h*((v1(1:1,1))^2+(v2(1:1,1))^2);
eplus(1:1,1)=h*((e1(1:2,1))'*e1(1:2,1)+(e2(1:2,1))'*e2(1:2,1));
ev(1:1,1)=eplus(1:1,1)/vplus(1:1,1);

for count=1:finaltime/h
t=count*h;

h11(:,count)=abs(0.637*ee1*x1(:,count));
h12(:,count)=1-h11(:,count);

h21(:,count)=abs(0.637*ee1*x2(:,count));
h22(:,count)=1-h21(:,count);

  v1(:,count)=(sin(4*pi*t))*exp(-0.5*t);
  v2(:,count)=(sin(4*pi*t))*exp(-0.5*t);

  ttau1(:,count)=0.2*sin(t)-0.1;

  ttau2(:,count)=0.2*sin(t)-0.1;

  if t<=ttau1(:,count)

dx1(:,count)=(h11(:,count)*h11(:,count)*A111+h11(:,count)*h12(:,count)*A112+h12(:,count)*h11(:,count)*A121+h12(:,count)*h12(:,count)*A122)*x1(:,count)+(h11(:,count)*h11(:,count)*Atau111+h11(:,count)*h12(:,count)*Atau112+h12(:,count)*h11(:,count)*Atau121+h12(:,count)*h12(:,count)*Atau122)*phi1+M21*x2(:,count)+(h11(:,count)*h11(:,count)*C111+h11(:,count)*h12(:,count)*C112+h12(:,count)*h11(:,count)*C121+h12(:,count)*h12(:,count)*C122)*v1(:,count);
x1(:,count+1)=dx1(:,count)*h+x1(:,count);
else


dx1(:,count)=(h11(:,count)*h11(:,count)*A111+h11(:,count)*h12(:,count)*A112+h12(:,count)*h11(:,count)*A121+h12(:,count)*h12(:,count)*A122)*x1(:,count)+(h11(:,count)*h11(:,count)*Atau111+h11(:,count)*h12(:,count)*Atau112+h12(:,count)*h11(:,count)*Atau121+h12(:,count)*h12(:,count)*Atau122)*x1(:,count-fix(ttau1(:,count)/h))+M21*x2(:,count)+(h11(:,count)*h11(:,count)*C111+h11(:,count)*h12(:,count)*C112+h12(:,count)*h11(:,count)*C121+h12(:,count)*h12(:,count)*C122)*v1(:,count);
x1(:,count+1)=dx1(:,count)*h+x1(:,count);
  end


  if t<=ttau2(:,count)

dx2(:,count)=(h21(:,count)*h21(:,count)*A211+h21(:,count)*h22(:,count)*A212+h22(:,count)*h21(:,count)*A221+h22(:,count)*h22(:,count)*A222)*x2(:,count)+(h21(:,count)*h21(:,count)*Atau211+h21(:,count)*h22(:,count)*Atau212+h22(:,count)*h21(:,count)*Atau221+h22(:,count)*h22(:,count)*Atau222)*phi2+M12*x1(:,count)+(h21(:,count)*h21(:,count)*C211+h21(:,count)*h22(:,count)*C212+h22(:,count)*h21(:,count)*C221+h22(:,count)*h22(:,count)*C222)*v2(:,count);
x2(:,count+1)=dx2(:,count)*h+x2(:,count);
else


dx2(:,count)=(h21(:,count)*h21(:,count)*A211+h21(:,count)*h22(:,count)*A212+h22(:,count)*h21(:,count)*A221+h22(:,count)*h22(:,count)*A222)*x2(:,count)+(h21(:,count)*h21(:,count)*Atau211+h21(:,count)*h22(:,count)*Atau212+h22(:,count)*h21(:,count)*Atau221+h22(:,count)*h22(:,count)*Atau222)*x2(:,count-fix(ttau2(:,count)/h))+M12*x1(:,count)+(h21(:,count)*h21(:,count)*C211+h21(:,count)*h22(:,count)*C212+h22(:,count)*h21(:,count)*C221+h22(:,count)*h22(:,count)*C222)*v2(:,count);
x2(:,count+1)=dx2(:,count)*h+x2(:,count);
  end

e1(:,count)=(h11(:,count)*h11(:,count)*G111+h11(:,count)*h12(:,count)*G112+h12(:,count)*h11(:,count)*G121+h12(:,count)*h12(:,count)*G122)*x1(:,count);
e2(:,count)=(h21(:,count)*h21(:,count)*G211+h21(:,count)*h22(:,count)*G212+h22(:,count)*h21(:,count)*G221+h22(:,count)*h22(:,count)*G222)*x2(:,count);
vplus(:,count+1)=h*((v1(:,count))^2+(v2(:,count))^2)+vplus(:,count);
eplus(:,count+1)=h*((e1(:,count))'*e1(:,count)+(e2(:,count))'*e2(:,count))+eplus(:,count);
ev(:,count+1)=eplus(:,count+1)/vplus(:,count+1);

end

hold on
plot([0:h:finaltime],x1(2,)
plot([0:h:finaltime],x1(2,-e1(,'--')
hold off

xlabel('Time t')

??? Index exceeds matrix dimensions.
本人是初学者，实在不会，求各位大神帮帮忙

赞一下

回复此楼

4楼2014-06-02 19:16:45

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

查看全部 6 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[论文投稿] 好运来！ +3	jibeidaxue 2024-10-15	3/150	2024-10-18 07:46 by lizhengke06
[有机交流] 产物呈油状怎么处理呀 5+3	许嵩YYDS 2024-10-15	3/150	2024-10-18 06:14 by xiaomei1031
[基金申请] 广东面上基金状态可查 +5	堡主大人 2024-10-17	8/400	2024-10-17 23:29 by 19033
[硕博家园] 25级博士申请（仅限同济，汽车机械交通计算机）个人原因，勿喷 +4	亚历山大1 2024-10-17	8/400	2024-10-17 23:05 by 亚历山大1
[基金申请] 科研人在困境中挣扎：如何在新政策下保住辛苦成果，成功变更项目依托单位？ 520+3	雨花阁主 2024-10-12	15/750	2024-10-17 21:40 by asdew124
[硕博家园] 工作不适应 +7	青萍之沫 2024-10-17	10/500	2024-10-17 17:34 by zeolitess
[基金申请] 博士后有基金我可以理解，但是博士，连本科生都申请国自然，我确实理解不了 +19	bychangming 2024-10-12	20/1000	2024-10-17 17:18 by ypp5181
[论文投稿] 老师压着论文不让投咋办？ +12	nvizqxuf 2024-10-16	20/1000	2024-10-17 16:49 by jurkat.1640
[基金申请] 武大14名本科生获自然科学基金青年学生项目资助 +18	babu2015 2024-10-13	20/1000	2024-10-17 16:41 by LNP@mRNA
[考博] 博士自荐 +3	栗子栗子冰 2024-10-16	3/150	2024-10-17 14:47 by 香麦奶茶
[论文投稿] 中文核心期刊。要7000发表非 .各位大佬，是不是不值得？ 6+8	rfmems 2024-10-16	9/450	2024-10-17 13:31 by chuecust
[硕博家园] 研究生避雷川大低碳/建环刘某某 +11	Allllice 2024-10-12	15/750	2024-10-17 09:38 by Allllice
[有机交流] 有什么分析手段可以辨别这两个物质 10+3	YUNQIQI 2024-10-14	8/400	2024-10-16 22:14 by templefire
[基金申请] 中国青年女科学家奖得主，优青，青拔，西北大学刘建妮被基金委查处 +10	babu2015 2024-10-14	13/650	2024-10-16 11:06 by gwjxiaolang
[考博] 学生自荐 +6	Heartwarming 2024-10-12	16/800	2024-10-16 08:40 by 824282658
[电化学] 微型电化学传感器 +4	528104109 2024-10-14	8/400	2024-10-15 16:42 by 528104109
[论文投稿] 2024.10.12 TIMC reject +4	香瓜木香 2024-10-14	6/300	2024-10-15 16:18 by yuanjijoy
[论文投稿] Applied Microbiology and Biotechnology 现在要收费了吗？记得以前不收啊 4+3	jizhihong 2024-10-13	5/250	2024-10-14 19:41 by jizhihong
[教师之家] 。。。。。。。。。 (金币-5) +19	lice_1987 2024-10-14	33/1650	2024-10-14 17:54 by lice_1987
[有机交流] 求助反应过程 +3	xianfeng2895 2024-10-11	4/200	2024-10-12 17:15 by xiaobenpu

24小时热门版块排行榜

summer_ziye

[求助] 求各位大神帮忙 我的程序老是提示错误 如何修改已有2人参与

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

summer_ziye

yangter92

【答案】应助回帖

summer_ziye

summer_ziye

[求助] 求各位大神帮忙我的程序老是提示错误如何修改已有2人参与