24小时热门版块排行榜

返回列表

当前主题已经存档。

愤怒的葡萄8833

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 44.5
帖子: 23
在线:
虫号: 568725
注册: 2008-06-04
专业: 光电工程

[交流] 【求助】关于分段函数作图问题

各位高人！
这几天我做一个分段函数的图，分别画出函数实部和虚部的图，因为涉及到归一化的问题，就有两个矩阵相除的问题，我用了循环语句，但得出图的结果和要求有一定的差距。由于本人初学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 ]

回复此楼

» 猜你喜欢

1楼 2008-07-01 21:04:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主愤怒的葡萄8833 的主题更新

返回列表