| 查看: 432 | 回复: 0 | ||
[求助]
求助高斯积分的程序
|
|
图中的式子如果用高斯积分计算最后的结果，程序运行一直错误，各位大神帮我看一下该怎么改程序 function y=gaussf(x) m=0; sym x; for j=1:inf m= m+(1/j)*(sin(3.14*j)-sin(2.86*j))*cos(3.14*j)*exp(-0.00001247*x*j^2 end; y=exp(-0.21*x-((15060-0.127*x)^2)/(2.286*x))*(0.089+0.6369*m))/(x^(-1.5)); function s=IntGauss(a,b) ta=(b-a)/2; tb=(a+b)/2; s=ta*(0.2955242247*subs(sym(gaussf),findsym(sym(gaussf)),ta*0.148874339+tb)+0.2955242247*subs(sym(gaussf),findsym(sym(gaussf)),-ta*0.148874339+tb)+0.2692667193*subs(sym(gaussf),findsym(sym(gaussf)),ta*0.4333953941+tb)+0.2692667193*subs(sym(gaussf),findsym(sym(gaussf)),-ta*0.4333953941+tb)+0.2190863625*subs(sym(gaussf),findsym(sym(gaussf)),ta*0.6794095683+tb)+0.2190863625*subs(sym(gaussf),findsym(sym(gaussf)),-ta*0.6794095683+tb)+0.1494513492*subs(sym(gaussf),findsym(sym(gaussf)),ta*0.8650633667+tb)+0.1494513492*subs(sym(gaussf),findsym(sym(gaussf)),-ta*0.8650633667+tb)+0.0666713443*subs(sym(gaussf),findsym(sym(gaussf)),ta*0.9739065285+tb)+0.0666713443*subs(sym(gaussf),findsym(sym(gaussf)),-ta*0.9739065285+tb)); 运行结果的程序： for i=1:1000 t=365*i; s(i)=IntGauss(0,t); d(i)=349903.14*s(i); end 错误提示： ??? Input argument "x" is undefined. Error in ==> gaussf at 5 y=exp(-0.21*x-((15060-0.127*x)^2)/(2.286*x))*(0.089+0.6369*(m+(1/j)*(sin(3.14*j)-sin(2.86*j))*cos(3.14*j)*exp(-0.00001247*x*j^2)))/(x^(-1.5)); Error in ==> IntGauss at 4 s=ta*(0.2955242247*subs(sym(gaussf),findsym(sym(gaussf)),ta*0.148874339+tb)+0.2955242247*subs(sym(gaussf),findsym(sym(gaussf)),-  Z6$0HC1X}X4_4{`XWU6R683.jpg |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
找到一些相关的精华帖子，希望有用哦~
求助高斯文件输入
已经有5人回复
- 高斯方程积分编程matlab或C语言
已经有15人回复
- 一个函数与高斯函数的乘积的积分问题
已经有8人回复
- 求助高斯输出文件中坐标系怎么定义的？
已经有11人回复
- 求助高斯拟合的问题
已经有7人回复
- 新手求助高斯输入文件写法，求大神帮忙
已经有8人回复
- 关于delta函数的高斯积分问题
已经有32人回复
- 新手—高斯仿真红外光谱
已经有31人回复
- 多重高斯积分fortran程序
已经有4人回复
- 【求助】帮忙解释一下高斯执行计算的命令
已经有3人回复
- 【求助】高斯计算需要温度的反应
已经有12人回复
- 【求助】高斯运行问题
已经有14人回复
- 【求助】高斯chk文件如何用Gaussiview 打开
已经有20人回复
- 【求助】高斯频率计算死掉的原因
已经有7人回复
- 【求助】高斯计算出错
已经有11人回复
- 【求助】高斯b3lyp/VDZ的输入文件写法
已经有14人回复
- 【求助】求高斯计算磁矩？
已经有5人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
40