24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

小指甲123

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 126.3
散金: 20
红花: 1
帖子: 67
在线: 32.8小时
虫号: 2926013
注册: 2014-01-11
专业: 化工热力学和基础数据

[求助] 提示矩阵维数必须一致，不清楚在哪里编程出现了错误，求指教已有1人参与

function ff4
clc
Data=[ 303.15    0.7766    0.02522
313.15    0.7702    0.02548
323.15    0.7605    0.02683
333.15    0.7534    0.02740];
T=Data(:,1);
x11=Data(:,2);
x12=Data(:,3);
[A,fval]=fminsearch(@func3,[1,1],[],T,x11,x12)
x11cal=func1(A,T,x11,x12)
x12cal=func2(A,T,x11,x12)
end

function y=func1(A,T,x11,x12)
A12=A(1);
A21=A(2);
alfa=0.2;
R=8.314;
tao12=A12./R./T;
tao21=A21./R./T;
G12=exp(-alfa*tao12);
G21=exp(-alfa*tao21);
x21=1-x11;
x22=1-x12;
gama11=exp(x21.^2.*(tao21.*(G21./(x11+x21.*G21)).^2+tao12.*G12./(x21+x11.*G12).^2));
gama12=exp(x22.^2.*(tao21.*(G21./(x12+x22.*G21)).^2+tao12.*G12./(x22+x12.*G12).^2));
gama21=exp(x11.^2.*(tao12.*(G12./(x21+x11.*G12)).^2+tao21.*G21./(x11+x21.*G21).^2));
gama22=exp(x12.^2.*(tao12.*(G12./(x22+x12.*G12)).^2+tao21.*G21./(x12+x22.*G21).^2));
y=(gama12.*gama22)/(gama11.*gama22-gama12.*gama21)-(gama12.*gama21)/(gama11.*gama22-gama12.*gama21);
end
function y=func2(A,T,x11,x12)
A12=A(1);
A21=A(2);
alfa=0.2;
R=8.314;
tao12=A12./R./T;
tao21=A21./R./T;
G12=exp(-alfa*tao12);
G21=exp(-alfa*tao21);
x21=1-x11;
x22=1-x12;
gama11=exp(x21.^2.*(tao21.*(G21./(x11+x21.*G21)).^2+tao12.*G12./(x21+x11.*G12).^2));
gama12=exp(x22.^2.*(tao21.*(G21./(x12+x22.*G21)).^2+tao12.*G12./(x22+x12.*G12).^2));
gama21=exp(x11.^2.*(tao12.*(G12./(x21+x11.*G12)).^2+tao21.*G21./(x11+x21.*G21).^2));
gama22=exp(x12.^2.*(tao12.*(G12./(x22+x12.*G12)).^2+tao21.*G21./(x12+x22.*G21).^2));
y=(gama11.*gama22)/(gama11.*gama22-gama12.*gama21)-(gama11.*gama21)/(gama11.*gama22-gama12.*gama21);
end
function y=func3(A,T,x11,x12)
x11cal=func1(A,T,x11,x12);
x12cal=func2(A,T,x11,x12);
sigma=sum((x11-x11cal).^2+(x12-x12cal).^2);
y=sigma;
format long g
end

回复此楼

» 猜你喜欢

22408 调剂材料已经有6人回复
285求调剂已经有12人回复
求调剂已经有6人回复
085600材料与化工301分求调剂院校已经有19人回复
277工科求调剂已经有11人回复
277求调剂数一104分已经有13人回复
304求调剂已经有10人回复
336求调剂，一志愿中科大已经有6人回复
071000生物学，一志愿深圳大学296分，求调剂已经有3人回复
一志愿北京化工085600 310分求调剂已经有19人回复

» 本主题相关价值贴推荐，对您同样有帮助:

矩阵A*B=0 如果已知矩阵B和矩阵A的维数，如何求矩阵A 已经有11人回复
求助一个简单的矩阵行列式问题~~~ 已经有12人回复
Mathematica求30乘30矩阵的行列式为何得不到结果，求指教已经有4人回复

1楼 2014-11-01 22:51:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hytao2012

铁杆木虫 (正式写手)

木头虫子

应助: 53 (初中生)
金币: 6343.8
散金: 115
红花: 11
帖子: 479
在线: 206.6小时
虫号: 2050091
注册: 2012-10-08
性别: GG
专业: 固体力学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
小指甲123: 金币+30, ★★★★★最佳答案 2014-11-02 12:55:13

貌似是因为你的func3中x11和和x11cal的形式不同，不能直接相减。比如一个1*2的向量减去1*3的向量就会出现这种问题

赞一下

回复此楼