How can I use "fsolve" for a complex function that is changing in each iteration?

Question

Ozgur 2017-6-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/344759-how-can-i-use-fsolve-for-a-complex-function-that-is-changing-in-each-iteration

回答： Ozgur 2017-6-15

code.m
H.mat

在 MATLAB Online 中打开

Hi everyone,

I have complex nonlinear function as below:

When used with fsolve, I find a solution for each iteration, but the values seem not right. Here is the code:

clear
clc
load('H.mat');
z0=1+1i;
x1=0.08;
x2=0.06;
c=340;
fr=200:2:6398;
k1=2*pi*fr/c;
options=optimset('TolFun',1e-12,'TolX',1e-12);
for m=1:length(fr);
f=@(z)cos(z*x2)/cos(z*x1)-(Hp1(m)*exp(Hp3(m)));
[z(m),fval(m)]=fsolve(f,z0,options);
end

(Please notice that H is written in polar form and j is included in the phase part to have analytical form for the fsolve)

When inserted in to function, the values are not resulting zero. That seems I have something wrong going on with the complex solution to this problem. Any advises will be welcomed.

Cheers,

Palazzo