| 查看: 111 | 回复: 0 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】关于分段函数作图问题
|
|||
|
各位高人！ 这几天我做一个分段函数的图，分别画出函数实部和虚部的图，因为涉及到归一化的问题，就有两个矩阵相除的问题，我用了循环语句，但得出图的结果和要求有一定的差距。由于本人初学matlab，所做程序均参照书上拼凑编写，不知正确与否，程序我粘在下面了，方程见附件，请帮忙指正！万分感谢！ 程序：clc; clear all; j=sqrt(-1); w=2*pi*8*10.^9;e0=8.85*10.^(-12);u0=4*pi*10.^(-7); u1=sqrt(e0*u0);u2=sqrt(6*u0*e0);u3=sqrt(6*u0*e0); R=0;T=1;C11=0.5-0.61*j;C12=0.5+0.61*j; C21=1.76+2.15*j;C22=0.14-0.17*j; %参数赋值 r=[]; s=[]; for z0=-0.003:0.00001:0 A=exp(-j*w*u1*z0); h=(sqrt(u0/e0))*A^2; %h为矩阵了 AA=(sqrt(u0/e0))*A^2; s01=real(AA/h); s02=imag(AA/h); r=[r s01]; s=[s s02]; end %第一段函数表达式完 m=[]; n=[]; for z1=0:0.00001:0.003 BB=sqrt((2*u0)/(j*(-3*e0)))*(C11^2*exp(-2*w*u2*z1)-C12^2*exp(2*w*u2*z1)); for z0=-0.003:0.00001:0 h=(sqrt(u0/e0))*exp(-2*j*w*u1*z0); s11=real(BB/h); %归一化 s12=imag(BB/h); %归一化 break end m=[m s11]; n=[n s12]; end %第二段函数表达式完 p=[]; q=[]; for z2=0.003:0.00001:0.006 CC=sqrt((-2*u0)/(j*(3*e0)))*(C21^2*exp(-2*w*u3*(z2-0.003))-C22^2*exp(2*w*u3*(z2-0.003))); for z0=-0.003:0.00001:0 h=(sqrt(u0/e0))*exp(-2*j*w*u1*z0); s21=real(CC/h); %归一化 s22=imag(CC/h); %归一化 break end p=[p s21]; q=[q s22]; end %第三段函数表达式完 x=[]; y=[]; for z3=0.006:0.00001:0.01 DD=sqrt(u0/e0)*exp(-2*j*w*u1*(z3-0.006)); for z0=-0.003:0.00001:0 h=(sqrt(u0/e0))*exp(-2*j*w*u1*z0); s31=real(DD/h); %归一化 s32=imag(DD/h);%归一化 break end x=[x s31]; y=[y s32]; end %第四段函数表达式完 z0=-0.003:0.00001:0; s10=[r]; s20=[s]; z1=0:0.00001:0.003; s11=[m]; s21=[n]; z2=0.003:0.00001:0.006; s12=[p]; s22=[q]; z3=0.006:0.00001:0.01; s13=[x]; s23=[y]; plot(z0,s10,'g-',z0,s20,'r--'); hold on; plot(z1,s11,'g-',z1,s21,'r--'); hold on; plot(z2,s12,'g-',z2,s22,'r--'); hold on; plot(z3,s13,'g-',z3,s23,'r--'); legend('-Re{Hy}','--im{Hy}'); [ Last edited by kuhailangyu on 2008-12-18 at 16:38 ] |
回复此楼» 猜你喜欢
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
-
300求调剂，材料科学英一数二
已经有8人回复
-
招收生物学/细胞生物学调剂
已经有5人回复
-
070305高分子化学与物理 304分求调剂
已经有7人回复
-
289求调剂
已经有13人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有9人回复
-
081200-11408-276学硕求调剂
已经有6人回复
-
调剂求院校招收
已经有5人回复
-
调剂310
已经有8人回复
6