版块导航: 正在加载中...

登录注册

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

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

回复此楼

» 猜你喜欢

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

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

【答案】应助回帖

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

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

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

小指甲123

2楼2014-11-02 09:49:16

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

小指甲123

新虫 (小有名气)

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

送红花一朵

引用回帖:

2楼: Originally posted by hytao2012 at 2014-11-02 09:49:16
貌似是因为你的func3中x11和和x11cal的形式不同，不能直接相减。比如一个1*2的向量减去1*3的向量就会出现这种问题

还是不太清楚，因为理论上x11cal计算都是X11一一对应求出的，不知在哪个地方写的有问题，它内部计算不是按照我们想的那样计算，得到的维数就不一样了呢，很纠结，还麻烦你如果有时间给仔细瞧一眼，谢谢啦

赞一下

回复此楼

3楼2014-11-02 10:57:12

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

hytao2012

铁杆木虫 (正式写手)

木头虫子

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

引用回帖:

3楼: Originally posted by 小指甲123 at 2014-11-02 10:57:12
还是不太清楚，因为理论上x11cal计算都是X11一一对应求出的，不知在哪个地方写的有问题，它内部计算不是按照我们想的那样计算，得到的维数就不一样了呢，很纠结，还麻烦你如果有时间给仔细瞧一眼，谢谢啦...

如果吧fun3改成下面的（当然只是为了找问题，这样得到的结果应该不是你要的），程序就可以运行了，说明的确是你求的x11cal大小有问题。
运行以后从结果里可以看出，x11cal是4*4的。
问题应该就是这个了，因为不知道你程序的目的，看不懂。。。所以不会改。。。
function y=func3(A,T,x11,x12)
x11cal=func1(A,T,x11,x12);
x12cal=func2(A,T,x11,x12);
sigma=sum((x11).^2+(x12).^2);
y=sigma;
format long g
end

赞一下

回复此楼

4楼2014-11-02 11:33:02

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