Matlab常用信号产生函数演示实例——产生方波、三角波、锯齿波及正弦序列波信号

Matlab常用信号产生函数演示实例——产生方波、三角波、锯齿波及正弦序列波信号

依次产生均匀分布的随机序列、高斯白噪声随机序列、方波信号序列、三角波信号序列、正弦波信号序列，以及信噪比SNR为10dB的加性高斯白噪声正弦信号。

clc; % 清空命令行窗口
clear; %清空工作区
close all;

% Matlab常用信号产生函数演示实例

% 依次产生均匀分布的随机序列、高斯白噪声随机序列、
% 方波信号序列、三角波信号序列、正弦波信号序列，
% 以及信噪比SNR为10dB
% 的加性高斯白噪声正弦信号。

% E2_1_BasicWave.m文件源代码
% 产生方波、三角波及正弦序列波信号
% 定义参数
Ps = 10; % 正弦信号功率为 10dBW
Pn = 1;  % 噪声信号功率为 1dBW
f = 100; % 信号频率为100Hz
Fs = 1000;   % 采样频率为 1kHz
width = 0.5; % 函数sawtooth()的尺度参数为0.5
width2 = 1.0; % 函数sawtooth()的尺度参数为1.0
duty = 50;   % 函数square()的尺度参数为50

%产生信号
t = 0:1/Fs:0.1;
c = 2*pi*f*t;
sq = square(c,duty); %产生方波
tr = sawtooth(c,width); %产生三角波
st = sawtooth(c,width2); %产生锯齿波
si = sin(c); %产生正弦波

%产生随机序列信号
noi = rand(1,length(t));    %产生均匀分布的随机序列
noise = randn(1,length(t)); %产生高斯白噪声序列

%产生带有加性高斯白噪声的正弦信号序列
sin_noise = sqrt(2*Ps)*si + sqrt(Pn)*noise;
sin_noise = sin_noise/max(abs(sin_noise)); %归一化处理

%画图
subplot(421); 
plot(t, noi);
axis([0 0.1 -1.1 1.1]);
xlabel('时间(s)','FontSize',8,'Position',[0.08,-1.3,0]);
ylabel('幅度(v)','FontSize',8);
title('均匀随机信号','FontSize',8);

subplot(422);
plot(t, noise);
axis([0 0.1 -max(abs(noise)) max(abs(noise))]);
xlabel('时间(s)','FontSize',8,'Position',[0.08,-3.2,0]);
ylabel('幅度(v)','FontSize',8);
title('高斯白噪声信号','FontSize',8);

subplot(423);
plot(t, sq);
axis([0 0.1 -1.1 1.1]);
xlabel('时间(s)','FontSize',8,'Position',[0.08,-1.3,0]);
ylabel('幅度(v)','FontSize',8);
title('方波信号','FontSize',8);

subplot(424);
plot(t, tr);
axis([0 0.1 -1.1 1.1]);
xlabel('时间(s)','FontSize',8,'Position',[0.08,-1.3,0]);
ylabel('幅度(v)','FontSize',8);
title('三角波信号','FontSize',8);

subplot(425);
plot(t, si);
axis([0 0.1 -1.1 1.1]);
xlabel('时间(s)','FontSize',8,'Position',[0.08,-1.3,0]);
ylabel('幅度(v)','FontSize',8);
title('正弦波信号','FontSize',8);

subplot(426);
plot(t, sin_noise);
axis([0 0.1 -1.1 1.1]);
xlabel('时间(s)','FontSize',8,'Position',[0.08,-1.3,0]);
ylabel('幅度(v)','FontSize',8);
title('SNR=10dB的正弦波信号','FontSize',8);

subplot(427);
plot(t, st);
axis([0 0.1 -1.1 1.1]);
xlabel('时间(s)','FontSize',8,'Position',[0.08,-1.3,0]);
ylabel('幅度(v)','FontSize',8);
title('锯齿波信号','FontSize',8);