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
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
darova 2021년 3월 17일
@ramya k comment's moved here
we have to write numerically and we should get plot for that
please verify the attached file.. in that electric field equation available.. we have to execcute that equation.

댓글을 달려면 로그인하십시오.

 채택된 답변

darova
darova 2021년 3월 17일

0 개 추천

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
ramya k 2021년 3월 18일
tq for your idea, but i didnt get proper pattern, i got some scattered pattern when wrtten code according to your idea. i have uploaded required model pattern and what i got now. i am very new the matlab so please help me.
darova
darova 2021년 3월 18일
please attach the code
ramya k
ramya k 2021년 3월 19일
please find the code
darova
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)
The code and formulas are also not identical. but you assumed A=1 in your code
Can you collaborate?
ramya k
ramya k 2021년 3월 22일
A means amplitude, it can be anything like constant or circle or parabola. thats y i have taken constant as 1 in first execution.
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
ramya k 2021년 4월 7일
thank you so much for your help

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

도움말 센터File Exchange에서 Mathematics에 대해 자세히 알아보기

질문:

2021년 3월 15일

댓글:

2021년 4월 7일

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by