DSP基本序列MATLAB代码

代码参考《数字信号处理基础及MATLAB实现》冀振元
因为只用于学习使用所以就不放在这里了，课本pdf可私信我（免费）

1.单位冲激序列

function [x,n] = impseq(n0,n1,n2,figureflag)
%impseq(n0,n1,n2,figureflag)
%产生x(n) = delta(n-n0)；n1 <= n0 <= n2
%n0为冲击位置，n1为序列起点，n2为序列终点
%figureflag为画图标志，1为画，0则不画
if( (n0 < n1) | (n0 > n2) | (n1 > n2) ) error('参数必须满足 n1 <= n0 <= n2');
elseif( figureflag ~= 1 & figureflag ~= 0)error('figureflag 只能为1或者0');
end
n = [n1:n2];
x = [(n-n0) == 0];
if(figureflag == 1)
stem(n,x),axis([n1-2,n2+2,-1,2]),xlabel('序列'),ylabel('序列值'),grid on;
end

2.单位阶跃序列

function [x,n] = stepseq(n0,n1,n2,figureflag)
%stepseq(n0,n1,n2,figureflag)
%产生x(n) = u(n-n0)；n1 <= n0 <= n2
%n0为阶跃位置，n1为序列起点，n2为序列终点
%figureflag为画图标志，1为画，0则不画
if( (n0 < n1) | (n0 > n2) | (n1 > n2) ) error('参数必须满足 n1 <= n0 <= n2');
elseif( figureflag ~= 1 & figureflag ~= 0)error('figureflag 只能为1或者0');
end
n = [n1:n2];
x = [(n-n0) >= 0];if (figureflag == 1)stem(n,x),axis([n1-2,n2+2,-1,2]),xlabel('序列'),ylabel('序列值'),grid on;end

3.矩形序列

function [RN,n] = RN(ns,nf,n1,n2,figureflag)
%RN(ns,nf,n1,n2,figureflag)
%ns = 矩阵始点，nf=矩阵终点，[n1,n2]=给出的坐标范围
%figureflag为画图标志，1为画，0则不画
if( figureflag ~= 1 & figureflag ~= 0)error('figureflag 只能为1或者0');
end
n = n1:n2;
RN = stepseq(ns,n1,n2,0) - stepseq(nf,n1,n2,0); %u(n-ns) - u(n-nf)
if (figureflag == 1)stem(n,RN),axis([n1-2,n2+2,-1,2]),xlabel('序列'),ylabel('序列值'),grid on;
end

这里用了之前定义的函数，可以根据自己定义的名字来更改这个函数。

4.实指数序列

function [x,n] = axp(a,ns,nf)
%axp(a,ns,nf)
%a为底 ns为序列开始，nf为序列结束
if (ns > nf)error('ns要小于等于nf');
end
n = [ns:nf];
x = a.^n;

5.复指数序列

clc;clear;
%产生复指数序列x(n) = exp((a+jw)n)
n0 = -1;n1 = 10;
n = n0:n1;
a = 0.4;
w = 0.6;
x = exp((a+i*w).*n);%i或者j都可以
figure(1)
subplot(211),stem(n,real(x),'filled'),axis([n0-1,n1+1,min(real(x))-1,max(real(x))+1]),title('复指数序列');
ylabel('实部'),grid on;
subplot(212),stem(n,imag(x),'filled'),axis([n0-1,n1+1,min(imag(x))-1,max(imag(x))+1]);
ylabel('虚部');
xlabel('n');
grid on;

6.正弦序列

clc;clear;close all;
n = [1:12];
x = 4*sin(0.3*pi.*n+pi/4) + 7*cos(0.7*pi.*n+pi/5);
stem(n,x),grid on;

DSP基本序列MATLAB代码相关推荐

matlab单位采样序列程序,常用序列的MATLAB代码（一）
常用序列的MATLAB代码(一) 常用序列的MATLAB代码(一) 1)单位冲激序列 function[x,n] = impseq(n0,ns,nf) % ns=序列的起点:nf=序列的终点:n0=序 ...
关于matlab生成M序列的代码
#关于matlab生成M序列的代码参考了相关书籍和网上matlab生产m序列的代码,我发现和书上的定义有些不一样(系统辨识原理及其运用-李彦俊)于是我就自己学了一个m文件来供大家参考.使用起来简单方 ...
基于动态时间规整算法(DTW)的语音识别技术研究-含Matlab代码
⭕⭕ 目录 ⭕⭕ ⏩ 一.引言 ⏩ 二.动态时间规整算法基本原理 ⏩ 三.语音识别实例分析 ⏩ 四.参考文献 ⏩ 五.Matlab代码获取 ⏩ 一.引言在语音识别技术的发展过程中,动态时间规整算法 ...
【数字信号处理】线性常系数差分方程 ( 使用 matlab 求解 “ 线性常系数差分方程 “ 示例 | A 向量分析 | B 向量分析 | 输入序列分析 | matlab 代码 )
文章目录一.使用 matlab 求解 " 线性常系数差分方程 " 示例 1.B 向量元素 : x(n) 参数 2.A 向量元素 : y(n) 参数 3.输入序列 4.matlab ...
【数字信号处理】卷积编程实现 ( Matlab 卷积和多项式乘法 conv 函数 | 使用 matlab 代码求卷积并绘图 )
文章目录一.Matlab 卷积和多项式乘法 conv 函数二.使用 matlab 代码求卷积并绘图一.Matlab 卷积和多项式乘法 conv 函数 Matlab 文档地址 : https:// ...
光流 | 基于光流的实时运动物体检测（MATLAB代码）
===================================================== github:https://github.com/MichaelBeechan CSDN: ...
数学规划模型总结（附MatLab代码）
数学规划模型总结(附MatLab代码) 概述定义一般形式分类线性规划(linear programming) 非线性规划(nonlinear programming) 整数规划(integer ...
【优化求解】基于NSGAII算法求解含约束多目标优化问题matlab代码
1 简介 NSGA(非支配排序遗传算法).NSGAII(带精英策略的非支配排序的遗传算法),都是基于遗传算法的多目标优化算法,都是基于pareto最优解讨论的多目标优化. 2 部分代码 clear ...
5G NR OFDM链路层仿真及Matlab代码实现（1）：LDPC信道编译码之5G Tollbox中相关函数使用介绍
目录 0. 写在前面 1. 信道编码 2. 经过信道(❤) 3. LDPC译码 0. 写在前面这是OFDM链路仿真系列的第一部分. 第一次写博客,望前辈们多多指教!!! 本来想着一口气写完整个过程, ...

DSP基本序列MATLAB代码

1.单位冲激序列

2.单位阶跃序列

3.矩形序列

4.实指数序列

5.复指数序列

6.正弦序列

DSP基本序列MATLAB代码相关推荐

最新文章

热门文章