| 查看: 1841 | 回复: 7 | |||
yj112358新虫 (小有名气)
|
[交流]
matlab计算带有复数的函数，最后求复数函数的模，结果里面却有复数标志i， 已有2人参与
|
|
matlab计算带有复数的函数，最后求复数函数的模，结果里面却有复数标志i，输入自变量，得到的结果确实是实数，不知道怎么回事; 代码如下： syms n2 k0 fa1 fa2 k2 real ; syms c v0 fw alpha real; syms sild kmax d kxmax kx delta sild real; syms v real ; syms Rsf Rs RA1 real; nord=1.46;next=1.57;n1=1;c=3*(10^8);alpha=pi/2;fa1=pi/3;n2=1.383;k2=0.0083;v0=2923.443;sild=0;kmax=0.6;d=2.17*10^(-9);fw=20.77174;k0=2*pi*v/1000; nx=next*(sin(sild))^2+nord*(cos(sild))^2; f=0.5*(3*(cos(sild))^2-1); kxmax=kmax*(f*(sin(alpha))^2/2+(1-f)/3); V=2*pi*c*fw; delta=2*pi*c*(v-v0); kx=kxmax*V^2/(4*delta^2+V^2); Nx=nx+kx*i; N2=n2+k2*i; I1=d*(Nx^2-N2^2); I2=k0/N2*sin(fa1)*I1; I3=I2*i; fa2=asin(n1/N2*sin(fa1)); b=-sin(fa1-fa2)+I3; b1=sin(fa1+fa2)-I3; rs=b/b1; rf=sin(fa1-fa2)/sin(fa2+fa1); Rsf=rf*conj(rf); Rs=rs*conj(rs); %Rsf=real(rf)*real(rf)+imag(rf)*imag(rf); %Rs=real(rs)*real(rs)+imag(rs)*imag(rs); RA1=-log10(Rs/Rsf) |
回复此楼» 猜你喜欢
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有200人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有40人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有1人回复
-
基金申请
已经有44人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB带有两重循环程序怎么改为向量话编程呢？
已经有6人回复
- MoS2 Castep计算能带结构的问题
已经有3人回复
- 高分求将我的MATLAB程序改写为FORTRAN程序（1000BB）
已经有13人回复
- MATLAB.Demystified.David.McMahon.（《MATLAB 揭秘》中文版—郑碧波译）
已经有82人回复
- matlab求解积分方程组
已经有2人回复
- Matlab程序求助，版本R2011b,程序运行到quadgk函数总是显示找不到指定模块
已经有2人回复
- matlab 中能不能进行矢量间的点乘和叉乘
已经有3人回复
- Matlab 实现DGGE电泳图谱的 泳道多样性指数，相似性指数，泳道聚类分析
已经有26人回复
- matlab for循环中带有矩阵运算 简化问题
已经有2人回复
- matlab对图像加噪
已经有2人回复
- 【求助】求助一些用matlab做物理方面计算的基础文献
已经有11人回复
- 【原创】matlab_2009b
已经有8人回复
- 【分享】Matlab 2009b win32/64 光盘镜像（4.7GB）
已经有24人回复
- 【分享】推荐一本好书《精通MATLAB最优化计算》
已经有880人回复
2楼2015-09-06 08:25:18
yj112358
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 82.4
- 散金: 16
- 帖子: 96
- 在线: 103.5小时
- 虫号: 2837253
- 注册: 2013-11-29
- 专业: 胶体与界面化学
3楼2015-09-06 08:40:21
yj112358
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 82.4
- 散金: 16
- 帖子: 96
- 在线: 103.5小时
- 虫号: 2837253
- 注册: 2013-11-29
- 专业: 胶体与界面化学
4楼2015-09-06 09:38:08
1314168apple
金虫 (知名作家)
- 应助: 68 (初中生)
- 金币: 677
- 红花: 12
- 帖子: 6872
- 在线: 1462.3小时
- 虫号: 287760
- 注册: 2006-10-21
- 专业: 色谱分析
5楼2015-09-06 15:10:39
yj112358
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 82.4
- 散金: 16
- 帖子: 96
- 在线: 103.5小时
- 虫号: 2837253
- 注册: 2013-11-29
- 专业: 胶体与界面化学
6楼2015-09-06 16:09:31
yj112358
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 82.4
- 散金: 16
- 帖子: 96
- 在线: 103.5小时
- 虫号: 2837253
- 注册: 2013-11-29
- 专业: 胶体与界面化学
7楼2015-09-06 16:10:41
1314168apple
金虫 (知名作家)
- 应助: 68 (初中生)
- 金币: 677
- 红花: 12
- 帖子: 6872
- 在线: 1462.3小时
- 虫号: 287760
- 注册: 2006-10-21
- 专业: 色谱分析
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
2013a版本下： 1、试试用factor命令：factor(expand(Rs/Rsf)).这一道题这个方法最好。 2、再用simplify命令的深度化简选项：simplify(expand(Rs/Rsf),'steps',300)。你的问题可以。提高深度值都不可以的话（比较耗时），只能用方法3。1,2有时可以叠用。 3、可以参考 http://blog.sina.com.cn/s/blog_3ebc3c2901018cnu.html |
8楼2015-09-06 19:20:11
15