arrayMATLAB
R1 = 10000;R2 = 10000;C = 1*10e-6;f = 0:1000;f1 = 1/(2*pi*R1*C);f2 = (R1 + R2)/(2 * pi * R1 * R2 * C);K = 1;%lgspA = logspace(0,5,1000);
Hp = K * (1 + (f./f1))/(1 + (f./f2));plot(log10(f),log10(Hp));
c = 3e8;h = 50e-9;ramd = 633.0d-9;omega = 2.*pi/ramd*c;e2 = 1;enpr = 1.723;e1 = -11.67345024 + 1i*1.1809292;theta = 35*pi/180: 0.001 : 45*pi/180;rpr1 = (cos(theta)./enpr - sqrt((e1-enpr.^2.*sin(theta).^2)./e1)) ./ (cos(theta)./enpr + sqrt(e1-enpr.^2.*sin(theta).^2)./e1);r12 = (sqrt(e1-enpr.^2.*sin(theta).^2)./e1 - sqrt(e2-enpr.^2.*sin(theta).^2)./e2) ./ (sqrt(e1-enpr.^2.*sin(theta).^2)./e1 + sqrt(e2-enpr.^2.*sin(theta).^2)./e2);alpha = 2.*omega./c.*h.*sqrt(e1-enpr.^2.*sin(theta).^2).*1i;R = (rpr1 + r12.*exp(alpha))./(1 + rpr1.*r12.*exp(alpha));t = theta./pi*180;plot(t,abs(R));r1 = [0 0.1 0.2 0.5 1 2];f1 = 0:0.1:1;f2 = 1-f1;for j=1:length(r1) if r1(j)<1 r2(j) = 1-r1(j); else r2(j)=r1(j); end for i=1:11 F1(i) = ((r1(j).*f1(i).^2)+(f1(i).*f2(i)))./((r1(j).*f1(i).^2)+(2*f1(i).*f2(i))+(r2(j).*f2(i).^2)); end plot(f1,F1);hold onend
Best Answer