| 查看: 897 | 回复: 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) |
回复此楼» 猜你喜欢
球磨粉体时遇到了大的问题，请指教！
已经有15人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有5人回复
-
江汉大学解明教授课题组招博士研究生/博士后
已经有3人回复
3楼2018-06-29 09:08:20
waterque
金虫 (正式写手)
- 应助: 59 (初中生)
- 金币: 1989.1
- 散金: 96
- 红花: 10
- 帖子: 759
- 在线: 191.7小时
- 虫号: 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