I can't get beam pattern plot of the code. The plot should similar to image attached
2 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
%d=lamda/2
clear all; clc;
A11=2*pi;
f1 = @(x) exp(-1i*pi/2*(sqrt(3)*sin(x)+cos(x))) ; A12 = integral(f1,0,2*pi);
f2 = @(x) exp(-1i*pi/2*(sqrt(3)*sin(x)-cos(x))) ; A13 = integral(f2,0,2*pi);
f3 = @(x) exp(1i*pi/2*(sqrt(3)*sin(x)+cos(x))) ; A21 = integral(f3,0,2*pi);
A22 = 2*pi;
f4 = @(x) exp(1i*pi*cos(x)) ; A23 = integral(f4,0,2*pi);
f5 = @(x) exp(1i*pi/2*(sqrt(3)*sin(x)-cos(x))) ; A31 = integral(f5,0,2*pi);
f6 = @(x) exp(-1i*pi*cos(x)) ; A32 = integral(f6,0,2*pi);
A33 = 2*pi;
f7 = @(x) (exp(-1i*pi/2*sqrt(3)*sin(x))) .* sinc(x-(pi/3)) ; b1 = integral(f7,0,2*pi);
f8 = @(x) (exp(1i*pi/2*cos(x))).* sinc(x-(pi/3)) ; b2 = integral(f8,0,2*pi);
f9 = @(x) (exp(-1i*pi/2*cos(x))).* sinc(x-(pi/3)) ; b3 = integral(f9,0,2*pi);
A = [A11 A12 A13; A21 A22 A23; A31 A32 A33];
B = [b1; b2; b3];
W = inv(A)*B;
W1 = W(1) W2 = W(2) W3 = W(3)
%Beam Pattern
theta = 0:0.01:2*pi; phi=0:0.01:2*pi; y1=exp(1i*pi/2*sqrt(3)*sin(theta)); y2=exp(-1i*pi/2*cos(theta)); y3=exp(1i*pi/2*cos(theta));
E = W1*y1+W2*y2+ W3*y3; F = (cos((pi/2)*cos(phi)))/sin(phi); Radiation= E*F; modeval= abs(Radiation); power = modeval.^2;
2 comentarios
Dyuman Joshi
el 27 de Ag. de 2023
I don't see a command for plotting anything in your code.
What is the reference for the image shown in the comment above?
Respuestas (0)
Ver también
Categorías
Más información sobre Analysis, Benchmarking, and Verification en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!