| 查看: 136 | 回复: 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 ] |
回复此楼» 猜你喜欢
国社科又开始会评了，不知道这次命运如何
已经有13人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
30