| 查看: 323 | 回复: 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 |
回复此楼» 猜你喜欢
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
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.5小时
- 虫号: 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