永发信息网

有个MATLAB程序,但是不知道里面是什么意思,不知道里面到底显示的是什么公式和算法,求解释

答案:2  悬赏:0  手机版
解决时间 2021-03-04 12:13
function [c,m,v]=norm1(C)
%对数据进行标准化处理
[n,s]=size(C);
for i=1:n
for j=1:s
c(i,j)=(C(i,j)-mean(C(:,j)))/sqrt(cov(C(:,j)));
end
end
m=mean(C);
for j=1:s
v(1,j)=sqrt(cov(C(:,j)));
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function w=maxdet(A)
%求矩阵的最大特征值
[v,d]=eig(A);
[n,p]=size(d);
d1=d*ones(p,1);
d2=max(d1);
i=find(d1==d2);
w=v(:,i);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [t,q,w,wh,f0,FF]=fun717(px,py,C)
% px自变量的输入个数
% py因变量的输入个数。
% C输入的自变量和因变量组成的矩阵
% t提取的主成分
% q为回归系数。
% w最大特征值所对应的特征向量。
% wh处理后的特征向量
% f0回归的标准化的方程系数
% FF原始变量的回归方程的系数
c=norm1(C); %norm1为标准化函数
%y=c(:,px+1:px+py); %截取标准化的因变量
E0=c(:,1:px); %截取标准化的自变量
F0=c(:,px+1:px+py); %截取标准化的因变量
A=E0'*F0*F0'*E0; %关联矩阵
w(:,1)=maxdet(A); % maxdet为求取最大特征值的函数
t(:,1)=E0*w(:,1); %第一主成分
p(:,1:px)=(E0'*t(:,1)/(t(:,1)'*t(:,1)))'; % p参
E(:,1:px)=E0-t(:,1)*(E0'*t(:,1)/(t(:,1)'*t(:,1)))'; % 残差矩阵
%进入循环
for i=0:px-2
B(:,px*i+1:px*i+px)=E(:,px*i+1:px*i+px)'*F0*F0'*E(:,px*i+1:px*i+px);
w(:,i+2)=maxdet(B(:,px*i+1:px*i+px));
t(:,i+2)=E(:,px*i+1:px*i+px)*w(:,i+2);
p(:,px*i+px+1:px*i+2*px)=(E(:,px*i+1:px*i+px)'*t(:,i+2)/(t(:,i+2)'*t(:,i+2)))';
E(:,px*i+px+1:px*i+2*px)=E(:,px*i+1:px*i+px)-t(:,i+2)*(E(:,px*i+1:px*i+px)'*t(:,i+2)/(t(:,i+2)'*t(:,i+2)))';
end

for s=1:px
q(:,s)=p(1,px*(s-1)+1:px*s)';
end

[n,d]=size(q);

for h=1:px
iw=eye(d);
for j=1:h-1
iw=iw*(eye(d)-w(:,j)*q(:,j)');
end
wh(:,h)=iw*w(:,h);
end

for j=1:py
zr(j,:)=(regress(F0(:,j),t))';
end

for j=1:px
for i=1:py
w1=wh(:,1:j);
zr1=(zr(i,1:j))';
f0(i,:,j)=(w1*zr1)'; %生成标准化变量的方程的系数矩阵
end
[c,m,v]=norm1(C);
ccxx=ones(py,1)*m(1,1:px);
ccy=(v(1,px+1:px+py))'*ones(1,px);
ccx=ones(py,1)*(v(1,1:px));
ff=ccy.*f0(:,:,j)./ccx;
fff=-(sum((ccy.*ccxx.*f0(:,:,j)./ccx)')-m(1,px+1:px+py))';
FF(:,:,j)=[fff,ff]; %生成原始变量方程的常数项和系数矩阵
end

以上就是代码,每块内容有备注,我是想知道里面具体是怎么个意思,求解答,好的可以加分~
最佳答案
是主成分分析法
全部回答
matlab有个共享程序的网站,去搜搜吧
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
有一首歌是女声的,有点忧伤,里面有一句we m
【蚂蚁有眼睛吗】蚂蚁有眼睛吗?能看到什么样
罗马利奥陶瓷怎么去啊,有知道地址的么
电站锅炉单台给水泵故障怎样处理
刚才骂了一个电话销售的 他又打回来找我 他会
72.56轮毂能装72.6的车上吗
休闲坊在什么地方啊,我要过去处理事情
双店小学二年级老师群二班主任语文老师王巍巍
无线网络安全关键字是什么东西啊 求解决
【绿色和蓝色混合是什么颜色】...的RGB颜色值
堤头村村地址在什么地方,想过去办事
十大佛教菩萨是哪些?
非常老实、实在的人在中国的生存空间大,还是
关牌有什么技巧吗?
看演奏吉他的视频时,看到琴头夹的夹子是干嘛
推荐资讯
白花额这个地址在什么地方,我要处理点事
如东 的历史名人有哪些
蒙妮莎养生馆怎么去啊,有知道地址的么
兴隆街地址在什么地方,想过去办事
【风尾竹】家养风尾竹有什么寓意
患儿男,8个月。因化脓性脑膜炎入院。现患儿
求都市类的小说
solidworks安装后,simulation插件加载后为何
你好,森雅S80是正时皮带还是正时链条谢谢
杭州下沙清雅院到商贸城公交怎么座
出差时西服怎么带
布朗熊身体是什么颜色
正方形一边上任一点到这个正方形两条对角线的
阴历怎么看 ?