| 查看: 506 | 回复: 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. 这个问题怎么解决？？？？ |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有80人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有14人回复
-
感觉是下周放榜了
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
filecode=后面第一个是大写字母
已经有8人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
【答案】应助回帖
感谢参与，应助指数 +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; |
» 本帖已获得的红花（最新10朵）
寻常2楼2016-08-15 23:59:13
3楼2016-08-16 10:25:32
25