| 查看: 496 | 回复: 3 | ||
[求助]
用matlab求解积分问题
|
|
用matlab编了一段程序无法运行，是不是求积分语句不对能，求高手指点。 clear; clc; w=10.0*1e-9; syms Theta fid=fopen('Sw.dat','w'); E=30.0; for n=1:1:800; w=w+1*1e-9 wa=w; wb=w; ky=E*sin(Theta)*1.5156e6; ka=(E-60.0)*1.5156e6; kb=E*1.5156e6; qa=-sqrt(ka^2-ky^2); qb=sqrt(kb^2-ky^2); T1=[cos(qa*wa+Theta)/cos(-Theta),i*sin(qa*wa)/cos(-Theta);i*sin(qa*wa)/cos(-Theta),cos(qa*wa-Theta)/cos(-Theta)]; T2=[cos(qb*wb-Theta)/cos(Theta),i*sin(qb*wb)/cos(Theta);i*sin(qb*wb)/cos(Theta),cos(qb*wb+Theta)/cos(Theta)]; T=T2*T1*T2*T1*T2*T1*T2*T1*T2*T1*T2*T1*T2*T1*T2*T1*T2*T1*T2*T1*T2*T1*T2*T1; y=2.0*cos(Theta)/(T(2,2)*exp(-i*Theta)+T(1,1)*exp(i*Theta)-T(1,2)-T(2,1)); Trans=abs(y)^2; Tx=Trans*cos(Theta); fun=inline(eval(Tx)); Conduct=quad(fun,-0.49*pi,0.49*pi); x=[n,Conduct]; fprintf(fid,'%f %f\n',x); end fclose(fid); |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有27人回复
-
所感
已经有3人回复
-
要不要辞职读博？
已经有7人回复
-
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 由入门到精通--matlab求解高等数学问题课件
已经有100人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 【求助】matlab求解一个简单的积分（急！！！）
已经有7人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【讨论】CFD软件是求解PDE的，matlab中也可以求解PDE
已经有6人回复
- 【求助】用matlab编过有限差分法的请进
已经有9人回复
枫桥夜过
木虫 (小有名气)
铁杆会员
- 应助: 0 (幼儿园)
- 金币: 2017.3
- 帖子: 187
- 在线: 59.2小时
- 虫号: 2121417
- 注册: 2012-11-12
- 性别: GG
- 专业: 岩土与基础工程
2楼2013-04-02 22:12:30
枫桥夜过
木虫 (小有名气)
铁杆会员
- 应助: 0 (幼儿园)
- 金币: 2017.3
- 帖子: 187
- 在线: 59.2小时
- 虫号: 2121417
- 注册: 2012-11-12
- 性别: GG
- 专业: 岩土与基础工程
3楼2013-04-02 22:19:41
枫桥夜过
木虫 (小有名气)
铁杆会员
- 应助: 0 (幼儿园)
- 金币: 2017.3
- 帖子: 187
- 在线: 59.2小时
- 虫号: 2121417
- 注册: 2012-11-12
- 性别: GG
- 专业: 岩土与基础工程
4楼2013-04-02 22:20:07