| 查看: 443 | 回复: 2 | ||
[求助]
MATLAB编程问题 已有1人参与
|
|
dt=0.0001; t=0:dt:0.5; R=1000000;S=58.0644*10^(-4);B=8.8541878*10^(-12);v=0.1;d=125/3400000;A=10*10^(-6);a=1;n=0.05; z1=exp((1/R*S*B)*(d.*t+1/6*a.*t.^3)); z2=exp((1/R*S*B)*(d*t+0.5*sqrt(a*n).*t.^2-1/6*a.*t.^3)); I1=-(A*d/(R*B))+((A*(d+0.5*a.t^2))/(R*B))*exp((-1/(R*S*B))*(d.*t+(1/6)*a.*t.^3))+(A*d/(R*B))*((d+0.5*a.*t.^2)/(R*B*S))*exp(-(1/(R*S*B))*(d*t+(1/6)*a.*t.^3)).*cumtrapz(z1)*dt; I2=-(A*d/(R*B))+((A*(d+sqrt(a*n).*t-0.5*a.*t.^2))/(R*B))*exp((-1/(R*S*B))*(d.*t+0.5*sqrt(a*n).*t^2-(1/6)*a.*t.^3))+(A*d/(R*B))*((d+sqrt(a*n).*t-0.5*a.*t.^2)/(R*B*S))*exp(-(1/(R*S*B))*(0.5*sqrt(a*n).*t^2-(1/6)*a.*t.^3)).*cumtrapz(z2)*dt; y=I1.*(t<=sqrt(n/a)&t>=0)+I2.*(t<=2*sqrt(n/a)&t>sqrt(n/a)); plot(t,y) Attempt to reference field of non-structure array. 这个问题怎么解决？？？？ |
回复此楼» 猜你喜欢
初试 317
已经有3人回复
-
279分求调剂 一志愿211
已经有15人回复
-
286求调剂
已经有4人回复
-
279求调剂
已经有4人回复
-
求调剂
已经有3人回复
-
求调剂
已经有3人回复
-
299求调剂
已经有3人回复
-
330求调剂0854
已经有3人回复
-
299求调剂
已经有6人回复
-
286分人工智能专业请求调剂愿意跨考！
已经有5人回复
【答案】应助回帖
感谢参与，应助指数 +1
|
I1，I2里面很多格式错误，重点搞清楚什么时候要矩阵点运算，下面代码是修改过的I1，I2： I1=-(A*d/(R*B))+((A*(d+0.5*a.*t.^2))/(R*B)).*exp((-1/(R*S*B))*(d.*t+(1/6)*a.*t.^3))+(A*d/(R*B))*((d+0.5*a.*t.^2)/(R*B*S)).*exp(-(1/(R*S*B))*(d*t+(1/6)*a.*t.^3)).*cumtrapz(z1)*dt; I2=-(A*d/(R*B))+((A*(d+sqrt(a*n).*t-0.5*a.*t.^2))/(R*B)).*exp((-1/(R*S*B))*(d.*t+0.5*sqrt(a*n).*t.^2-(1/6)*a.*t.^3))+(A*d/(R*B))*((d+sqrt(a*n).*t-0.5*a.*t.^2)/(R*B*S)).*exp(-(1/(R*S*B))*(0.5*sqrt(a*n).*t.^2-(1/6)*a.*t.^3)).*cumtrapz(z2)*dt; |
寻常2楼2016-08-15 23:59:13
3楼2016-08-16 10:25:32
