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

求解含参方程组

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
求解图里的含参方程组。) t& {; U& `7 ^# L* O7 @
参数:i(虚数)、k、e(自然对数)、L1和L2。
6 V8 T7 c; i% L' o) ?, X* U6 @8 G需要求出a、c1、c2、d1、d2和g含有上面参数的表达式,一共6个方程,6个未知数,可以求解,但是应该怎么才能求出含参的表达式呢?4 e5 P0 K3 ^" j5 }6 o
, ?. j& `) \- P: J
1 Z& f5 B+ e% _# @

该用户从未签到

2#
发表于 2020-7-13 15:16 | 只看该作者
建议先做数学推导,依次找到未知变量间的关系,直到找到一个方程Y(可能是超越的),解出该方程Y即可依次代入其他关系方程,一个一个求出其他未知量。方程Y可以使用solve():
, d3 ~8 e8 o, K/ j+ ?clear;clc;" c8 O# g% p" a+ j% c; A: R4 ~$ H
syms x y/ V! m; u+ D' N0 t+ d
y=cos(x)+x;
" ~( l' @* }% S- A0 Y) u* lx0=vpa(solve(y))%cos(x)+x=0

该用户从未签到

3#
发表于 2020-7-13 15:17 | 只看该作者
参考如下代码片段,全局搜索一组方程的解:
/ V0 q  U* j6 \; [3 N1 Ay(x)=a*exp(-c*(x-b)^2)+d;
4 U# ]8 a  q3 m3 Y% n6 t6 \yp(a,b,c,d,x)=y(x);3 Y9 K4 g6 B  c+ D! e& b" Z
dy(x)=diff(y,x); dyp(a,b,c,d,x)=dy(x);+ m+ E! A- n* B( c1 S6 O- M7 l
func1(a,b,c,d)=y(p1(1))-p1(2);
* \) ~: a6 B% D. ]% R# c8 ^func2(a,b,c,d)=y(p2(1))-p2(2);
- h# v. V, P; F5 `3 e! vfunc3(a,b,c,d)=dy(p1(1))-k1;( C9 |  S3 b8 `/ p# ?0 J& z
func4(a,b,c,d)=dy(p2(1))-k2;
) F1 `# B3 R2 v5 h' M: e$ }funcs=[func1(a,b,c,d);func2(a,b,c,d);func3(a,b,c,d);func4(a,b,c,d)];4 ?, i- m3 p# P; Z3 \' l2 m4 T
vs=[a,b,c,d];
1 A1 E2 P, L$ I" a2 cvs0=[h,L,1.5*L/2,h];) H- D$ q! h( q2 F0 B* P# l
funcHd=matlabFunction(funcs,'vars',{[a,b,c,d]});0 t4 z" e+ H6 o9 _! n3 s6 A
n=4;m=20;
( M! f# N+ y( ~[vss,fval,exitFlag] = GlobalSolve(funcHd,n);

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-8-19 11:51 , Processed in 0.125000 second(s), 26 queries , Gzip On.

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

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

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