24小时热门版块排行榜    

CyRhmU.jpeg
查看: 515  |  回复: 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的虚部部应该为负值的,我自己感觉应该是程序出错,但是又难以发现问题,望各位指点,谢谢。
回复此楼
astudent
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanldcr

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我用你的代码算的是正确的啊
没有问题
人在江湖混,遇事要镇定,让子弹再飞一会
2楼2012-07-23 08:18:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pink01

铁杆木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
内容已删除

» 本帖已获得的红花(最新10朵)

天道酬勤!
3楼2012-07-23 12:18:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sonjinww

金虫 (正式写手)

送鲜花一朵
内容已删除
astudent
4楼2012-07-23 18:02:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sonjinww 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见