function [yn,y,W,A,Q,d,B,Devals]=ICA_normal(X,cont)
[m,n]=size(X); %获取输入矩阵的行/列数,行数为观测数据的数目,列数为采样点数
%---------白化/球化------
Cx = cov(X',1); %计算协方差矩阵Cx
[evecs,evals,V] = svd(Cx); %计算Cx的特征值和特征向量
Q=evals^(-1/2)*evecs'; %白化矩阵
z=Q*X; %正交矩阵
Evals=diag(evals);
for i=1:m
ICa(i)=sum(Evals(1:i))/sum(Evals);
if ICa(i)>=cont
d=i;
break;
end;
end;
Devals=diag(Evals(1:d));
Devecs=evecs(:,1:d);
%---------正交矩阵B------
B=[eye(d,d);zeros(m-d,d)];% 初始化列向量B的寄存矩阵,B=[b1 b2 ... bd]
for r=1:d % 迭代求取每一个独立元
b=B(:,r);
i=1;maxIterationsNum=1000;j=1; % 设置最大迭代次数(即对于每个独立分量而言迭代均不超过此次数)
while i<=maxIterationsNum+1
if i==maxIterationsNum % 循环结束处理
fprintf('\n第%d分量在%d次迭代内并不收敛。',r,maxIterationsNum);
break;
end
b0=b;
t=z'*b;
g=tanh(t);
dg=1-tanh(t).^2;
b=z*g/n-mean(dg)*b; % 核心公式,见参考文献[qinsizhao]
%%%%%%%%%%%% 对b正交化
for k=1:r-1
b=b-(b'*B(:,k))*B(:,k);
end
b=b/norm(b);
%%%%%%%%%%% 收敛性判断
if abs(abs(b'*b0)-1)<1e-6 % 如果收敛,则保存b
B(:,r)=b;
break;
end
i=i+1;
end
end
yn=B'*z; % ICA分离的独立元
y=sqrt(Devals)*yn;
W=sqrt(Devals)*B'*Q; % X=W*S
A=evecs*sqrt(evals)*B*Devals^(-1/2);
没有合适的资源?快使用搜索试试~ 我知道了~
故障分析基于matlab ICA故障监测【含Matlab源码 1591期】.zip
共15个文件
m:9个
jpg:4个
dat:2个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
5星 · 超过95%的资源 3 下载量 31 浏览量
2021-12-12
13:29:21
上传
评论 1
收藏 301KB ZIP 举报
温馨提示
CSDN海神之光上传的代码均可运行,亲测可用,直接替换数据即可,适合小白; 1、代码压缩包内容 主函数:GDPSO.m; 调用函数:其他m文件;无需运行 运行结果效果图; 2、代码运行版本 Matlab 2019b;若运行有误,根据提示修改;若不会,可私信博主; 3、运行操作步骤 步骤一:将所有文件放到Matlab的当前文件夹中; 步骤二:双击打开GDPSO.m文件; 步骤三:点击运行,等程序运行完得到结果; 4、仿真咨询 如需其他服务,可私信博主或扫描博主博客文章底部QQ名片; 4.1 CSDN博客或资源的完整代码提供 4.2 期刊或参考文献复现 4.3 Matlab程序定制 4.4 科研合作 功率谱估计: 故障诊断分析: 雷达通信:雷达LFM、MIMO、通信、成像、定位、干扰、检测、信号分析、脉冲压缩 滤波估计:SOC估计 目标定位:WSN定位、滤波跟踪、目标定位 生物电信号:肌电信号EMG、脑电信号EEG、心电信号ECG 通信系统:DOA估计、编码译码、变分模态分解、管道泄漏、滤波器、数字信号处理+传输+分析+去噪、数字信号调制、误码率、信号估计、信号检测信号识别融合
资源推荐
资源详情
资源评论
收起资源包目录
【故障分析】基于matlab ICA故障监测【含Matlab源码 1591期】.zip (15个子文件)
【故障分析】基于matlab ICA故障监测【含Matlab源码 1591】
ICA_modified
ICA_monitor.m 1KB
运行结果4.jpg 24KB
main_n.m 1KB
d10.dat 390KB
运行结果3.jpg 16KB
variable_c.m 187B
d10_te.dat 781KB
运行结果1.jpg 48KB
control_limit.m 199B
ComCon.m 1KB
sort_IC_negentropy.m 293B
sort_IC_cum.m 313B
ICA_normal.m 2KB
sort_IC_VRE.m 293B
运行结果2.jpg 54KB
共 15 条
- 1
资源评论
- m0_749616532023-12-12发现一个宝藏资源,资源有很高的参考价值,赶紧学起来~
- 未来↣2022-05-10用户下载后在一定时间内未进行评价,系统默认好评。
- 766289362024-02-25资源质量不错,和资源描述一致,内容详细,对我很有用。
海神之光
- 粉丝: 3w+
- 资源: 2094
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功