| 查看: 123 | 回复: 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 ] |
回复此楼» 猜你喜欢
反应很差，大量原料没有反应
已经有8人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有13人回复
-
评审有感
已经有23人回复
-
提交了我也来说说感想
已经有6人回复
-
Sci. Bull. 悲剧经验
已经有7人回复
-
今年审到国自然15份，谈谈感受
已经有20人回复
-
国自然评分
已经有4人回复
-
如果工作了想读博，可以边工作边读全日制嘛？
已经有5人回复
-
如何实现卤原子转化
已经有8人回复
-
河北省自然科学基金
已经有3人回复
30