| 查看: 356 | 回复: 2 | |||
[交流]
一个参数改变后方程组为什么没有解了？？？ 已有2人参与
|
|
小弟有如下一个matlab程序，12个未知数，12个方程。当W=320000时，运行后有解，且正确。当把W微调成W=320001，运行后就cannot find explicit solution了。请大神们指点一下！谢谢大家！ clear;clc syms A1 B1 C1 D1 phi1 A2 B2 C2 D2 phi2 k Wm theta f=125000.0; W=320000.0; z1=61.0; z2=z1*3; a=9.525; f1=k-2*pi*f/Wm; f2=W-Wm/(1-(phi2-phi1)*Wm/(2*pi*f*(z2-z1))); f3=A1-1+(4/pi)*(C1*cos(2*pi*f*z1/Wm)+D1*sin(2*pi*f*z1/Wm)); f4=B1-(4/pi)*(D1*cos(2*pi*f*z1/Wm)-C1*sin(2*pi*f*z1/Wm)); f5=C1-int(cos((2*pi*f/Wm)*(z1^2+4*a^2*(cos(theta))^2)^0.5)*(sin(theta))^2,0,pi/2); f6=D1-int(sin((2*pi*f/Wm)*(z1^2+4*a^2*(cos(theta))^2)^0.5)*(sin(theta))^2,0,pi/2); f7=phi1-atan(B1/A1); f8=A2-1+(4/pi)*(C2*cos(2*pi*f*z2/Wm)+D2*sin(2*pi*f*z2/Wm)); f9=B2-(4/pi)*(D2*cos(2*pi*f*z2/Wm)-C2*sin(2*pi*f*z2/Wm)); f10=C2-int(cos((2*pi*f/Wm)*(z2^2+4*a^2*(cos(theta))^2)^0.5)*(sin(theta))^2,0,pi/2); f11=D2-int(sin((2*pi*f/Wm)*(z2^2+4*a^2*(cos(theta))^2)^0.5)*(sin(theta))^2,0,pi/2); f12=phi2-atan(B2/A2); S=vpasolve(f1,f2,f3,f4,f5,f6,f7,f8,f9,f10,f11,f12); % s1=S.A1 % s2=S.B1 % s3=S.C1 % s4=S.D1 % s5=S.phi1 % s6=S.A2 % s7=S.B2 % s8=S.C2 % s9=S.D2 % s10=S.phi2 % s11=S.k % s12=S.Wm |
回复此楼» 猜你喜欢
281求调剂（0805）
已经有23人回复
-
265求调剂
已经有9人回复
-
279分求调剂 一志愿211
已经有10人回复
-
招收调剂硕士
已经有12人回复
-
296求调剂
已经有8人回复
-
招收博士1-2人
已经有3人回复
-
085700资源与环境308求调剂
已经有11人回复
-
274求调剂
已经有8人回复
-
290求调剂
已经有4人回复
-
一志愿西南交通 专硕 材料355 本科双非 求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 拟合求解参数值
已经有4人回复
- comsol multiphysics 软件的系数型偏微分方程系数设置？？？
已经有9人回复
- 动力学方程组参数求解问题
已经有9人回复
- 师兄师姐们，帮小弟看看这个 Laplace变换后的二阶常微分方程怎么解啊
已经有23人回复
- PDE偏微分方程系数设置问题
已经有3人回复
- COMSOL的PDE方程如何参数设置菜能实现求解
已经有8人回复
- 如何把一个方程的时域数值解变换到频域数值解？
已经有5人回复
- 用1stopt软件做微分方程组的参数估计
已经有13人回复
- 状态方程变换后的一个矩阵是否可以写成左边部分列满秩矩阵的形式？
已经有3人回复
- 用matlab 求解参数方程组
已经有8人回复
- COMSOL中PDE方程的变量个数的设置问题
已经有10人回复
- matlab中参数改变后怎样求得曲线？？？
已经有6人回复
- matlab如何求解带参数的线性方程组
已经有6人回复
- 求助Matlab求解含已知变量的微分方程组（变量随时间变化的离散值已知）
已经有12人回复
- 如何求解带有固定参数的超定方程组？
已经有21人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 微分方程组参数拟合的问题
已经有11人回复
- 改变世界的五个方程(Five Equations That Changed the World)【转载】
已经有87人回复
- 怎么改变参数类型conjg(z)
已经有7人回复
- 二阶偏微分方程Hankel变换后的结果
已经有3人回复
- 关于Jade拟合R值问题。
已经有14人回复
- lsqnonlin函数拟合微分方程组参数拟合问题
已经有10人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- 请教Matlab高手解多元非线性方程组中的参数
已经有6人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2015-04-01 17:55:24
CelestialCYJ
木虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 2482.6
- 红花: 12
- 帖子: 144
- 在线: 43.3小时
- 虫号: 3277506
- 注册: 2014-06-16
- 专业: 结构工程
3楼2015-04-03 20:51:02
50