field radiation pattern generation
显示 更早的评论
please help me in writing code for attached file. i have written code like below but it is taking many hours for execution. so please help me for fast execution of equation in the attached file.
code i have written is
theta=[-90:0.01:90];
syms x
L=5;
A(x)=sym(1-(2*x/L)^2);
A(x)=(1-(2*x/L)^2);
phi(x)=sym(0)
phi(x)=0;
u=sin(theta);
fun=(A(x)*(exp((j*pi*10)*((u*x)+phi(x)))));
E= int(fun,x,-1,1)
Edb=20*log10(E);
plot(u,Edb);
grid on
title("radiation pattern");
help me what modifications i have to do for fast execution of this code
2 个评论
Image Analyst
2021-3-15
Does it have to be symbolic, or can we do it numerically where you define some range for your input?
采纳的回答
darova
2021-3-17
Use trapz
t = -90:90; % i choosed step 1
x = -1:0.1:1);
[T,X] = meshgrid(t,x);
u = sin(T);
A = ...%
fun = A.*exp(...);
E = trapz(x,fun,1);
plot(t,E)
7 个评论
ramya k
2021-3-18
darova
2021-3-18
please attach the code
ramya k
2021-3-19
darova
2021-3-19
Range of X axis is [-1 1] here. BUt it's impossible since it's range of 
(x is integration variable, it disappears after integration)
(x is integration variable, it disappears after integration)
The code and formulas are also not identical. 
but you assumed A=1 in your code
but you assumed A=1 in your codeCan you collaborate?
ramya k
2021-3-22
darova
2021-3-22
i made some correction to your code
clc,clear
t = -90:90;
x = -1:0.1:1
[T,X] = meshgrid(t,x);
u = sin(T);
A = 1;
phi = 0;
fun = A.*exp(1j*pi*10*(u.*X+phi));
E = trapz(X,fun,1);
Edb = 20*log10(E);
plot(t,E)
ramya k
2021-4-7
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Mathematics 的更多信息
Translated by 
