24小时热门版块排行榜

>有奖起名 (21736)
>公派出国 (10143)
>找工作 (9746)
>博后之家 (9539)
>考研 (8924)
>海外博后 (7069)
>留学生活 (6240)
>考博 (5956)
>教师之家 (5471)
>留学DIY (5461)
>职场人生 (4614)
>论坛更新日志 (4447)
>访问学者 (4255)
>虫友互识 (3518)
>文学芳草园 (3510)
>育儿交流 (3106)

返回列表

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

回复此楼

» 猜你喜欢

重要消息，中午系统在维护已经有7人回复
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练已经有5人回复
欢迎发来filecode的Mz6后的代码验证其规律已经有102人回复
时间戳又变了8-15 已经有26人回复
今天维护系统维护祝所有人高中已经有9人回复
哪位老哥知道今年的国自然具体哪一天放榜？已经有16人回复
2027广东省杰青已经有8人回复
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。已经有5人回复
投稿咨询已经有5人回复
时间戳变了，能看出什么问题？已经有6人回复

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

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

【答案】应助回帖

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

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

赞一下

回复此楼