应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 程序出现矩阵矢量长度不匹配
131/2返回列表12下一页
查看: 3178  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yuyi2009

铜虫 (小有名气)

[求助] 程序出现矩阵矢量长度不匹配 已有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)
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-08-30 12:24:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuyi2009

铜虫 (小有名气)
2楼2014-08-30 13:44:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
逐行调试; 逐个追踪变量, 并disp ( size ( x ))
赞一下 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
3楼2014-08-30 13:47:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuyi2009

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by cooooldog at 2014-08-30 13:47:01
逐行调试; 逐个追踪变量, 并disp ( size ( x ))

这个，不大会弄啊~
赞一下 回复此楼
4楼2014-08-30 13:57:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuyi2009

铜虫 (小有名气)
就没人看一下么~
赞一下 回复此楼
5楼2014-08-30 19:17:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hytao2012

铁杆木虫 (正式写手)

木头虫子
你A等是向量，表达Cp 的时候，T貌似是一个值。这样写应该不行吧。建议修改一下试试
赞一下 回复此楼
6楼2014-08-30 20:11:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuyi2009

铜虫 (小有名气)
引用回帖:
6楼: Originally posted by hytao2012 at 2014-08-30 20:11:31
你A等是向量，表达Cp 的时候，T貌似是一个值。这样写应该不行吧。建议修改一下试试

怎么修改呢~
回复此楼
7楼2014-08-30 22:41:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hytao2012

铁杆木虫 (正式写手)

木头虫子
引用回帖:
7楼: Originally posted by yuyi2009 at 2014-08-30 22:41:07
怎么修改呢~...

因为不知道你要达到什么目的，所以我也不好说。是不是可以Cp的表达式里不含向量可以吗？
赞一下 回复此楼
8楼2014-08-30 23:00:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuyi2009

铜虫 (小有名气)
引用回帖:
8楼: Originally posted by hytao2012 at 2014-08-30 23:00:50
因为不知道你要达到什么目的，所以我也不好说。是不是可以Cp的表达式里不含向量可以吗？...

不含向量那就要一个数一个数的输入，比较麻烦~
赞一下 回复此楼
9楼2014-08-30 23:56:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hytao2012

铁杆木虫 (正式写手)

木头虫子
引用回帖:
9楼: Originally posted by yuyi2009 at 2014-08-30 23:56:01
不含向量那就要一个数一个数的输入，比较麻烦~...

但是，A和T的长度貌似不匹配啊
赞一下 回复此楼
10楼2014-08-31 09:46:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yuyi2009 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 北科281学硕材料求调剂 +3 tcxiaoxx 2026-03-20 3/150 2026-03-20 19:50 by zhukairuo
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 广西大学家禽遗传育种课题组2026年硕士招生（接收计算机专业调剂） +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +21 rare12345 2026-03-18 21/1050 2026-03-20 14:31 by 无懈可击111
[考研] 一志愿武汉理工材料工程专硕调剂 +7 Doleres 2026-03-19 7/350 2026-03-20 10:39 by EBSD
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 14/700 2026-03-20 09:30 by kkcoco25
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 生物学调剂招人！！！ +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by ！本暗一次！
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 一志愿福大288有机化学，求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 一志愿中海洋材料工程专硕330分求调剂 +7 小材化本科 2026-03-18 7/350 2026-03-19 10:46 by Linda Hu
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
信息提示
关闭
请填处理意见
关闭