本帖最后由 西皮南子 于 2017-3-6 14:10 编辑

%  论文中的卷积混合模型

% 如果出现错误就在matlab安装路径下D:\Program Files\MATLAB\R2012b\toolbox\symbolic\symbolic\@sym中加入一个函数kroneckerDelta.m

% function f = kroneckerDelta(m,n)

% if m==n

%     f=1;

% else

%     f=0;

% end

%==================================================================================================================

clear all;

close all;

clc;

syms t z;

format long;

s1=sin(0.04*pi*t);

s2=sin(0.5*pi*t);

s3=sin(5*pi*t);

% s1=sin(40*t);

% s2=sin(640*t);

% s3=sin(240*t);

for ii=1:1:500

S11(ii)=sin(0.04*pi*ii);

S12(ii)=sin(0.1*pi*ii);

S13(ii)=sin(0.2*pi*ii);

end

figure

subplot(3,1,1);plot(S11);title('混合信号X');grid on;

subplot(3,1,2);plot(S12);title('混合信号Y');grid on;

subplot(3,1,3);plot(S13);title('混合信号Z');xlabel('Time/ms');grid on;

S1=ztrans(s1);

S2=ztrans(s2);

S3=ztrans(s3);

S=[S1;S2;S3];

A12=0.2026+0.4508*z^(-1)-0.1932*z^(-2)-0.1501*z^(-3)-0.4830*z^(-4);

A13=-0.0527+0.0999*z^(-1)+0.3012*z^(-2)-0.1252*z^(-3)-0.2334*z^(-4);

A21=0.1241+0.4901*z^(-1)+0.0801*z^(-2)-0.2537*z^(-3)-0.2939*z^(-4);

A31=0.3061-0.4048*z^(-1)-0.0692*z^(-2)-0.1318*z^(-3)+0.2656*z^(-4);

A23=-0.4200-0.3447*z^(-1)+0.1322*z^(-2)-0.1208*z^(-3)-0.2313*z^(-4);

A32=-0.0488-0.4316*z^(-1)+0.1648*z^(-2)-0.1555*z^(-3)+0.1081*z^(-4);

AZ=[1,A12,A13;A21,1,A23;A31,A32,1];

X=AZ*S;

x1=X(1,:);

x2=X(2,:);

x3=X(3,:);

y1=iztrans(x1,t);% 得到的是y1(t)函数

y2=iztrans(x2,t);

y3=iztrans(x3,t);

% figure

% subplot(311);ezplot(y1);

% subplot(312);ezplot(y2);

% subplot(313);ezplot(y3);

%将y1(t)，y2(t)，y3(t)离散化并赋值

for i=1:1:2000

y11(i)=double(subs(y1,t,i));%将y1(t)中的t换成i，结果强制转换成double型

y12(i)=double(subs(y2,t,i));

y13(i)=double(subs(y3,t,i));

end

Y=[y11;y12;y13];

figure

subplot(3,1,1);plot(Y(1,:));title('混合信号X');grid on;

subplot(3,1,2);plot(Y(2,:));title('混合信号Y');grid on;

subplot(3,1,3);plot(Y(3,:));title('混合信号Z');xlabel('Time/ms');grid on;

2017-3-6 14:10 上传