在MATLAB中求离散信号的编程主要涉及以下几个步骤：
生成离散信号
使用`sin`函数生成正弦波信号，例如：
```matlab
N = 100； % 信号长度
t = 0:N-1； % 时间向量
x = sin（2*pi*t/N）； % 离散正弦波信号
```
使用`stem`函数生成冲击信号，例如：
```matlab
n = -5:5；
x = n == 0；
stem（n, x, 'filled'）；
axis（[-5 5 0 1.1 *max（x）]）；
xlabel（'时间（n）'）；
ylabel（'幅度（n）'）；
```
使用`ones`和`zeros`函数生成阶跃信号和单位脉冲信号，例如：
```matlab
N = 20；
x = ones（1,N）；
x（1:3） = 0；
n = [-3:N-4]；
stem（n,x）；
title（'单位阶跃序列'）；
```
绘制离散信号
使用`plot`函数绘制信号，例如：
```matlab
plot（t,x）； % 绘制离散正弦波信号
xlabel（'Time（s）'）；
ylabel（'Amplitude'）；
title（'Discrete Sine Wave'）；
```
基本操作
对离散信号进行加减乘除等基本操作，例如：
```matlab
y = x + 2； % 将离散正弦波信号加上2
```
卷积运算
使用`conv`函数计算两个离散信号的卷积，例如：
```matlab
h = [1, 2, 1, -2, 1, 0, -4, 0, 3]； % 脉冲响应
x = [1, -2, 3, -4, 3, 2, 1]； % 输入序列
y = conv（h, x）；
n = 0:14；
subplot（2,1,1）；
stem（n,y）；
xlabel（'Time index n'）；
ylabel（'Amplitude'）；
title（'Output Obtained by Convolution'）；
```
离散傅里叶变换（DFT）
手动计算DFT，例如：
```matlab
N = 8； % 信号长度
x = [1, 2, 3, 4, 5, 6, 7, 8]； % 时域信号
X = zeros（1, N）； % 频域信号
for k = 0:N-1
X（k+1） = 0；
for n = 0:N-1
X（k+1） = X（k+1） + x（n+1） * exp（-1j * 2 * pi * k * n / N）；
end
end
disp（'时域信号 x（n）:'）；
disp（x）；
disp（'频域信号 X（k）:'）；
disp（X）；
```
计算带宽
使用`fft`函数计算傅里叶变换，然后根据奈奎斯特抽样定理计算带宽，例如：
```matlab
function bandwidth = calcBandwidth（signal, Fs）
N = length（signal）； % 信号长度
f = （0:（N-1））*（Fs/N）； % 构造频率向量
Y = fft（signal）； % 计算傅里叶变换
P = abs（Y/N）.^2； % 计算功率谱密度
P = P（1:N/2+1）； % 只取正频段
P（2:end-1） = 2*P（2:end-1）； % 对于复数信号要把负频部分翻转复制给正频部分
maxP = max（P）； % 计算功率谱密度的峰值
bw3dBIdx = find（P >= maxP/2, 1, 'last'）； % 找到功率谱密度下降3dB的位置
bw = f（bw3dBIdx）； % 带宽为该位置处的频率值
bandwidth = bw*2； % 由于只取了正频段，因此要将带宽数值乘2
end
```
这些步骤涵盖了