常用序列的MATLAB代码(一)

常用序列的MATLAB代码(一)

1)单位冲激序列

function[x,n] = impseq(n0,ns,nf)

% ns=序列的起点；nf=序列的终点；n0=序列在n0处有一个单位脉冲。

% x=产生的单位采样序列；n=产生序列的位置信息

n = [ns:nf];

x = [(n-n0)==0];

2)单位阶跃序列

function [x,n] = stepseq(n0,n1,n2)

% Generates x(n) = u(n-n0); n1 <= n,n0 <= n2

% ------------------------------------------

% [x,n] = stepseq(n0,n1,n2)

%

if ((n0 < n1) | (n0 > n2) | (n1 > n2))

error('arguments must satisfy n1 <= n0 <= n2')

end

n = [n1:n2];

%x = [zeros(1,(n0-n1)), ones(1,(n2-n0+1))];

x = [(n-n0) >= 0];

1)加法函数

function [y,n] = sigadd(x1,n1,x2,n2)

% implements y(n) = x1(n)+x2(n)

% -----------------------------

% [y,n] = sigadd(x1,n1,x2,n2)

% y = sum sequence over n, which includes n1 and n2

% x1 = first sequence over n1

% x2 = second sequence over n2 (n2 can be different from n1)

%

n = min(min(n1),min(n2)):max(max(n1),max(n2)); % duration of y(n)

y1 = zeros(1,length(n)); y2 = y1;

% initialization

y1(find((n>=min(n1))&(n<=max(n1))==1))=x1; % x1 with duration of y

y2(find((n>=min(n2))&(n<=max(n2))==1))=x2; % x2 with duration of y

y = y1+y2;

% sequence addition

2)反折函数

function [y,n]=sigfold(x,n)

% implements y(n) = x(-n)

% -----------------------------

% [y,n] = sigfold(x,n)

y=fliplr(x);

n=-fliplr(n);

3)乘法函数

function [y,n] = sigmult(x1,n1,x2,n2)

% implements y(n) = x1(n)*x2(n)

% -----------------------------

% [y,n] = sigmult(x1,n1,x2,n2)

% y = produce sequence over n, which includes n1 and n2

% x1 = first sequence over n1

% x2 = second sequence over n2 (n2 can be different from n1)

%

n = min(min(n1),min(n2)):max(max(n1),max(n2)); % duration of y(n)

y1 = zeros(1,length(n)); y2 = y1;

% initialization

y1(find((n>=min(n1))&(n<=max(n1))==1))=x1; % x1 with duration of y

y2(find((n>=min(n2))&(n<=max(n2))==1))=x2; % x2 with duration of y

y = y1.*y2;

% sequence multiplication

4)移位函数

function [y,n]=sigshift(x,m,n0)

% implements y(n) = x(n-n0)

% -----------------------------

% [y,n] = sigshift(x,m,n0)

n=m+n0;

y=x;

///************************************************************

例1：

代码：

%1.1

clc

clear

close all

n = [0:25];

x=0;

for m=0:10

x=x+(m+1)*(impseq(2*m,0,25)-impseq(2*m+1,0,25));

end

stem(n,x);

xlabel('n');ylabel('x(n)');

%ylim([-2,3]);

运行结果：

例2：

代码：

clc

clear

close all

n = [-10:10];

x1 = n.*n.*(stepseq(-5,-10,10) - stepseq(6,-10,10));

x2=10*impseq(0,-10,10);

x3 = 20*0.5.^n.*(stepseq(4,-10,10) - stepseq(10,-10,10));

x = x1+x2+x3;

stem(n,x);

xlabel('n');ylabel('x(n)');

%ylim([0,50]);

运行结果：

例3：

代码：

clc

clear

close all

n = [0:20];

x =0.9.^n.*(cos(0.2 * pi * n+pi/3));

stem(n,x,'r');

xlabel('n');

ylabel('x(n)');

运行结果：

注：

1、clc：清除命令窗口的内容，对工作环境中的全部变量无任何影响

2、clear：清除工作空间的所有变量

3、close all:关闭所有的Figure窗口

4、stem(n,x,‘r’);函数的第三个入口参数为绘制线的颜色，“r”表示红色

“k”为黑色，默认浅蓝色。

常用序列的MATLAB代码(一)相关教程

基于MATLAB下巴特沃斯IIR数字滤波器的实现

基于MATLAB下巴特沃斯IIR数字滤波器的实现 这里写目录标题 实验目的及要求 实验步骤 实验结果及分析 总结 实验目的及要求 1.学生自己运用MATLAB设计IIR数字低通滤波器。方法不限； 2.实现信号的滤波 参考巴特沃斯滤波器的实验原理《数字信号处理》第五章巴特

Hadoop实例学习(六)Hadoop序列化

Hadoop实例学习(六)Hadoop序列化 目录 序列化相关概念 什么是序列化与反序列化 序列化的目的 Hadoop序列化的特点 常用数据序列化类型 定义序列化接口Bean 序列化实例 解决的问题 Map和Reduce分析 Map阶段 Reduce阶段 编写实例 添加log文件 编写Bean对象 编写M

excel常用操作

excel常用操作 首先我们先下载好演示的excel表格 1.设计规则 2.怎样填对单元格内容 四大规则： 1.一个单元格只对应一个数据属性 2.相同的数据的填写要完全一致 那么我们该如何避免： 在录入前，对表格进行设置，控制录入的内容 若我们要保证输入的数据在0-100

Delaunay三角剖分，MATLAB代码

Delaunay三角剖分，MATLAB代码 clear all;clc;close all;data = load('255440.txt');x = data(:,1);y = data(:,2);z = data(:,3);%load seamount% [x,y,z]=peaks;plot3(x,y,z);%plot(x,y,'.','markersize',10)xlabel('x'), ylabel('y')grid ontri = delaunay(

深度学习常用技巧之改变矩阵大小和归一化

深度学习常用技巧之改变矩阵大小和归一化 在深度学习中我们常用numpy模块中的np.shape和np.reshape来查询矩阵数组的形状或者重新定义矩阵的形状(需保证改变形状前后，数组的元素总个数相等) 当我们用一个算法处理一幅图像时，我们需要将输入的矩阵的形状变为

python及pip中常用命令，经常总结

python及pip中常用命令，经常总结 python基本命令 查看python版本 1.未进入python shell python --version 2.进入python shell，有两种方法 方法一： help() 方法二： import syssys.version pip基本的命令： 查看pip版本 pip --version 安装 sudo easy_insta

字符串操作常用函数三

字符串操作常用函数三 字符串操作常用函数三 函数原型；char *strchr(const char *str,char c) str – 要被检索的 c 字符串 c – 在 str 中要搜索的字符 功能； 在参数str所指向的字符串中搜索第一次出现字符c(一个无符号字符)的位置 返回值； 返回一个指向

C语言求不同的子序列

C语言求不同的子序列 C语言求不同的子序列 代码如下： #include stdio.h#include string.hint c[10];// 存入母串里与子串先等的字符 int flag=0;// C数组的下标 void fun(char a[],char b[]){int i;for (i=0;istrlen(a);i++){if(a[i]==b[flag])// 如果a[i]=子