matlab复数方程的根,matlab解超越方程的复数根

刺骨的言语ヽ痛彻心扉 2023-01-13 13:57 169阅读 0赞

本帖最后由 zhaocb 于 2015-7-18 16:22 编辑

各位大神，我需要解一个超越方程，'0=(r^2+(50\*m\*(2\*pi\*k\*C)^2))\*tanh(r\*0.04)+i\*r\*(50-m)\*2\*pi\*k\*C'其中r为未知量。但结果只有0，但实际情况是有复数根，而且复数根才是我要的结果；我也试了fsolve解，但初始值不知道是多少，也得不到正确的答案。我也怀疑方程本身没有复数解，但是有大牛用其他代码可以得到正确的复数解，但我有看不懂那是什么语言。。。求助

Matlab

clear;

clc;

S0=load ('C:\\Users\\Zhang\\Desktop\\20150716z\\S40.txt');

emp=load ('C:\\Users\\Zhang\\Desktop\\20150716z\\emp.txt');

sfilm=0.23025851.\*S0(:,2)+i\*pi.\*(1-S0(:,3)/180);

semp=0.23025851.\*emp(:,2)+i\*pi.\*(1-emp(:,3)/180);

w=2\*pi.\*emp(:,1);

sp=sfilm./semp;

A=50\*(1+sp)./(1-sp);

l0=0.04;

ZL=50;

for a=1:801;

m=A(a);

k=w(a);

C=5.8155E-11;

Lext=0.74332929;

syms r

equ=sprintf('0=(r^2+(50\*m\*(2\*pi\*k\*C)^2))\*tanh(r\*0.04)+i\*r\*(50-m)\*2\*pi\*k\*C');

equ=subs(equ);

r=subs(solve(equ,r))；

end

save C:\\Users\\Zhang\\Desktop\\20150716z\\111.txt ZZ -ascii;

不明代码(求大神告知这是什么语言)

\*\*\*\*\*\* Mathematica script to solve eq (4)\*\*\*\*\*\*\*\*

Clear\[fun,f,zexp,g,gout,z,Cap,g0,goini,nf,nc,k,m\];

zexp=Import\["zexp\_r.dat","Table"\]+\\\[ImaginaryI\]Import\["zexp\_i.dat","Table"\];

f=Import\["frec.dat","Table"\]\[\[All,1\]\];

nf=Dimensions\[f\]\[\[1\]\];

nc=Dimensions\[zexp\]\[\[2\]\];

\\!\\(\\(Cap = 4.6366\*10\\^\\(-11\\);\\)\\)

l=0.015;

g0ini=10+15\\\[ImaginaryI\];

\\!\\(\\(fun\[g\_, f\_, z\_\]:= \\((g\*g + z\*50\*\\((Cap\*2\*Pi\*f)\\)\\^2)\\)\*

Tanh\[g\*l\] + \\ \\\[ImaginaryI\]\\ \\((50 -z)\\)\*Cap\*2\*Pi\*f\*g;\\)\\)

gout=Table\[0,\{k,1,nf\},\{m,1,nc\}\];

Do\[g0=g0ini;

Do\[gout\[\[k,m\]\]=g/.FindRoot\[fun\[g,f\[\[k\]\],zexp\[\[k,m\]\]\]==0,\{g,g0\}\];

g0=gout\[\[k,m\]\],\{k,1,nf\}\],\{m,1,nc\}\]

stmp=OpenWrite\["g\_real.dat",FormatType\\\[Rule\]TableForm\];

SetOptions\["g\_real.dat",PageWidth\\\[Rule\]Infinity\];

Write\[stmp,TableForm\[Re\[gout\],TableSpacing->\{0,1\}\]\];

Close\[stmp\];

stmp=OpenWrite\["g\_imag.dat",FormatType\\\[Rule\]TableForm\];

SetOptions\["g\_imag.dat",PageWidth\\\[Rule\]Infinity\];

Write\[stmp,TableForm\[Im\[gout\],TableSpacing->\{0,1\}\]\];

Close\[stmp\];

发表评论取消回复

表情：

评论列表（有 0 条评论，169人围观）

还没有评论，来说两句吧...

相关阅读

相关 Matlab：解方程的仿真

Matlab：解方程的仿真在科学和工程领域，解方程是一项常见的任务。Matlab是一种功能强大的数值计算软件，通过其丰富的工具箱和编程功能，我们可以轻松地进行方程的仿真和求

小鱼儿/ 2024年03月26日 23:32/ 0 赞/ 56 阅读

相关用Python求方程根

1 问题如何利用python求二元一次方程的根？ ![da1f4375b4ad701b8ae74daec15bb37d.jpeg][] 2 方法通过代码输入

我不是女神ヾ/ 2024年03月22日 22:54/ 0 赞/ 97 阅读

相关使用MATLAB求解方程的根

使用MATLAB求解方程的根 MATLAB是一种功能强大的数值计算软件，它提供了各种函数和工具，用于求解各种数学问题，包括方程的根。在MATLAB中，我们可以使用不同的方法来

向右看齐/ 2023年10月15日 15:49/ 0 赞/ 65 阅读

相关 matlab复数方程的根,matlab解超越方程的复数根

本帖最后由 zhaocb 于 2015-7-18 16:22 编辑各位大神，我需要解一个超越方程，'0=(r^2+(50\m\(2\pi\k\C)^2))\tanh(r\0

刺骨的言语ヽ痛彻心扉/ 2023年01月13日 13:57/ 0 赞/ 170 阅读

相关射频功放学习之复数圆方程

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ub

喜欢ヅ旅行/ 2022年12月16日 13:25/ 0 赞/ 214 阅读

相关二分法求解方程的值 matlab

首先定义一个二分求根的函数： function root=bisect(fun,a,b,eps) n=1+round((log(b-a)-log(ep

柔光的暖阳◎/ 2022年07月26日 00:21/ 0 赞/ 41 阅读

相关初步学习二分法求方程的根

二分法求方程的根 [http://icpc.ahu.edu.cn/OJ/Problem.aspx?id=482][http_icpc.ahu.edu.cn_OJ_Proble

谁借莪１个温暖的怀抱￠/ 2022年07月15日 07:53/ 0 赞/ 196 阅读

相关 matlab非线性数值方程的求解

![Image 1][] ![Image 1][] 对于二分法的matlab代码： myfunction.m function y=myfunction(

川长思鸟来/ 2022年05月29日 01:22/ 0 赞/ 274 阅读

相关二分法求方程的根

二分法的原理其实还是零点定理的应用。首先必须要有一个初试区间，而且区间端点的函数值异号，由零点定理可以知道在次区间内至少存在一个零点。然后取区间的中间值，这样就把一个区间

短命女/ 2022年05月15日 10:24/ 0 赞/ 213 阅读

相关 11.3 二分法求方程的根

![70][] ![70 1][] include<iostream> include<cmath> define EPS 1e-6 usi

短命女/ 2022年05月14日 00:05/ 0 赞/ 186 阅读