matlab如何查找输入输出信号 matlab2018a版本中接地如何找?
matlab2018a版本中接地如何找?
你找是其它信号的接地。仿真设计(尤其是供电系统的仿真)接地,在:SimPowerSystem——elements——ground也可以收索:ground只希望对你有帮助~~
matlab中sigshift函数用法?
function[你out]sigexpand(d,M)
Nlength(d)
outzeros(M,N)
太out(1,:)d
outreshape(out,1,M*N)
可以使用MATLAB对离散时间信号信号并且乘法运算,除了相加,小位移,交叉相乘包括自由变化等。(TheuseforMATLABforcomputingdiscretesignals,includingtheadd,shift,multiplybecausetransform.)
跪求非负独立成分分析的matlab程序?
您好,这样的:
一、FastICA算法的都差不多步骤:
1.对观测到数据接受中心化,使它的均值为0;
2.对数据进行白化,。
3.你选不需要肯定的分量的个数,设迭代次数
4.你选择一个精灵权矢量(随机的)。
5.令,非线性函数的选定见前文。
6.。
7.令。
8.若果不收敛的话,前往第5步。
9.令,要是,赶往第5步。
二.MATLAB源程序及说明:
%下程序为ICA的调用函数,输入为观察的信号,控制输出为解混后的信号
functionZICA(X)
%-----------去均值---------
[M,T]size(X)%获取键入矩阵的行/列数,行数为观测数据的数目,列数为样本采集点数
cumulativemean(X)%均值
whilei1:M
X(i,:)X(i,:)-cumulative(i)*ones(1,T)
end
%---------白化/球化------
Cxcov(X,1)%换算协方差矩阵Cx
[eigvector,eigvalue]eig(Cx)%可以计算Cx的特征值和特征向量
Weigvalue^(-1/2)*eigvector%白化矩阵
ZW*X%正交矩阵
%----------迭代-------
Maxcount10000%的最迭代次数
Critical0.00001%判断是否需要收敛
mM%需要估记的分量的个数Wrand(m)
whilen1:m
WPW(:,n)%精灵召唤权矢量(不可以)
%YWP*Z
%GY.^3%G为非线性函数,不可行y^3等
%haha3*Y.^2%G的导数
count0
LastWPzeros(m,1)
W(:,n)W(:,n)/norm(W(:,n))
whileabs(WP-LastWP)ampabs(WP LastWP)gtCritical
countcount1%迭代次数
LastWPWP%前段时间迭代的值
%WP1/T*Z*((LastWP*Z).^3)-3*LastWP
fori1:m
WP(i)means(Z(i,:).*(tanh((LastWP)*Z)))-(means(1-(tanh((LastWP))*Z).^2)).*LastWP(i)
end
WPPzeros(m,1)
forj1:n-1
WPPWPP(WP*W(:,j))*W(:,j)
end
WPWP-WPP
WPWP/(norm(WP))
ifcountMaxcount
printf(未不能找到相应的信号)
return
end
end
W(:,n)WPendZW*Z
%以下为主程序,主要注意为遗留下来信号的产生,仔细的观察信号和解混信号的作图
clearallclc
N200n1:N%N为样本采集点数
s12*sin(0.02*pi*n)%正弦函数信号
t1:Ns22*square(100*t,50)%方波信号
alinspace(1,-1,25)s32*[a,a,a,a,a,a,a,a]%锯齿信号
s4rand(1,N)%洗技能噪声
S[s1s2s3s4]%信号分成4*N
Arand(4,4)
XA*S%仔细的观察信号
%源信号波形图
figure(1)subplot(4,1,1)plot(s1)axis([0N-5,5])title(源信号)
subplot(4,1,2)plot(s2)axis([0N-5,5])
subplot(4,1,3)plot(s3)axis([0N-5,5])
subplot(4,1,4)plot(s4)xlabel(Time/xd)
%仔细的观察信号(混合信号)波形图
figure(2)subplot(4,1,1)plot(X(1,:))title(仔细的观察信号(混合信号))
subplot(4,1,2)plot(X(2,:))
subplot(4,1,3)plot(X(3,:))subplot(4,1,4)plot(X(4,:))
ZICA(X)
figure(3)subplot(4,1,1)plot(Z(1,:))title(解混后的信号)
subplot(4,1,2)plot(Z(2,:))
subplot(4,1,3)plot(Z(3,:))
subplot(4,1,4)plot(Z(4,:))xlabel(Time/ms)
版权声明:本文内容由互联网用户自发贡献,本站不承担相关法律责任.如有侵权/违法内容,本站将立刻删除。