|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我用matlab编了一个程序,目的很简单,就是先求出矩阵M中的四个矩阵元M11,M12,M21,M22,然后求出反射系数R,自变量是波长ramuda,现在的问题是我的R不随着ramuda改变,程序最后只有一个R值,请大虾帮忙看看,小弟新手,程序可能有点乱!恳请赐教,急!
clc;clear all
global sita n1 n2 a b n
sita=0;n1=1.46;n2=2.35;a=0.00000015;b=0.00000015;n=4;
H=ones(2,91);
ramuda=0.0000001:0.00000001:0.000001;
H(1,:)=ramuda;
R=myjuzhen(H);
disp(R);
figure(1);
plot(R(1,:),R(2,:));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y=myjuzhen(x)
global sita n1 n2 a b n
y=ones(2,91);
y(1,:)=x(1,:);
delta1=(2.*3.142./x(1,:)).*n1.*a;
delta2=(2.*3.142./x(1,:)).*n2.*b;
M1=[cos(delta1) (i/n1)*sin(delta1);i*n1*sin(delta1) cos(delta1)];
M2=[cos(delta2) (i/n2)*sin(delta2);i*n2*sin(delta2) cos(delta2)];
M0=M1*M2;
M=(M0)^n;
y(2,:)=((M(1,1)-M(2,2))^2+(M(1,2)-M(2,1))^2)/((M(1,1)+M(2,2))^2+(M(1,2)+M(2,1))^2);
不胜感激!!!!! |
|