| 查看: 3266 | 回复: 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) |
回复此楼» 猜你喜欢
26年申博自荐-计算机视觉
已经有4人回复
-
考博
已经有4人回复
-
药化及相关博士的申请
已经有3人回复
-
一篇MDPI论文改变了学习工作和生活
已经有4人回复
-
一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇
已经有4人回复
-
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业
已经有3人回复
-
收到国自然专家邀请后几年才会有本子送过来评
已经有3人回复
11楼2014-08-31 10:07:55
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 (初中生)
- 金币: 6350.3
- 散金: 115
- 红花: 11
- 帖子: 479
- 在线: 206.8小时
- 虫号: 2050091
- 注册: 2012-10-08
- 性别: GG
- 专业: 固体力学
6楼2014-08-30 20:11:31