How to convert the attached equations in Matlab?

Question

Sadiq Akbar 2024-3-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2089406-how-to-convert-the-attached-equations-in-matlab

评论： Sadiq Akbar 2024-3-3

equations.png

I have two equations which I have attached in the attachment here. I want to convert them in Matlab code. But how? I tried but in vain. I don't know why it cannot be coded? The required desired vector u is given belwo:

u=[1 2 0.1 0.1 3 4 30 40 50 60];% u=[a1 a2 r1 r2 f1 f2 theta1 theta2 phi1 phi2];
fmax=10;
m=1:5;
n=m;

As you can see the values of ap in given equations have 1st two values inside u namely 1 and 2. Likewise, the values of rp in given equations have 3rd and 4th values inside u. Similarly, the fp in given equations have values 3 and 4 inside u. Likewise, the values of thetap are 30 and 40 inside u and the values of phip are 50 and 60 inside u.

I want to find the values of xo and yo in MATLAB?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2024-3-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2089406-how-to-convert-the-attached-equations-in-matlab#answer_1420021

编辑：Torsten 2024-3-2

在 MATLAB Online 中打开

I'm not sure if you mean m.^2.*f.^2./16.*r or m.^2.*f.^2./(16.*r) in your equations.

I assumed that your equations.png is correct.

u=[1 2 0.1 0.1 3 4 30 40 50 60];
a = u(1:2);
r = u(3:4);
f = u(5:6);
theta = u(7:8);
phi = u(9:10);
fmax=10;
m=(1:5).';
xo = sum(a.*exp(-1i*((pi/fmax).*(-m.*f/2).*sind(theta).*cosd(phi)+m.^2.*f.^2./16.*r).*(1-sind(theta).^2.*cosd(phi).^2)),2)
xo = 
   2.9880 + 0.2678i
   2.9975 + 0.0768i
   2.9297 - 0.5680i
   2.4816 - 1.5650i
   1.3000 - 2.4949i
yo = sum(a.*exp(-1i*((pi/fmax).*(-m.*f/2).*sind(theta).*sind(phi)+m.^2.*f.^2./16.*r).*(1-sind(theta).^2.*sind(phi).^2)),2)
yo = 
   2.9647 + 0.4489i
   2.9507 + 0.5263i
   2.9889 + 0.2358i
   2.9693 - 0.4278i
   2.6463 - 1.4077i

7 个评论
显示 5更早的评论隐藏 5更早的评论

Sadiq Akbar 2024-3-3

编辑：Sadiq Akbar 2024-3-3

在 MATLAB Online 中打开

@Torsten

I did it like this:

clear;clc
u=[1 2 0.1 0.1 3 4 30 40 50 60];b=u;
a1 = u(1:2);
r1 = u(3:4);
f1 = u(5:6);
theta1 = u(7:8);
phi1 = u(9:10);
fmax1=10;
m=(1:5).';
n=m;
% for b
a2 = b(1:2);
r2 = b(3:4);
f2 = b(5:6);
theta2 = b(7:8);
phi2 = b(9:10);
fmax2=10;
m=(1:5).';
n=m;
xo = sum(a1.*exp(-1i*((pi/fmax1).*(-m.*f1/2).*sind(theta1).*cosd(phi1)+m.^2.*f1.^2./16.*r1).*(1-sind(theta1).^2.*cosd(phi1).^2)),2)
yo = sum(a1.*exp(-1i*((pi/fmax1).*(-n.*f1/2).*sind(theta1).*sind(phi1)+n.^2.*f1.^2./16.*r1).*(1-sind(theta1).^2.*sind(phi1).^2)),2)
xe = sum(a2.*exp(-1i*((pi/fmax2).*(-m.*f2/2).*sind(theta2).*cosd(phi2)+m.^2.*f2.^2./16.*r2).*(1-sind(theta2).^2.*cosd(phi2).^2)),2)
ye = sum(a2.*exp(-1i*((pi/fmax2).*(-n.*f2/2).*sind(theta2).*sind(phi2)+n.^2.*f2.^2./16.*r2).*(1-sind(theta2).^2.*sind(phi2).^2)),2)
%%%%%%%%%%%%%%%%%%
% MSE
%%%%%%%%%%%%%%%%%%
%e=norm(xo-xe).^2/(M);
errx=norm(xo-xe).^2/(m);
erry=norm(yo-ye).^2/(n);
err=errx+erry

But in this the final err is a row vector instead of a single zero value. I don't know why?

Torsten 2024-3-3

编辑：Torsten 2024-3-3

在 MATLAB Online 中打开

But in this the final err is a row vector instead of a single zero value. I don't know why?

m and n are 5x1 vectors. Dividing a scalar by a 5x1 vector gives a 1x5 vector:

x = 0;
y = (1:5).';
x/y
ans = 1×5
     0     0     0     0     0

Sadiq Akbar 2024-3-3

@Torsten

Thanks a lot for your guiadance.

请先登录，再进行评论。

How to convert the attached equations in Matlab?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

7 个评论
显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

How to convert the attached equations in Matlab?

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

7 个评论 显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

7 个评论
显示 5更早的评论隐藏 5更早的评论