| 查看: 3213 | 回复: 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) |
回复此楼» 猜你喜欢
291分调剂
已经有9人回复
-
调剂求收留
已经有34人回复
-
291 求调剂
已经有38人回复
-
22408 312求调剂
已经有17人回复
-
一志愿华中农业071010，320求调剂
已经有6人回复
-
290调剂生物0860
已经有41人回复
-
291求调剂
已经有3人回复
-
211本科材料化工求调剂
已经有23人回复
-
山东省基金2026
已经有9人回复
-
药学求调剂
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- opencv中如何利用匹配的特征点求解2帧图像间的运动矢量矩阵，新人求助~
已经有10人回复
- 矩阵的本征矢量
已经有6人回复
- 《计算机应用》怎么在文章后面标注矢量，向量或矩阵
已经有5人回复
- 考博中矢量怎么写？？书中是以黑体表示出来的比如矢量矩阵A A是黑体的！
已经有4人回复
4楼2014-08-30 13:57:24
cooooldog
铁杆木虫 (著名写手)
ส็็็
- 应助: 237 (大学生)
- 金币: 6101.9
- 散金: 1114
- 红花: 39
- 帖子: 1380
- 在线: 553.8小时
- 虫号: 506699
- 注册: 2008-02-18
- 专业: 模式识别
3楼2014-08-30 13:47:01
hytao2012
铁杆木虫 (正式写手)
木头虫子
- 应助: 53 (初中生)
- 金币: 6345.8
- 散金: 115
- 红花: 11
- 帖子: 479
- 在线: 206.6小时
- 虫号: 2050091
- 注册: 2012-10-08
- 性别: GG
- 专业: 固体力学
6楼2014-08-30 20:11:31
7楼2014-08-30 22:41:07
5