| 查看: 635 | 回复: 3 | |||
sonjinww金虫 (正式写手)
|
[交流]
懂MATLAB的进！！！！！ 已有2人参与
|
|
clc; x=0.5*1e9:0.00001*1e9:18*1e9; er_r00=14.991;er_rs=22.921;t=5.096*1e-11;d=0.3; y=er_r00+(er_rs-er_r00)./(1+(i*2*pi*x*t).^(1-d)); plot(x,real(y),'b'); hold on; plot(x,imag(y),'r'); hold on; save_1=[x',real(y)']; save permitivity-real.txt save_1 -ascii; save_2=[x',imag(y)']; save permitivity-imag.txt save_2 -ascii; grid on 如题，y函数是一个对x的复数，分解y为实部和虚部，是本程序的目的。其中d在0-1之间。但是问题出现了，按照程序（如上），虚部竟然是复数！！，这是不可能的，因为根据y 的分解式，见附件虚部等式。当er_rs大于er_r00的条件下，d又在0-1之间，y的虚部部应该为负值的，我自己感觉应该是程序出错，但是又难以发现问题，望各位指点，谢谢。  |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有50人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于动力学参数拟合的matlab程序的一些问题（懂动力学拟合的大神进）
已经有1人回复
- matlab高手请进
已经有16人回复
- 大家看看，我该如何投稿到哪个期刊？《从金庸的武学视角谈计算机编程思维的培养》
已经有10人回复
wanldcr
木虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 2108.7
- 散金: 25
- 红花: 1
- 帖子: 553
- 在线: 274.5小时
- 虫号: 498791
- 注册: 2008-02-04
- 性别: GG
- 专业: 无机非金属类高温超导与磁
2楼2012-07-23 08:18:38
pink01
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 7188.7
- 散金: 20
- 帖子: 729
- 在线: 918.2小时
- 虫号: 891598
- 注册: 2009-11-02
- 性别: GG
- 专业: 能源化工
sonjinww3楼2012-07-23 12:18:38
sonjinww
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2837.2
- 散金: 454
- 红花: 2
- 帖子: 516
- 在线: 250.3小时
- 虫号: 779692
- 注册: 2009-05-25
- 性别: GG
- 专业: 智能材料
4楼2012-07-23 18:02:52
