24小时热门版块排行榜    

查看: 2591  |  回复: 20
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

swenson

新虫 (初入文坛)

[交流] 【求助】 求MAtlab求解程序!!!已有4人参与

小弟是化学专业的,在验证化学公式时由于数学基础不好,不知道怎么弄。能否用matlab编程算出下面的三个方程组,此外要算出置信度或者线性之类的东西,因为我要验证哪个公式符合化学规律,谢谢虫子们帮忙啊!!!
方程组1
3.98*a+59.74*b+1/9*c+9d=67.81
21.76*a+46.39*b+3/7*c+7/3*d=81.19
51.25*a+28.84*b+c+d=93.49
75.31*a+14.43*b+7/3*c+3/7*d=86.71
93.25*a+6.68*b+9*c+1/9*d=97.45

方程组2
3.98*a+59.74*b+1/9*c=67.81
21.76*a+46.39*b+3/7*c=81.19
51.25*a+28.84*b+c=93.49
75.31*a+14.43*b+7/3*c=86.71
93.25*a+6.68*b+9*c=97.45

方程组3
1/9*a+9b=4.09
3/7*a+7/3*b=13.04
a+b=13.4
7/3*a+3/7*b= -3.03
9*a+1/9*b=-2.48
请一定要算出置信度,因为我要凭此看这三个方程组哪个有有效解。谢谢帮助,不胜感激!!!

[ Last edited by bluesine on 2010-5-25 at 13:16 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saladin983

铁杆木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
javeey(金币+2):谢谢参与交流 2010-05-30 22:46:11
引用回帖:
Originally posted by swenson at 2010-05-30 12:32:55:

具体点好吗?

看这个吧:http://www.mathworks.co.kr/acces ... it/bq_5ka6-1_1.html。这是曲线拟合的工具箱的一个说明档,虽然前面说到这个工具箱跟你的模型有点出入,不过都是回归分析,衡量模型好坏的准则是差不多的,而且似乎背后理论的东西不多了,这些量的定义来自于直观理解。所以你可以考虑对拟合后的残量进行分析,定义几个类似的量,比如我之前定义的那个,然后比较不同模型下得到的数值,应该是可以分析出来差异的,毕竟都是数字,总有大小。至于最小二乘法求解的问题,我们之前讨论过的几种方法得到的应该都是一样的解,随便哪个方法都可以解决问题。

[ Last edited by saladin983 on 2010-5-30 at 16:46 ]
21楼2010-05-30 19:23:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答

absoluty

金虫 (正式写手)


Doctorcbw(金币+1):有的不行,有的求不出精确解 2010-05-25 10:26:32
[a,b,c,d]=solve('等式1','等式2'...,'等式8'),你用这个能求取有效解,应该你试试
2楼2010-05-25 10:22:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

absoluty

金虫 (正式写手)

★ ★
Doctorcbw(金币+2):谢谢 2010-05-25 12:06:15
swenson(金币+1):5 2010-05-25 18:44:19
试试这个求解线性方程组
solve,linsolve
例:
A=[5 0 4 2;1 -1 2 1;4 1 2 0;1 1 1 1];
%矩阵的行之间用分号隔开,元素之间用逗号或空格
B=[3;1;1;0]
X=zeros(4,1);%建立一个4元列向量
X=linsolve(A,B)
diff(fun,var,n):对表达式fun中的变量var求n阶导数。

例如:F=sym('u(x,y)*v(x,y)'); %sym()用来定义一个符号表达式
diff(F); %matlab区分大小写
pretty(ans) %pretty():用习惯书写方式显示变量;ans是答案表达式
3楼2010-05-25 11:06:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫


Doctorcbw(金币+1):谢谢 2010-05-25 12:06:07
第一个方程,5方程4未知数,条件过剩。
第一个方程,5方程3未知数,条件过剩。
第一个方程,5方程2未知数,条件过剩。
板凳要做十年冷文章不发一个字
4楼2010-05-25 11:26:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 中南大学 粉末冶金国家重点实验室 陈超教授课题组拟招收 1~2位博士研究生通知 +6 中南大学-金属材 2024-11-14 9/450 2024-11-15 04:02 by 求带走hh
[论文投稿] chemical science和advanced science哪个好一点 +5 yly150 2024-11-12 5/250 2024-11-15 01:09 by ca0yan9
[硕博家园] 和导师撕破脸了 +15 圆满123456 2024-11-13 19/950 2024-11-14 19:10 by yudaoqian88
[教师之家] 北大教授何怀宏曾如此描述他的同行 +11 zju2000 2024-11-09 11/550 2024-11-14 18:14 by frks
[硕博家园] reject后小感 +5 sioc-sunj 2024-11-14 7/350 2024-11-14 18:09 by m_zyj
[论文投稿] angew三个小修,返修过去已经20多天了没动静,有朋友遇到这种情况吗? +3 迟迟未到场 2024-11-14 3/150 2024-11-14 17:53 by TopEdit
[精细化工] 同一个反应相同的反应条件,是不是反应结果相差不大? +4 青霉素 2024-11-11 4/200 2024-11-14 17:33 by zinedine
[基金申请] 博后面上太难了 +10 不安定因素 2024-11-14 12/600 2024-11-14 16:59 by passion_43
[基金申请] 广东省自然科学基金-面上项目消息有了么? +4 pachang 2024-11-14 5/250 2024-11-14 16:41 by owenyaa
[基金申请] 变性人申请基金咋算性别? +20 leicahfs 2024-11-11 20/1000 2024-11-14 15:00 by wanghuawei
[硕博家园] 大龄已婚想读博如何 +15 米娅阳 2024-11-11 18/900 2024-11-14 14:18 by xiaomi0401
[教师之家] 处在人生职业的分水岭 +4 otani 2024-11-13 4/200 2024-11-14 14:17 by mddzwo
[找工作] 咨询一下江西的高校待遇,人文氛围怎么样? +5 akslis2024 2024-11-09 5/250 2024-11-14 13:53 by 啄木鸟、
[基金申请] 博后面上出结果了 +5 wl87925139 2024-11-14 5/250 2024-11-14 10:26 by kingmax18996
[基金申请] 第76批博士后面上大概什么时间公示 +6 探际者 2024-11-11 7/350 2024-11-14 10:06 by Foxicut88
[论文投稿] 要不要撤稿另投 10+4 wangzhesd 2024-11-09 8/400 2024-11-14 09:09 by LIU_V
[论文投稿] 核心期刊那个阶段最难 +3 公共茅斯 2024-11-11 3/150 2024-11-13 15:27 by 13333105938
[论文投稿] 求JACS和ACIE的最新文章模板 10+4 Jane060390 2024-11-10 4/200 2024-11-13 09:13 by 北京莱茵润色
[教师之家] 评正教授需要两个国家级项目,有人用子课题糊弄,结果在评审前资格公示时被举报拿下了 +20 瞬息宇宙 2024-11-12 29/1450 2024-11-13 08:49 by Quakerbird
[论文投稿] 投稿求助 40+3 liukun11 2024-11-12 4/200 2024-11-12 19:13 by liukun11
信息提示
请填处理意见