利用matlab求解转移矩阵(A矩阵)转换为散射矩阵(S矩阵)

参考《简明微波》书上的公式:

举个例题:

matlab程序如下:

clear all;
clc;
A1=[cos(pi/5) 1i*sin(pi/5);1i*sin(pi/5) cos(pi/5)];%网络1的A矩阵
A2=[2^(1/4) 0;0 (1/2)^(1/4)];%网络2的A矩阵
A3=[0 j;j 0];%网络3的A矩阵
A4=[(1/2)^(1/4) 0;0 2^(1/4)];%网络4的A矩阵
A5=[1/sqrt(2) j/sqrt(2);j/sqrt(2) 1/sqrt(2)];%网络5的A矩阵
A=A1*A2*A3*A4*A5;%整个网络的A矩阵
C=sum(A(:)); %矩阵所有元素求和
B1=sum(A(:,1));   %矩阵第1列元素的和
B2=sum(A(:,2));    %矩阵第2列元素的和
D1=sum(A(1,:));   %矩阵第1行元素的和
D2=sum(A(2,:));   %矩阵第2行元素的和

S11=(1/C)*(D1-D2);%求S11
abs(S11)%求S11模值
angle(S11)*180/pi%求S11相角
S12=(1/C)*2*det(A)%求S12
S21=(1/C)*2%求S21
abs(S21)%求S21模值
S22=(1/C)*(B2-B1);%求S22
abs(S22)%求S22模值
angle(S22)*180/pi%求S22相角

计算结果如下:

|S11|=  0.3333;φ11=-72°;

|S22|=  0.3333;φ22=-90°;

|S21|=  0.9428

答案仅为个人想法,或许有误,欢迎指正!

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值