| 查看: 497 | 回复: 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. 这个问题怎么解决？？？？ |
回复此楼» 猜你喜欢
面上的水很深
已经有21人回复
-
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有15人回复
-
国社科通讯评审
已经有5人回复
-
5个%2F，非固定段22个字母
已经有9人回复
-
生命口C13面上会评时间
已经有11人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
伯胺上甲基求助
已经有7人回复
-
跨出版社商投稿
已经有5人回复
-
生生命学部会评
已经有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