电子产业一站式赋能平台

PCB联盟网

搜索
查看: 2012|回复: 5
收起左侧

一直提示串联的矩阵维度不一致,新手求教!

[复制链接]

578

主题

858

帖子

4821

积分

四级会员

Rank: 4

积分
4821
发表于 2023-1-12 16:51:55 | 显示全部楼层 |阅读模式
一直提示串联的矩阵的维度不一致,求大神帮助

clear;clc;
muX=[235;125;0.414;1.093];sigmaX=[39.98;9.25;0.268;0.256]
sLn=sqrt(log(1+(sigmaX(1)/muX(1))^2));mLn=log(muX(1))-sLn^2/2;   
aEv(1)=sqrt(6)*sigmaX(3)/pi;uEv(1)=-psi(1)*aEv(1)-muX(3);   
aEv(2)=sqrt(6)*sigmaX(4)/pi;uEv(1)=-psi(1)*aEv(2)-muX(4);
muX1=muX;sigmaX1=sigmaX;   
x=muX;normX=eps;
while abs(norm(x)-normX)/normX>1e-6
   normX=norm(x);
   g=x(1)-x(2)-x(3)-x(4);
   gX=[1;-0.9;-1.08;-0.756];
   cdfX=[logncdf(x(1),mLn,sLn);1-evcdf(-   x(3),uEv,aEv);1-evcdf(-x(4),uEv,aEv)];
   pdfX=[lognpdf(x(1),mLn,sLn);evpdf(-x(3),uEv,aEv);evpdf(-x(4),uEv,aEv)];
   nc=norminv(cdfX);
   sigmaX1(1:2:3:4)=normpdf(nc)./pdfX;
   muX1(1:2:3:4)=[x(1:2:3:4)-nc.*sigmaX1(1:2:3:4)];
   gs=gX.*sigmaX1;alphaX=-gs/norm(gs);
   bbeta=(g+gX'*(muX1-x))/norm(gs)
   x=muX1+bbeta*sigmaX1.*alphaX;
end
回复

使用道具 举报

591

主题

1774

帖子

5769

积分

四级会员

Rank: 4

积分
5769
发表于 2023-1-12 16:52:30 | 显示全部楼层
急!求大神帮助!
回复 支持 反对

使用道具 举报

591

主题

1774

帖子

5769

积分

四级会员

Rank: 4

积分
5769
发表于 2023-1-12 16:52:39 | 显示全部楼层
还是没有人帮忙看看么?
回复 支持 反对

使用道具 举报

2851

主题

4614

帖子

2万

积分

官网认证

Rank: 3Rank: 3

积分
25789

活跃会员热心会员推广达人优秀版主

发表于 2023-1-12 16:53:22 | 显示全部楼层
解决了吗
回复 支持 反对

使用道具 举报

591

主题

1774

帖子

5769

积分

四级会员

Rank: 4

积分
5769
发表于 2023-1-12 16:53:47 | 显示全部楼层
从中华人民共和国住房和城乡建设部(http://www.mohurd.gov.cn/index.html)上下载中国城市建设统计年鉴(http://www.mohurd.gov.cn/xytj/tjzljsxytjgb/jstjnj/index.html)现有所有年限中中国县级市以上城市人口和城市建设用地面积(从各年全国城市人口和建设用地(按城市分列)列表中获取)数据,构建中国县级市以上城市人口、城市建设用地面积长时间序列数值列表。
求大佬支招
回复 支持 反对

使用道具 举报

0

主题

829

帖子

622

积分

二级会员

Rank: 2

积分
622
发表于 2023-1-12 16:54:28 | 显示全部楼层
发题目出来
回复 支持 反对

使用道具 举报

发表回复

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

关闭

站长推荐上一条 /1 下一条


联系客服 关注微信 下载APP 返回顶部 返回列表