资源描述:
《drude模型复数显示》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、1、有关matlab对复数的显示问题: 复数画图当plot(z)中的为复数单变量时,即相当于plot(real(z),imag(z))。即z=a+j*b,那么绘出的图x轴为a,y轴为b,y轴即为j虚轴;Forexample:>>t=0:0.1:2;>>y=2*t.^2+j*3*t;>>plot(y); 此图横坐标为real(y),竖坐标为imag(z)>>2*t.^2 ans= Columns1through11 0 0.0200 0.0800 0.1800 0.3200 0.5000 0.7200
2、 0.9800 1.2800 1.6200 2.0000 Columns12through21 2.4200 2.8800 3.3800 3.9200 4.5000 5.1200 5.7800 6.4800 7.2200 8.0000 >>y y= Columns1through7 0 0.0200+0.3000i 0.0800+0.6000i 0.1800+0.9000i 0.3200+1.2000i 0.5000+1.5000i 0.7200
3、+1.8000i Columns8through14 0.9800+2.1000i 1.2800+2.4000i 1.6200+2.7000i 2.0000+3.0000i 2.4200+3.3000i 2.8800+3.6000i 3.3800+3.9000i Columns15through21 3.9200+4.2000i 4.5000+4.5000i 5.1200+4.8000i 5.7800+5.1000i 6.4800+5.4000i 7.2200+5.7000i 8.0000+6.0000i如果是双变量
4、如plot(t,z),则z中的虚数部分将被丢弃。>>plot(t,y); 此图横坐标为t,纵坐标为y, y=t.^2;即此图t对应y的实部。 要说明复数的运算,先从解以下的二次方程式的复数根谈起上式的根有实部(-2)及虚部(±3),我们就这个复数的表示法来说明MATLAB的复数功能。MATLAB是以i或j字元来代表虚部,其它的复数相关函数有real,imag,conj,abs,angle等等,详见线上说明lookforcomplex。如果复数表示为x=a+bi共轭复数=,复数大小r=,复数向量的夹角θ=tan-1(b/a)复数实部a=rcosθ,复数虚
5、部b=rsinθ,复数指数表示法x=rei上述各函数对应MATLAB的复数指令为a=real(x),b=imag(x),=conj(x),r=abs(x),=angle(x),x=r*exp(i*angle(x))以下是几个复数表示式的例子:>>x=1-2*i;%注意是2*i不是2i>>real(x)%列出实部ans=1>>imag(x)%列出虚部ans=-2>>conj(x)%计算共轭复数ans=1.0000+2.0000i>>abs(x)%计算复数的大小ans=2.2361>>angle(x)%计算复数向量的夹角(以径度表示)ans=-1.1071
6、>>a=1;b=4;c=13;>>x1=(-b+sqrt(b^2-4*a*c))/(2*a)%以解二次方程式根的公式计算复数根x1=-2.0000+3.0000i>>x2=(-b-sqrt(b^2-4*a*c))/(2*a)x2=-2.0000-3.0000i>>y=exp(i)%以复数指数方式表示一个复数y=0.5403+0.8415i>>y=exp(i*pi*0.75)y=-0.7071+0.7071i和复数有关的图以极座标来表示会比一般的卡氏座标要合适,polar指令可以将数据以极座标方式加以绘图,其语法为polar(theta,r),(thet
7、a,r)分别代表极座标上的角度及半径值。以下的例子说明了polar用法:>>t=0:0.01:2*pi;>>r=sin(2*t).*cos(2*t);>>polar(t,r)>>title('Polarplotofsin(2t)cos(2t)')>>angle=0:2*pi/100:2*pi;>>r=angle/(2*pi);>>polar(angle,r)>>title('Polarplot')>>grid