| 查看: 2992 | 回复: 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 ] |
回复此楼» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
saladin983
铁杆木虫 (正式写手)
- 数学EPI: 4
- 应助: 6 (幼儿园)
- 贵宾: 0.2
- 金币: 9197.9
- 红花: 3
- 帖子: 687
- 在线: 115.7小时
- 虫号: 448168
- 注册: 2007-11-01
- 专业: 计算数学与科学工程计算
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
javeey(金币+2):谢谢给出帮助，很详细。有兴趣担任数学版的专家么？ 2010-05-27 21:06:15
小木虫(金币+0.5):给个红包，谢谢回帖交流
javeey(金币+2):谢谢给出帮助，很详细。有兴趣担任数学版的专家么？ 2010-05-27 21:06:15
|
看了你给的这个例子，问题类型不太一样，这是个一元的非线性回归，而你的问题是多元线性回归。我看了一下Matlab曲线拟合工具箱的文档，没有看到任何关于多元回归的说明，所有的library model都是一元的。所以我认为这个工具箱在你的问题上没有用处。 网上稍微查了一下，regress函数是做多元线性回归最常用的。关于有效度的问题，我翻了这个文档（http://em.hzu.edu.cn/uploadfile/2005113002830244.ppt），现学了一下多元线性回归的内容，里面倒是有几种检验方法。但是你的模型里面有一个很大的差别，也就是前面说到的，没有常数项（注意：在regress说明档里面就提到了），这个跟一般的线性回归模型都不同，所以他们的一些校验参数似乎也用不上。 看起来这个有效度的参数只能自己定义了，我的设想是norm(Xb-y,2)^2/length(y)（变量定义同13楼），看起来还算是合理，但是我不知道这样做是否有什么理论基础。 |
17楼2010-05-27 20:12:59
absoluty
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.1
- 金币: 1606.9
- 散金: 1548
- 红花: 6
- 帖子: 430
- 在线: 192.1小时
- 虫号: 869225
- 注册: 2009-10-12
- 性别: GG
- 专业: 计算机应用技术
2楼2010-05-25 10:22:09
absoluty
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.1
- 金币: 1606.9
- 散金: 1548
- 红花: 6
- 帖子: 430
- 在线: 192.1小时
- 虫号: 869225
- 注册: 2009-10-12
- 性别: GG
- 专业: 计算机应用技术
★ ★
Doctorcbw(金币+2):谢谢 2010-05-25 12:06:15
swenson(金币+1):5 2010-05-25 18:44:19
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
bluesine
铁杆木虫 (职业作家)
科苑小木虫
- 数学EPI: 5
- 应助: 132 (高中生)
- 贵宾: 1.991
- 金币: 9496.3
- 散金: 89
- 红花: 19
- 帖子: 3534
- 在线: 363.8小时
- 虫号: 869544
- 注册: 2009-10-12
- 性别: GG
- 专业: 数学物理
4楼2010-05-25 11:26:13