24小时热门版块排行榜

>论坛更新日志 (3386)
>学术会议 (282)
>导师招生 (156)
>虫友互识 (156)
>找工作 (94)
>考博 (54)
>硕博家园 (35)
>文献求助 (27)
>休闲灌水 (19)
>论文投稿 (14)
>基金申请 (13)
>博后之家 (10)
>公派出国 (10)
>职场人生 (9)
>海外博后 (9)
>教师之家 (9)

返回列表

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

回复此楼

» 猜你喜欢

售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有3人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有4人回复
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有4人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有3人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有3人回复
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有3人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有3人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有3人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有3人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有3人回复

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

矩阵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 (初中生)
金币: 7297.3
散金: 115
红花: 11
帖子: 610
在线: 209小时
虫号: 2050091
注册: 2012-10-08
性别: GG
专业: 固体力学

【答案】应助回帖

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

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

赞一下

回复此楼