24小时热门版块排行榜

返回列表

小指甲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

回复此楼

» 猜你喜欢

[复试调剂]西南科技大学国防/材料导师推荐已经有6人回复
化学工程321分求调剂已经有12人回复
211本，11408一志愿中科院277分，曾在中科院自动化所实习已经有4人回复
材料专硕326求调剂已经有5人回复
东南大学364求调剂已经有5人回复
国自科面上基金字体已经有7人回复
药学383 求调剂已经有4人回复
286求调剂已经有5人回复
085601求调剂已经有3人回复
302求调剂已经有5人回复

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

矩阵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 (初中生)
金币: 6341.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的向量就会出现这种问题

赞一下

回复此楼