| 查看: 1772 | 回复: 3 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
wahrrysu木虫 (正式写手)
|
[求助]
求助用matlab或者（1st0pt)编程解方程组
|
|||
|
方程组 x*y=10^(5.17-9810/t+(540.6/t-0.1651)*0.33) ; (0.021-x)/(0.33-y)=92.9064/12.011; 其中x的区间为[0,0.021],y的区间为[0，0.33] 给定一个t值，能够计算出X与y的值 求解当t从800增加到1500时，x,y的解，并做t-x的图 谢谢 |
回复此楼» 猜你喜欢
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有184人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
-
3月 中国石油和化学工业联合会科学技术奖 马上开始，申获奖要点
已经有15人回复
-
中国石油和化学工业联合会科学技术奖的评审标准是什么？
已经有5人回复
-
经费如何使用，花不完怎么办？
已经有5人回复
-
兼职
已经有11人回复
-
海外Gustau Catalan课题组情况了解
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助帮帮我编写一个MATLAB程序求解下面一个方程组！！多谢了
已经有5人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- 有偿求助：matlab怎么解联立方程组？
已经有11人回复
- matlab解非线性方程组解不开
已经有3人回复
- 求助matlab 非线性方程组全部解！
已经有6人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- matlab解一阶常微分方程组
已经有4人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- matlab求解非线性方程组，求助！
已经有6人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- matlab中怎么取方程组的一组解并画图
已经有5人回复
- 求助MATLAB解方程组
已经有6人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab解微分方程组
已经有15人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab解方程组
已经有8人回复
- 【求助】MATLAB 解方程组【已解决】
已经有7人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
wahrrysu
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 4221.3
- 散金: 499
- 红花: 2
- 帖子: 346
- 在线: 206.8小时
- 虫号: 889613
- 注册: 2009-10-31
- 专业: 金属结构材料
3楼2013-07-18 17:28:12
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68183.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.5小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+2, 谢谢 2013-07-18 17:56:04
wahrrysu: 金币+15, ★★★★★最佳答案 2013-07-19 11:49:19
感谢参与，应助指数 +1
csgt0: 金币+2, 谢谢 2013-07-18 17:56:04
wahrrysu: 金币+15, ★★★★★最佳答案 2013-07-19 11:49:19
|
function solve_nonl clear all;clc format long global t t=800; x0=[0.12 0.00001];%y=x(1);x=x(2); n=8; for i=1:n tt(i)=t; xresult=fsolve(@nonfun,x0); t=t+100; u(i)=xresult(1); v(i)=xresult(2); end tt [u' v'] plot(tt,v,'b-o'),axis([750 1550 -0.005 max(v)+0.01]) function f=nonfun(x) global t %y=x(1); %x=x(2); f(1)= x(2)*x(1)-10^(5.17-9810/t+(540.6/t-0.1651)*0.33); f(2)= (0.021-x(2))/(0.33-x(1))-92.9064/12.011; 结果： 0.327285153357280 0.000000363988946 0.327286130401890 0.000007921535971 0.327297142969628 0.000093104953106 0.327375471784637 0.000698986912854 0.327769442659100 0.003746394760137 0.329282927902684 0.015453368820073 0.333947896050736 0.051537408179756 0.345785289192641 0.143100940125555  附图1.jpg |
2楼2013-07-18 17:01:41
wahrrysu
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 4221.3
- 散金: 499
- 红花: 2
- 帖子: 346
- 在线: 206.8小时
- 虫号: 889613
- 注册: 2009-10-31
- 专业: 金属结构材料
4楼2013-07-19 11:49:36
