CMFM

ddm

clc
clear all
n=3
M=[3 2 1 10;0 1 6 20;0 0 4 12]
x(n)=M(n,n+1)./M(n,n)
for i=n-1:-1:1
s=M(i,n+1)
for j=i+1:1:n
s=s-M(i,j)*x(j)
end
x(i)=s/M(i,i)
end

clear all
clc
fi1=[0:5:30]
fi2=atand((sind(fi1)+0.40)./(cosd(fi1)-1.38))
plot(fi1,fi2,'DisplayName','fi1 vs fi2','XDataSource','fi1','YDataSource','fi2');figure(gcf)
cb=((cosd(fi1)-1.38)./(cosd(fi2)))
u21=((sind(fi1).*tand(fi2)+cosd(fi1))./(cb.*(sind(fi2).*tand(fi2)+cosd(fi1))))
plot(fi1,u1,'DisplayName','fi1 vs u1','XDataSource','fi1','YDataSource','u1');figure(gcf)
Vqcb=-(sind(fi1)-(cosd(fi1).*tand(fi2))./((sind(fi2).*tand(fi2))+cosd(fi2)))
plot(fi1,Vqcb,'DisplayName','fi1 vs Vqcb','XDataSource','fi1','YDataSource','Vqcb');figure(gcf)
Vqd=0.7*u1
plot(fi1,Vqd,'DisplayName','fi1 vs Vqd','XDataSource','fi1','YDataSource','Vqd');figure(gcf)
Vqs2=0.35*u21
plot(fi1,Vqs2,'DisplayName','fi1 vs Vqs2','XDataSource','fi1','YDataSource','Vqs2');figure(gcf)
Vqs3=0.4*u21
plot(fi1,Vqs3,'DisplayName','fi1 vs Vqs3','XDataSource','fi1','YDataSource','Vqs3');figure(gcf)