| 查看: 453 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有19人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有24人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有10人回复
找到一些相关的精华帖子，希望有用哦~
求助高斯文件输入
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
50