[交流] 懂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的虚部部应该为负值的，我自己感觉应该是程序出错，但是又难以发现问题，望各位指点，谢谢。

回复此楼