| 查看: 985 | 回复: 2 | ||
[求助]
matlab中想对一个公式进行积分，但是公式中还含有其它变量 已有1人参与
|
|
如题，公式[z(j)./(1*Uz(j,t))./(1+50.*n(t).*z(j)./(2*pi*Uz(j,t))).^(5/3)].*exp(-[abs(dy)/60]-[abs(dz)/10])中好多变量，其中z(j)为坐标，Uz(j,t)为t时刻j点的速度，dy，dz为两点坐标差值，n(t)为频率，但我想把上式对n，即频率积分，以下是我的脚本 for i=1:20 for j=i:20 for t=1:500 dy=y(i)-y(j); dz=z(i)-z(j); R0(i,j,t)=int([200*1.76^2.*z(j)./(1*Uz(j,t))./(1+50.*n(t).*z(j)./(2*pi*Uz(j,t))).^(5/3)].*exp(-[abs(dy)/60]-[abs(dz)/10]),0,10);%H0后是频率上下限 end end end 但是一直报错，之前用的quadl函数也不行，下面是报错内容，上面积分那块好像是没有声明积分对象，只给了上下限，求大神指点 尝试引用非结构体数组的字段。 出错 isAllVars (line 9) res = strcmp(mupadmex('symobj::isAllVars',expr.s,0),'TRUE'); 出错 sym/int (line 150) if ~isscalar(x) || ~isAllVars(x) |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有14人回复
-
这样的filecode谁见过
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
waterque
金虫 (正式写手)
- 应助: 59 (初中生)
- 金币: 1989.1
- 散金: 96
- 红花: 10
- 帖子: 759
- 在线: 191.9小时
- 虫号: 2095273
- 注册: 2012-10-30
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2018-06-29 22:36:13
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2018-06-29 22:36:13
|
可能是你的积分表达式中有多个变量，所以需要说明对哪个变量积分，积分上下限是哪个变量的上下限。具体你可以参考matlab官网的说明： https://www.mathworks.com/help/symbolic/int.html |
2楼2018-06-26 09:40:31
3楼2018-06-29 09:08:20
20