[英]How can I write in MATLAB?
我解決了NRW(Nicholson-Ross-Weir轉換方法)轉換方法,發現epsilonr(Er)。
NRW轉換方法
clc
clear all
% yansima_genlik= input('Genliği giriniz =')
yansima_genlik= .856;
%disp(['Genlik: ' num2str(yansima_genlik) ' dir. '])
%yansima_faz= input('Fazı giriniz =')
yansima_faz= 163.2;
%disp(['Faz açısı: ' num2str(yansima_faz) ' dir. '])
s11 = yansima_genlik*cosd(yansima_faz)+i*yansima_genlik*sind(yansima_faz);
s22 = s11;
%gecis_genlik= input('Genliği giriniz =')
gecis_genlik= .609;
%disp(['Genlik: ' num2str(gecis_genlik) ' dir. '])
%gecis_faz= input('Fazı giriniz =')
gecis_faz= -140.5;
%disp(['Faz açısı: ' num2str(gecis_faz) ' dir. '])
s21 = gecis_genlik*cosd(gecis_faz)+i*gecis_genlik*sind(gecis_faz);
s12 = s21;
f= 8*10^9;
l=0.4; %örnek uzunluğu
fc=5.26*10^9; %kesim frekansı
lamda0 = 3.75;
lamdac = 5.703;
x = (s11^2-s21^2+1)/(2*s11)
yansima1 = x + sqrt(x^2-1)
yansima2 = x - sqrt(x^2-1)
iletim = (s11+s21-yansima1)/(1-(s11+s21)*yansima1)
a = log(1/iletim)
b = -(((1/(2*pi*l))*a)^2)
v = sqrt(1/b)
p = 1/v
Mr= (1+yansima1)/(v*(1-yansima1)*sqrt((1/lamda0)^2-(1/lamdac)^2))
%Mr=1;
Er= (lamda0^2/Mr)*(((1/lamdac)^2)+b);
Er_1= real(Er)
Er_2=imag(Er)
但是我的問題是這張照片。 我想在MATLAB中編寫“ NIST迭代轉換方法”。 我寫了部分命令。 但是我寫不了更多。 因為我不懂算法。
NIST迭代方法 (某些部分)
%%%NIST Iterative Yöntemi
e0 = 8.85*10^-12;
m0 = 4*pi*10^-7;
b=3;
l1 = 1;
l2 =1;
la = l1+l2+l;
w= 2*pi*f;
isik_hizi = 1/sqrt(e0*m0);
Mr = 1;
m=m0*Mr;
y = i*sqrt((((w^2*Mr*(Er_1+i*Er_2))/isik_hizi^2)-(2*pi/lamdac)^2));
y0 = i* sqrt((w/isik_hizi)^2-(2*pi/lamdac)^2);
yansima = ((y0/m0)-(y/m))/((y0/m0)+(y/m));
T = exp(-y*l);
fx = (s11*s22-s21*s12-(exp(-2*y0)*(la-l))*((T^2)-(yansima^2))/(1-(yansima^2*T^2)))
請幫我。 我想用MATLAB編寫。 但是我不知道命令。
如果要找到兩個根,可以使用fslove ,如下所示:
y = @(Er) i*sqrt((((w^2*Mr*(Er))/isik_hizi^2)-(2*pi/lamdac)^2));
y0 = i* sqrt((w/isik_hizi)^2-(2*pi/lamdac)^2);
yansima = @(Er) ((y0/m0)-(y(Er)/m))/((y0/m0)+(y(Er)/m));
T = @(Er) exp(-y(Er)*l);
fx = @(Er) (s11*s22-s21*s12-(exp(-2*y0)*(la-l))*((T(Er)^2)-(yansima(Er)^2))/(1-(yansima(Er)^2*T(Er)^2)))
options = optimset(optimset('fsolve'), 'TolFun', 1.0e-12, 'TolX',1.0e-12);
Er1 = fsolve(fx, x01, options);
Er2 = fsolve(fx, x02, options);
但就像我不知道NIST方法,我不能幫你找到x01
和x02
,但在fsolve將找到最近的解決方案x01
(或x02)
時, fx(Er1) ~= 0
,其中Er1
是最接近的解決方案x01
。 如果您希望在第一部分中找到更接近Er
的解決方案,則可以使用x01 = real(Er)
。
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.