| 查看: 3208 | 回复: 12 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
程序出现矩阵矢量长度不匹配 已有1人参与
|
||
|
程序如下,曝出以下红字 错误使用 quadl (line 65) 被积函数必须返回长度与输入矢量的长度相同的输出矢量。 出错 therm_calculation (line 17) H_R1 = H_R1_0 + quadl(@(T)Cp(2)+Cp(6)-Cp(1)-Cp(3)-2.*Cp(5),T0,T1) %-------------------------------------------- clear; clc; global T0 T1 A B C D E R T0 = 298.15;T1=453.15;R = 8.314; H_R1_0 = -175.733;H_R2_0 = -41.166; G_R1_0 = -175.733;G_R2_0 = -41.166; A = [4.714 4.396 3.912 3.259 2.8833 4.395]; B = [-6.986 0.628 -3.913 1.35 3.6807 -4.186]; C = [4.211 5.546 1.182 1.502 -0.7720 1.405]; D = [-4.443 -7.024 -1.302 -2.374 0.6915 -1.564]; E = [1.535 2.685 0.515 1.056 -0.2125 0.632]; Cp =@(T)A+B*10^-3.*T+C*10^-5.*T.^2+D*10^-8.*T.^3+E*10^-11.*T.^4; H_R1 = H_R1_0 + quadl(@(T)Cp(2)+Cp(6)-Cp(1)-Cp(3)-2.*Cp(5),T0,T1) H_R2 = H_R2_0 + quadl(@(T)Cp(4)+Cp(6)-Cp(3)-Cp(5),T0,T1) %--------------------------------------------------------------- G_R1 = -G_R1_0/T0/R + quadl(@(T)H_R1*T.^-2/R,T0,T1) G_R2 = -G_R2_0/T0/R + quadl(@(T)H_R2*T.^-2/R,T0,T1) |
回复此楼» 猜你喜欢
290求调剂
已经有22人回复
-
化工学硕294分，求导师收留
已经有32人回复
-
291 求调剂
已经有36人回复
-
食品与营养（0955）271求调剂
已经有16人回复
-
22408 312求调剂
已经有10人回复
-
一志愿鲁东大学071000生物学学硕初试分数276求调剂
已经有26人回复
-
药学求调剂
已经有6人回复
-
材料工程281还有调剂机会吗
已经有39人回复
-
求助调剂，跨调
已经有18人回复
-
08工学 309分求调剂
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- opencv中如何利用匹配的特征点求解2帧图像间的运动矢量矩阵，新人求助~
已经有10人回复
- 矩阵的本征矢量
已经有6人回复
- 《计算机应用》怎么在文章后面标注矢量，向量或矩阵
已经有5人回复
- 考博中矢量怎么写？？书中是以黑体表示出来的比如矢量矩阵A A是黑体的！
已经有4人回复
9楼2014-08-30 23:56:01
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
3楼2014-08-30 13:47:01
4楼2014-08-30 13:57:24
hytao2012
铁杆木虫 (正式写手)
木头虫子
- 应助: 53 (初中生)
- 金币: 6345.8
- 散金: 115
- 红花: 11
- 帖子: 479
- 在线: 206.6小时
- 虫号: 2050091
- 注册: 2012-10-08
- 性别: GG
- 专业: 固体力学
6楼2014-08-30 20:11:31
