找回密码
 注册
关于网站域名变更的通知
查看: 497|回复: 3
打印 上一主题 下一主题

求解含参方程组

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-7-13 14:27 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
求解图里的含参方程组。+ `# Y) H1 E) p+ I6 Y. \- y
参数:i(虚数)、k、e(自然对数)、L1和L2。
; k/ D7 A- X* S& R3 ]4 l需要求出a、c1、c2、d1、d2和g含有上面参数的表达式,一共6个方程,6个未知数,可以求解,但是应该怎么才能求出含参的表达式呢?
8 _6 M) W% _6 g9 o. b$ e
5 \- U  m; H7 y3 f' ?5 e ' \9 }9 s9 }" S% C0 Z

该用户从未签到

2#
发表于 2020-7-13 15:16 | 只看该作者
建议先做数学推导,依次找到未知变量间的关系,直到找到一个方程Y(可能是超越的),解出该方程Y即可依次代入其他关系方程,一个一个求出其他未知量。方程Y可以使用solve():
  v: O5 z$ W/ x( N4 G0 R5 ?clear;clc;
. q) n3 h3 e9 _" W. O/ |8 Q. b7 Z0 {syms x y
9 B6 J0 o% B- z$ ^  ny=cos(x)+x;' G. V( `! p" \$ B
x0=vpa(solve(y))%cos(x)+x=0

该用户从未签到

3#
发表于 2020-7-13 15:17 | 只看该作者
参考如下代码片段,全局搜索一组方程的解:
; F# U  u* o) g- Uy(x)=a*exp(-c*(x-b)^2)+d;7 w" Q" Z% j: ^1 a3 \4 N4 B
yp(a,b,c,d,x)=y(x);
8 R6 f# Y9 c6 @+ {8 ndy(x)=diff(y,x); dyp(a,b,c,d,x)=dy(x);' M* o6 ?1 T6 k. P
func1(a,b,c,d)=y(p1(1))-p1(2);
, v* V, [9 z5 M9 r, vfunc2(a,b,c,d)=y(p2(1))-p2(2);
2 J  v7 G) C4 ]4 r. k6 efunc3(a,b,c,d)=dy(p1(1))-k1;, V  p$ K9 h8 F' t( |% v0 @
func4(a,b,c,d)=dy(p2(1))-k2;
  N" A6 Q, Z+ l9 Qfuncs=[func1(a,b,c,d);func2(a,b,c,d);func3(a,b,c,d);func4(a,b,c,d)];
5 k) p( y3 a) \vs=[a,b,c,d];
0 \# P7 M& e. Q; xvs0=[h,L,1.5*L/2,h];$ w5 f) u3 D$ ?& ^$ p/ W) i
funcHd=matlabFunction(funcs,'vars',{[a,b,c,d]});
5 [- U! z3 z0 `3 Mn=4;m=20;
5 F9 A# j- {3 r" U/ o6 \) \, Z3 l" r' r- h[vss,fval,exitFlag] = GlobalSolve(funcHd,n);

该用户从未签到

4#
发表于 2020-7-13 15:17 | 只看该作者
楼上都太厉害了啊
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-24 13:40 , Processed in 0.156250 second(s), 27 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表