I can't get beam pattern plot of the code. The plot should similar to image attached
2 views (last 30 days)
Show older comments
%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 Comments
Dyuman Joshi
on 27 Aug 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?
Answers (0)
See Also
Categories
Find more on Analysis, Benchmarking, and Verification in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!