注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mr.Right

不顾一切的去想,于是我们有了梦想。脚踏实地的去做,于是梦想成了现实。

 
 
 

日志

 
 
关于我

人生一年又一年,只要每年都有所积累,有所成长,都有那么一次自己认为满意的花开时刻就好。即使一时不顺,也要敞开胸怀。生命的荣枯并不是简单的重复,一时的得失不是成败的尺度。花开不是荣耀,而是一个美丽的结束,花谢也不是耻辱,而是一个低调的开始。

网易考拉推荐

Matlab program for the voltage on a finite transmission line  

2012-09-17 21:00:59|  分类: 编程 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
% compare numbers for input impedance
 clear
 zg=50;       % generator impedanc
 z0=50;       % characteristic impedance
 zl=100+j*25;  % load impedance
 gg=(zg-z0)/(zg+z0);   % generator reflection coefficient
 gl=(zl-z0)/(zl+z0);   % load reflection coefficient
 % use 300 MHz so wave=1m and lengths will be in wavelengths
 wave=3e8/1e8;
 L=3.3;       % length of line in wavelengths
 b=2*pi/wave; % beta, phase constant
 a=0.2;         % alpha, attenuation constant
 gam=a+j*b;   % gamma
 vg=1;        % generator voltage
 vm=z0*vg/(z0+zg);  
 % enter value of z (z=0 is the input; z=L is the load)
 disp(' ')
 disp(['z=0 is the input; z= ',num2str(L),' is the load'])
 %---------- calculate power into the line (z=0) -------------------
 z=0;
 zp=L-z;
 V=vm*exp(-gam*z)*(1+gl*exp(-2*gam*zp))/(1-gg*gl*exp(-2*gam*L));
 I=vg/(z0+zg)*exp(-gam*z)*(1-gl*exp(-2*gam*zp))/(1-gg*gl*exp(-2*gam*L));
 Pin=.5*real(V*conj(I));
 disp(['power into the line using exact formulas: ',num2str(Pin)])
 
 %---------- calculate power to the load (z=L) ---------------------
 z=L;
 zp=L-z;
 V=vm*exp(-gam*z)*(1+gl*exp(-2*gam*zp))/(1-gg*gl*exp(-2*gam*L));
 I=vg/(z0+zg)*exp(-gam*z)*(1-gl*exp(-2*gam*zp))/(1-gg*gl*exp(-2*gam*L));
 PL=.5*real(V*conj(I));
 disp(['power to the load using exact formulas: ',num2str(PL)])
 
 % closed form
  zin=z0*(zl+z0*tanh(gam*L))/(z0+zl*tanh(gam*L));
  Vin=vg*zin/(zin+z0);
  Iin=Vin/zin;
  Pin=.5*real(Vin*conj(Iin));
  disp(['power to the load using Vin and Iin: ',num2str(Pin)])

% plot V(z)
iz=0;
for z=0:0.05:L
    iz=iz+1;
    Z(iz)=z;
    zp=L-z;
    V(iz)=vm*exp(-gam*z)*(1+gl*exp(-2*gam*zp))/(1-gg*gl*exp(-2*gam*L));
    I=vg/(z0+zg)*exp(-gam*z)*(1-gl*exp(-2*gam*zp))/(1-gg*gl*exp(-2*gam*L));
end
plot(Z,abs(V))
xlabel('z (wavelengths, z=0 is input; z=L is load)')
ylabel('voltage, V (Vg=1 V)')
axis([0,L,0,1])
title(['voltage on line for ZL=',num2str(zl),', line length =',num2str(L)])
 

  评论这张
 
阅读(418)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016