| 查看: 442 | 回复: 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. 这个问题怎么解决？？？？ |
回复此楼» 猜你喜欢
一志愿重庆大学085700资源与环境总分308求调剂
已经有6人回复
-
268求调剂
已经有7人回复
-
302求调剂
已经有11人回复
-
306求调剂
已经有4人回复
-
材料学学硕080502 337求调剂-一志愿华中科技大学
已经有4人回复
-
求调剂
已经有6人回复
-
313求调剂
已经有3人回复
-
南昌大学材料专硕311分求调剂
已经有6人回复
-
316求调剂
已经有6人回复
-
346求调剂[0856]
已经有7人回复
【答案】应助回帖
感谢参与，应助指数 +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
