| 查看: 1682 | 回复: 12 | |||
[求助]
如何使用fsolve求解经过微分计算生成的函数的解？不是微分方程 已有1人参与
|
|||
|
标题可能不清楚，意思就是fsolve求解一个方程我已经会了，但是方程必须是固定已知的，而不能是通过微分计算计算出来的。 比方说： x=fsolve('myfun',[0.5],optimset('Display','off')) function q=myfun(p) alpha1 = p; q =25*exp(5/(2*alpha1)) - 100/(9*exp(5/(2*(alpha1 - 1)))) - 250/(9*exp(5/(2*(alpha1 - 1)))*(alpha1 - 1)) - (125*exp(5/(2*alpha1)))/(2*alpha1) - 125/9; 是可以计算出结果的，而 x=fsolve('myfun',[0.5],optimset('Display','off')) function q=myfun(p) RT=2.5; h1=0.2; h2=0.3; r1=abs(h1)^2; r2=abs(h2)^2; [P1, P3, beta1]=solve ('alpha1*log(1+P1*r1/alpha1) = RT','beta1*log(1+P3*r2/beta1) = RT','alpha1 + beta1 = 1','P1','P3','beta1'); P1 = subs(P1); P3 = subs(P3); y = P1 + P3; q = diff (y); alpha1 = p; 就会报错，因为q是经过微分计算所得，alpha1总是无法将数值解带入。 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有78人回复
-
filecode=后面第一个是大写字母
已经有8人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有13人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用1stOpt解出非线性方程组的解作为初值用Matlab求解，解不出来
已经有4人回复
- 求助如何用matlab求非线性方程组的解
已经有8人回复
- 求大神帮忙解题啊！！！
已经有4人回复
- Matlab心得及学习方法（不断更新）
已经有26人回复
- 关于matlab求解非线性指数方程组出现问题
已经有3人回复
- 重金求助高手解决二元高次方程组MatlAB求解问题
已经有3人回复
- 求助matlab 非线性方程组全部解！
已经有6人回复
- 如何使用matlab求解非线性方程组的所有整数解？
已经有9人回复
- fslove命令如何求解这个复杂的四元一次方程组
已经有0人回复
- fslove命令如何求解这个复杂的四元一次方程组
已经有0人回复
- 【求助】求高手看一下
已经有2人回复
Mr__Right
专家顾问 (著名写手)
-
专家经验: +31 - 应助: 317 (大学生)
- 金币: 14456.3
- 散金: 500
- 红花: 54
- 帖子: 2716
- 在线: 950.6小时
- 虫号: 1972612
- 注册: 2012-09-04
- 性别: GG
- 专业: 应用数学方法
- 管辖: 程序语言
3楼2016-02-12 21:57:54
2楼2016-02-12 15:44:20
4楼2016-02-12 22:10:01
Mr__Right
专家顾问 (著名写手)
-
专家经验: +31 - 应助: 317 (大学生)
- 金币: 14456.3
- 散金: 500
- 红花: 54
- 帖子: 2716
- 在线: 950.6小时
- 虫号: 1972612
- 注册: 2012-09-04
- 性别: GG
- 专业: 应用数学方法
- 管辖: 程序语言
chain4165楼2016-02-12 22:26:22
6楼2016-02-14 11:23:52
7楼2016-02-14 19:08:38
8楼2016-02-16 09:28:19
9楼2016-02-16 14:49:32
10楼2016-02-17 18:34:22
25