matlab学习001-简单的矩阵输入及绘制信号曲线

目录

1,熟悉简单的矩阵输入

1.1,创建矩阵

1.2,在命令行调用文件中的变量

1.3,ones函数

1.4,who和whos的使用

2,绘制信号曲线 

2.1,实指数信号 

2.2,频率为50Hz的周期方波信号

2.3,信号卷积 *

2.4,产生从-4π到4π的Sa(t)函数

2.5,阶跃信号u(t)

2.6,g(t)=u(t)-u(t-1)

2.7,冲激信号δ(t)

1)stem函数

2.8,u(n-3)

2.9  冲击序列δ(n)


1,熟悉简单的矩阵输入

1.1,创建矩阵

在命令窗口输入矩阵A=[123;456;789]或A=[1,2,3;4,5,6;7,8,9]观察输出结果。

  

观察结果,可知:

经观察可以发现,在创建矩阵的时候如果不在每列数字中间加上逗号或者是空格,那么每一列都会很紧凑的挨在一起。如果不想要使用逗号的话,可以直接使用空格键,让每列分开。

用回车代替矩阵里面的分号,观察输出结果。

观察结果,可知:

输出依旧能够实现矩阵的样式,且更加方便,不用写分号。唯一不足就是定义矩阵的时候占用的行过多,不利于观察。我还是比较喜欢使用分号来分开每一行。

1.2,在命令行调用文件中的变量

选择Filenew菜单中的M-file,输入B=[9,8,7;6,5,4;3,2,1]

保存为B1.m文件,退出编辑环境。此时在命令窗口中使用B1命令就

可调出B矩阵。

需要注意的是,变量名和脚本名不能同名,不然就会无法识别。

如果B1.m文件不在project工程目录下,就会有如下结果:

这个时候可以直接使用window自带的拖拽功能,直接将文件拖拽到project目录下即可。

1.3,ones函数

输入x=ones(3,4),观察输出结果。

观察结果:

Ones有1的意思,在创建矩阵时,会用1来填充矩阵的所有行和列。

1.4,who和whos的使用

输入who和whos观察结果,了解其作用。

 观察结果并了解到其作用:

who只是简单的列出了同一个工程下所有文件的变量名。

whos相比于who列出的信息更详细,如:该变量的名字,大小,类型和所占的字节。

2,绘制信号曲线 

2.1,实指数信号 

%f(t)=3e-t;
%(t>0(取适当的△t)
dt=0.01;
t=0:dt:5;
y=3*exp(-t);
plot(t,y);

 产生3-e^(-t)的信号

%f(t)=3-e-t;
%(t>0(取适当的△t)
dt=0.01;
t=0:dt:5;
y=3-exp(-t);   % exp表示以常数e为底的指数函数
plot(t,y);

2.2,频率为50Hz的周期方波信号

dt=0.001;    % 采样间隔为0.001s
t=0:dt:1;
y=square(2*pi*50*t);
plot(t,y);
axis([0,0.1,-1.1,1.1]);
title('·½²方波');
xlabel('时间/s');
ylabel('幅度/db');
axis manual;

axis - 设置坐标轴范围和纵横比。此 MATLAB 函数 指定当前坐标区的范围。以包含 4 个、6 个或 8 个元素的向量形式指定范围。

2.3,信号卷积 *

%f(t)=e^(-t)  % 实指数信号
%sin(2лt)(0<t<3)
dt=0.01;
t=0:dt:5;
y=exp(-t).*sin(2*pi*t);  % exp表示以常数e为底的指数函数
plot(t,y);

2.4,产生从-4π到4π的Sa(t)函数

%f(t)=sa(t); 样本函数
%(t>0(取适当的△t)
dt=0.001;
t=-4*pi:dt:4*pi;
y=sin(t)./t;
plot(t,y);

2.5,阶跃信号u(t)

t=-1:0.01:3;
f=heaviside(t);
plot(t,f);
axis([-1 3 -0.1 1.2]);

在Matlab中,heaviside函数用于计算阶跃函数。这个函数以0为分界点,当函数的自变量小于0时返回0,当自变量等于0时返回0.5,当自变量大于0时返回1。例如,heaviside(2)返回1,heaviside(-3)返回0, heaviside(0)返回0.5。 

2.6,g(t)=u(t)-u(t-1)

依旧可以使用heaviside函数来生成u(t-1),如下👇 

t=0:0.01:3;
f1=heaviside(t);  % u(t)
f2=heaviside(t-1);% u(t-1)
g=f1-f2;          % u(t)-u(t-1)
plot(t,g);
axis([0 3 -2 2]);

 如果想要绘制u(t-2)的信号,可以用heaviside(t-2)。

2.7,冲激信号δ(t)

t = (-1:0.0001:1);  % 中间的0.0001是采样间隔,也可以理解为采样周期,采样周期越小,采样频率就会越大,得到的图像就越是我们想要的结果
impulse = t==0;     % 冲击信号只有在t=0时才有取值,其它部分值为0
plot(t,impulse);

1)stem函数

stem函数主要用于绘制离散序列,而plot用于绘制连续信号。如下,我先用stem函数来绘制离散的正弦信号,👇 

x=0:0.1:2*pi;
y=sin(x);
stem(y,"filled");

 对应的效果如下图:

 如果我用plot函数来代替stem函数,得到的就是下图:

2.8,u(n-3)

绘制u(x-3)延迟阶跃信号有两种方法,一种是使用heaviside函数,另外一种就是使用表达式,如下,都可以得到图中的信号图👇

t=0:0.001:6;
y=heaviside(t-3);  % 使用heaviside函数绘制阶跃信号
%y=(t>3);  % 阶跃信号只有在指定初始位置之后才有值,之前的值为0
plot(t,y);

如果觉得y轴太“矮”,想要自定义高度,也可以,我就需要使用到之前的axis坐标轴函数来设置x和y相应的坐标,如下👇

t=0:0.01:6;
%y=heaviside(t-3);  % 使用heaviside函数绘制阶跃信号
y=(t>=3);  % 阶跃信号只有在指定初始位置之后才有值,之前的值为0
plot(t,y);
axis([-1 6 0 2]);

如果想要离散的时间序列,可以将plot函数替换成stem函数,如下👇

n=0:0.001:6;
y=heaviside(n-3);  % 使用heaviside函数绘制阶跃信号
%y=(n>3);  % 阶跃信号只有在指定初始位置之后才有值,之前的值为0
%plot(n,y);
stem(n,y);  % 绘制离散序列

因为在对连续时间信号进行抽样时,我使用的是按照采样周期0.001来进行采样,所以,就会显示出有无数多个单位冲击序列组成的单位阶跃序列,如果放大蓝色区域就可以看到冲击序列,如下图👇

如果不想要那么密集的显示,可以将采样周期删除掉,如下图👇

2.9  冲击序列δ(n)

n=-1:3;
y=(n==0);
stem(n,y);

有问题请在评论区留言或者是私信我,一天8h在线。

相关推荐

最近更新

  1. TCP协议是安全的吗?

    2024-04-09 08:06:04       18 阅读
  2. 阿里云服务器执行yum,一直下载docker-ce-stable失败

    2024-04-09 08:06:04       19 阅读
  3. 【Python教程】压缩PDF文件大小

    2024-04-09 08:06:04       18 阅读
  4. 通过文章id递归查询所有评论(xml)

    2024-04-09 08:06:04       20 阅读

热门阅读

  1. 初探 Spring Boot 源码:揭秘其高效魔法

    2024-04-09 08:06:04       13 阅读
  2. ESP32和STM32的区别

    2024-04-09 08:06:04       13 阅读
  3. PDF锐化

    PDF锐化

    2024-04-09 08:06:04      13 阅读
  4. 光学透雾方法和图像处理算法透雾

    2024-04-09 08:06:04       15 阅读
  5. 修改huggingface 的cache缓存路径

    2024-04-09 08:06:04       17 阅读
  6. django怎么设置把logger.info的日志保存到本地文件

    2024-04-09 08:06:04       14 阅读