yaoqi826/matlab_data
07
1[
2 {
3 "Question": "用 MATLAB 命令生成单边指数衰减信号 f(t) = 2e^{-1.5t} u(t)，并绘制出 0 <= t <= 3 的波形图。",
4 "Response": "首先，定义时间范围 t = 0:0.01:3，并利用指数函数 exp(-1.5*t) 计算衰减部分，结合单位阶跃函数 heaviside(t) 来处理 t >= 0 的部分。最后，使用 plot() 绘制波形图，并添加合适的标签和网格以便观察。",
5 "Complex_CoT": "1. 定义时间范围 t，设置为 0 到 3 之间，间隔 0.01；2. 使用 exp() 计算指数部分，乘以 heaviside(t) 以处理单位阶跃函数；3. 使用 plot() 绘制结果；4. 使用 xlabel()、ylabel() 和 title() 设置图形标签；5. 使用 grid on 打开网格以更好地查看波形。"
6 },
7 {
8 "Question": "用 MATLAB 产生正弦信号 f(t) = 2sin(2πt + π/4)，并绘制出 0 <= t <= 3 的波形图。",
9 "Response": "首先，定义时间范围 t = 0:0.01:3，使用 sin() 函数计算正弦信号 f(t) = 2 * sin(2 * pi * t + pi/4)。然后，使用 plot() 绘制信号波形图，最后添加坐标轴标签和标题，启用网格以便观察波形。",
10 "Complex_CoT": "1. 定义时间范围 t，设置为 0 到 3 之间，间隔 0.01；2. 使用 sin() 函数计算正弦信号；3. 使用 plot() 绘制结果；4. 使用 xlabel()、ylabel() 和 title() 设置图形标签；5. 使用 grid on 打开网格以便更好地查看波形。"
11 },
12 {
13 "Question": "用 MATLAB 产生抽样信号 Sa(t)，并绘制出时间为 -6π <= t <= 6π 的波形图。",
14 "Response": "首先，定义时间范围 t = -6*pi:0.1:6*pi，并生成一个连续信号，例如正弦信号 Sa(t) = sin(t)。然后，使用 stem() 绘制离散抽样信号图，若需要绘制连续信号，可使用 plot() 函数。最后，添加坐标轴标签和标题，启用网格以便观察波形。",
15 "Complex_CoT": "1. 定义时间范围 t，从 -6*pi 到 6*pi，步长设置为 0.1；2. 使用 sin() 函数生成正弦信号 Sa(t)；3. 使用 stem() 函数绘制离散的抽样信号图；4. 使用 xlabel()、ylabel() 和 title() 设置图形标签；5. 使用 grid on 打开网格以便观察波形。"
16 },
17 {
18 "Question": "用 MATLAB 画出以 t = 2 为对称中心的短形脉冲信号的波形图。",
19 "Response": "首先，定义时间范围 t = 0:0.01:4，并设置脉冲信号的宽度 w。使用 `heaviside()` 函数生成脉冲信号，计算在 t = 2-w/2 到 t = 2+w/2 之间的信号值。然后，使用 `plot()` 绘制波形图，最后添加坐标轴标签和标题，启用网格以便观察波形。",
20 "Complex_CoT": "1. 定义时间范围 t，设置为 0 到 4，步长为 0.01；2. 使用 `heaviside()` 函数定义脉冲信号；3. 使用 `plot()` 函数绘制脉冲信号图；4. 使用 `xlabel()`, `ylabel()`, `title()` 添加图形标签；5. 使用 `grid on` 打开网格以便更好地观察波形。"
21 },
22 {
23 "Question": "用 MATLAB 产生频率为 10Hz，占空比为 30% 的周期方波信号。",
24 "Response": "首先，确定信号的频率为 10Hz，周期为 0.1秒，占空比为 30%，即高电平持续 0.03秒。定义时间向量 t，设置时间范围为 0 到 1秒，步长为 0.001。然后，使用 `square()` 函数生成周期方波信号。最后，使用 `plot()` 绘制信号波形图，并添加坐标轴标签和标题。",
25 "Complex_CoT": "1. 计算周期 T = 1/10 = 0.1 秒，设置占空比 duty = 30；2. 定义时间向量 t，范围从 0 到 1秒，步长为 0.001；3. 使用 `square(2*pi*f*t, duty)` 生成周期方波信号，f 为频率，duty 为占空比；4. 使用 `plot()` 绘制周期方波图；5. 使用 `xlabel()`、`ylabel()` 和 `title()` 添加标签，使用 `grid on` 打开网格以便更好地观察波形。"
26 },
27 {
28 "Question": "用 MATLAB 产生幅度为 1，宽度为 4，斜率为 -0.5 的非周期三角信号的波形图。",
29 "Response": "首先，定义时间向量 t，设置时间范围，如 -2 到 2，步长为 0.01。然后，使用 `tripuls()` 函数生成三角波信号，设定宽度为 4，幅度为 1。最后，使用 `plot()` 绘制波形图，并通过 `xlabel()`、`ylabel()` 和 `title()` 添加标签，使用 `grid on` 来显示网格。",
30 "Complex_CoT": "1. 定义时间范围 t，设置为 -2 到 2，步长为 0.01；2. 使用 `tripuls(t, width)` 生成三角波信号，设置宽度为 4，幅度为 1；3. 使用 `plot()` 绘制波形图；4. 使用 `xlabel()`、`ylabel()` 和 `title()` 设置标签；5. 使用 `grid on` 打开网格以便更好地观察波形。"
31 },
32 {
33 "Question": "用 MATLAB 产生峰值为 ±1，周期为 2 的周期三角信号波形图。",
34 "Response": "首先，定义时间向量 t，设置时间范围从 0 到 10，步长为 0.01。然后，使用 `sawtooth()` 函数生成三角波信号，调整周期为 2，幅度为 1。最后，使用 `plot()` 绘制波形图，并通过 `xlabel()`、`ylabel()` 和 `title()` 添加坐标轴标签和图形标题，使用 `grid on` 来显示网格。",
35 "Complex_CoT": "1. 定义时间范围 t，设置为 0 到 10，步长为 0.01；2. 使用 `sawtooth(2*pi*f*t, 0.5)` 生成周期性三角波，设定周期为 2，峰值为 ±1；3. 使用 `plot()` 绘制周期三角波图；4. 使用 `xlabel()`、`ylabel()` 和 `title()` 添加标签；5. 使用 `grid on` 打开网格以便观察波形。"
36 },
37 {
38 "Question": "用 MATLAB 命令画出指数序列 x(n) = (-0.6)^n 的波形图。",
39 "Response": "首先，定义时间索引向量 n，例如 n = 0:20，表示生成从 n = 0 到 n = 20 的序列。然后，使用 `x = (-0.6).^n` 计算对应的指数序列。最后，使用 `stem(n, x)` 绘制波形图，并通过 `xlabel()`、`ylabel()` 和 `title()` 添加坐标轴标签和标题，使用 `grid on` 显示网格。",
40 "Complex_CoT": "1. 定义时间向量 n，设置为 0 到 20；2. 使用 `x = (-0.6).^n` 计算指数序列；3. 使用 `stem(n, x)` 绘制离散波形图；4. 使用 `xlabel()`、`ylabel()` 和 `title()` 添加标签；5. 使用 `grid on` 打开网格以便更好地观察波形。"
41 },
42 {
43 "Question": "用 MATLAB 绘制正弦序列 x(n) = sin(nπ/6) 的波形图。",
44 "Response": "首先，定义时间索引向量 n，例如 n = 0:30，表示从 n = 0 到 n = 30 的序列。然后，使用 `x = sin(n*pi/6)` 计算对应的正弦序列。最后，使用 `stem(n, x)` 绘制波形图，并通过 `xlabel()`、`ylabel()` 和 `title()` 添加坐标轴标签和标题，使用 `grid on` 显示网格。",
45 "Complex_CoT": "1. 定义时间向量 n，设置为 0 到 30；2. 使用 `x = sin(n*pi/6)` 计算正弦序列；3. 使用 `stem(n, x)` 绘制离散波形图；4. 使用 `xlabel()`、`ylabel()` 和 `title()` 添加标签；5. 使用 `grid on` 打开网格以便更好地观察波形。"
46 },
47 {
48 "Question": "试用 MATLAB 命令画出下列信号的波形图：\n(1) y1(t) = e^{-2|t|}\n(2) y2(n) = (0.9)^n, (-10 ≤ n ≤ 10)\n(3) y3(n) = e^{jπn}, (-10 ≤ n ≤ 10)",
49 "Response": "首先，为每个信号定义时间向量。对于连续信号 y1(t)，使用 `t = -5:0.1:5`，然后计算 `y1 = exp(-2*abs(t))`，使用 `plot(t, y1)` 绘制波形。对于离散信号 y2(n)，设置 `n = -10:10`，通过 `y2 = (0.9).^n` 计算信号值，使用 `stem(n, y2)` 绘制离散图。对于复数信号 y3(n)，设置 `n = -10:10`，通过 `y3 = exp(1i*pi*n)` 计算复信号，使用 `real()` 和 `imag()` 绘制实部和虚部。",
50 "Complex_CoT": "1. 定义时间向量 t 为 -5 到 5，使用 `exp(-2*abs(t))` 生成信号；2. 使用 `plot(t, y1)` 绘制连续波形；3. 定义离散时间向量 n 为 -10 到 10，使用 `(0.9).^n` 生成离散信号 y2，使用 `stem(n, y2)` 绘制波形；4. 计算复数信号 y3，使用 `exp(1i*pi*n)`，通过 `real()` 和 `imag()` 分别提取实部和虚部；5. 使用 `xlabel()`、`ylabel()`、`title()` 添加标签，使用 `grid on` 打开网格。"
51 },
52 {
53 "Question": "试画出函数 H(Ω) 在 A_c=1，ω_c=0 和 π/2，K=3, 5, 7 时的绝对值图像，函数表达式为：\nH(Ω) = A_c / (2π) * (sin((ω_c - Ω) / 2 * K) / sin((ω_c - Ω) / 2) + sin((ω_c + Ω) / 2 * K) / sin((ω_c + Ω) / 2)) * e^(-jΩ(K - 1) / 2)。",
54 "Response": "首先，定义 Ω 的范围为 0 到 2π，使用 `Ω = 0:0.01:2*pi`。然后，设置参数 A_c = 1，ω_c 取值为 0 和 π/2，K 取值为 3, 5, 7。根据公式计算 H(Ω)，使用 `sin()` 和 `exp()` 函数计算正弦和指数项。最后，使用 `abs()` 函数计算绝对值，并用 `plot()` 绘制不同 K 值对应的绝对值图像。",
55 "Complex_CoT": "1. 定义 Ω 向量，设置为 0 到 2π，步长为 0.01；2. 设置参数 A_c = 1，ω_c 取 0 和 π/2，K 为 3, 5, 7；3. 使用 `sin()` 和 `exp()` 函数计算 H(Ω)；4. 使用 `abs()` 计算 H(Ω) 的绝对值；5. 使用 `plot()` 绘制图像，并使用 `legend()` 区分不同 K 值的图形。"
56 },
57 {
58 "Question": "试画出单位阶跃和冲激函数，使用 MATLAB 命令绘制 u(t-5) 和 δ(t-3) 的波形图，要求在 0<=t<=10 内取步长为 0.001。",
59 "Response": "首先，定义时间向量 t，范围为 0 到 10，步长为 0.001。然后，使用 `heaviside(t-5)` 绘制单位阶跃函数 `u(t-5)`，并使用 `dirac(t-3)` 绘制冲激函数 `δ(t-3)`。对于冲激函数，由于 MATLAB 不会显示无穷大，可以在 `t=3` 处手动标记冲击点，使用 `stem()` 函数来绘制离散信号。",
60 "Complex_CoT": "1. 定义时间向量 t，范围为 0 到 10，步长为 0.001；2. 使用 `heaviside(t-5)` 绘制单位阶跃函数 `u(t-5)`，并使用 `plot()` 绘制；3. 使用 `dirac(t-3)` 绘制冲激函数 `δ(t-3)`，并用 `stem()` 绘制离散点，手动标记 `t=3` 的冲击；4. 使用 `xlabel()`、`ylabel()` 和 `title()` 添加图形标签，使用 `grid on` 打开网格以便观察波形。"
61 },
62 {
63 "Question": "使用 plot 函数作图时，相邻两个点之间直接通过连线连接；而阶梯函数 stairs 在作图时可以使点与点之间形成阶梯形状。为体会这一点，在脚本文件中使用 plot 和 stairs 函数在同一幅图中分别画出从 0 到 2π 之间取值，步长为 0.1π 时的正弦函数图像。",
64 "Response": "首先，定义时间向量 x，从 0 到 2π，步长为 0.1π，即 `x = 0:0.1*pi:2*pi`。然后，使用 `sin(x)` 计算正弦函数值。接着，使用 `plot(x, y)` 绘制平滑曲线，使用 `stairs(x, y)` 绘制阶梯形状的曲线。最后，通过 `xlabel()`、`ylabel()`、`title()` 添加图形标签，使用 `legend()` 区分两种曲线，`grid on` 显示网格。",
65 "Complex_CoT": "1. 定义时间向量 x，设置为 0 到 2π，步长为 0.1π；2. 使用 `y = sin(x)` 计算正弦函数；3. 使用 `plot(x, y)` 绘制平滑曲线，使用 `stairs(x, y)` 绘制阶梯曲线；4. 使用 `xlabel()`、`ylabel()`、`title()` 添加标签，使用 `legend()` 区分两种曲线；5. 使用 `grid on` 打开网格以便更好地观察波形。"
66 },
67 {
68 "Question": "在脚本文件中绘制单位冲激函数和单位阶跃函数，要求在 -5 ≤ t ≤ 10 内取步长为 0.05，t_0 = 1.2，画出该冲激信号和阶跃信号图。",
69 "Response": "首先，定义时间向量 t，范围为 -5 到 10，步长为 0.05。然后，使用 `zeros()` 创建一个全零向量，通过调整 `t_0` 到 `t_0 + Δ` 范围内的值为 `1/Δ`，并用 `stairs()` 绘制冲激函数图。接着，使用 `ones()` 创建一个全为 1 的向量，将 `t_0` 之前的值设为 0，再用 `stairs()` 绘制阶跃函数图。",
70 "Complex_CoT": "1. 定义时间向量 t，设置为 -5 到 10，步长为 0.05；2. 使用 `zeros()` 创建全零向量，按 `t_0` 设置冲激函数值；3. 使用 `stairs()` 绘制冲激函数图；4. 使用 `ones()` 创建全为 1 的向量，修改 `t_0` 之前的值为 0，使用 `stairs()` 绘制阶跃函数图；5. 使用 `xlabel()`、`ylabel()` 和 `title()` 添加图形标签，使用 `grid on` 显示网格。"
71 },
72 {
73 "Question": "在脚本文件中绘制单位协变函数 r(t - t_0)，在 -5 ≤ t ≤ 10 内取步长为 0.05，t_0 = 1.9，使用 `for` 循环和 `if` 语句赋值，并通过 `plot()` 绘制图像。",
74 "Response": "首先，定义时间向量 `t = -5:0.05:10`，然后创建一个与 `t` 长度相同的零向量 `r`。接着，使用 `for` 循环遍历每个 `t`，并在 `if` 语句中检查 `t >= t_0`，若满足条件则设置 `r = t - t_0`，否则保持 `r = 0`。最后，使用 `plot(t, r)` 绘制图形。",
75 "Complex_CoT": "1. 定义时间向量 t，设置为 -5 到 10，步长为 0.05；2. 创建一个全零的向量 `r`，与 `t` 长度相同；3. 使用 `for` 循环遍历 `t`，在 `if` 语句中设置条件 `t >= t_0`，计算协变函数值；4. 使用 `plot()` 绘制图像，并使用 `xlabel()`、`ylabel()` 和 `title()` 添加标签；5. 使用 `grid on` 显示网格，帮助观察波形。"
76 },
77 {
78 "Question": "编写脚本文件和必要的函数以生成信号并画出以下函数的图形：-5 ≤ t ≤ 10，步长 0.05。\ny(t) = 3r(t+3) - 6r(t+1) + 3r(t) - 3u(t-3)",
79 "Response": "首先，定义时间向量 `t = -5:0.05:10`。然后，使用 `zeros()` 函数创建一个全零的向量，计算单位协变函数 `r(t)`，并使用 `heaviside()` 生成单位阶跃函数 `u(t)`。接着，使用公式计算 `y(t) = 3r(t+3) - 6r(t+1) + 3r(t) - 3u(t-3)`，最后通过 `plot(t, y)` 绘制结果图像。",
80 "Complex_CoT": "1. 定义时间向量 `t`，设置为 -5 到 10，步长为 0.05；2. 使用 `zeros()` 创建全零向量，计算 `r(t)` 并用 `heaviside()` 计算 `u(t)`；3. 使用公式 `y(t) = 3r(t+3) - 6r(t+1) + 3r(t) - 3u(t-3)` 计算结果；4. 使用 `plot(t, y)` 绘制图像，并使用 `xlabel()`、`ylabel()` 和 `title()` 添加标签；5. 使用 `grid on` 显示网格，帮助观察波形。"
81 },
82 {
83 "Question": "编写脚本文件生成信号并画出 AM 信号图，给定消息信号 y(t) 和载波信号 x(t) = sin(5πt)，计算 AM 信号 z(t) = (A₀ + y(t)) * sin(5πt)，并验证 AM 信号的包络与消息信号 y(t) 的关系。",
84 "Response": "首先，定义时间向量 t，范围为 -5 ≤ t ≤ 10，步长为 0.05，使用 `t = -5:0.05:10`。然后，使用 `heaviside()` 函数生成阶跃函数，计算消息信号 `y(t)`。接着，使用 `sin(5*pi*t)` 生成载波信号 `x(t)`，并通过 `z(t) = (A₀ + y(t)) * sin(5*pi*t)` 计算 AM 信号。最后，计算包络 `A₀ + y(t)`，并通过 `plot()` 绘制 `y(t)`、`z(t)` 和包络图，使用 `xlabel()`、`ylabel()`、`title()` 添加标签，使用 `grid on` 显示网格。",
85 "Complex_CoT": "1. 定义时间向量 `t`，设置为 -5 到 10，步长为 0.05；2. 使用 `heaviside()` 计算消息信号 `y(t)`；3. 使用 `sin(5*pi*t)` 生成载波信号 `x(t)`；4. 计算 AM 信号 `z(t) = (A₀ + y(t)) * sin(5*pi*t)`；5. 计算包络 `A₀ + y(t)`，使用 `plot()` 绘制图像；6. 使用 `xlabel()`、`ylabel()` 和 `title()` 添加图形标签，使用 `grid on` 显示网格。"
86 },
87 {
88 "Question": "已知信号 f(t) 是幅度为 1，宽度为 4，斜率为 -0.5 的非周期三角信号，试用 MATLAB 命令画出 f(t-2)，f(t-1) 和 f(3t) 的波形图，以观察时移、反折和尺度变换后的信号图像。",
89 "Response": "首先，定义时间向量 `t`，范围为 -5 到 5，步长为 0.1，即 `t = -5:0.1:5`。然后，使用 `tripuls(t, 4)` 生成原始三角波信号。接着，通过 `tripuls(t-2, 4)` 生成 `f(t-2)`，通过 `tripuls(t-1, 4)` 生成 `f(t-1)`，通过 `tripuls(3*t, 4)` 生成 `f(3t)`，分别绘制它们的波形图。",
90 "Complex_CoT": "1. 定义时间向量 `t`，设置为 -5 到 5，步长为 0.1；2. 使用 `tripuls(t, 4)` 生成原始三角波信号；3. 使用 `tripuls(t-2, 4)` 生成 `f(t-2)`，使用 `tripuls(t-1, 4)` 生成 `f(t-1)`，使用 `tripuls(3*t, 4)` 生成 `f(3t)`；4. 使用 `plot()` 绘制波形图，使用 `xlabel()`、`ylabel()` 和 `title()` 添加标签，使用 `grid on` 显示网格。"
91 },
92 {
93 "Question": "已知离散信号 f(n) = n, (0 ≤ n ≤ 6)，试用 MATLAB 命令画出 f(-n)、f(n) 和 f(3n) 的波形图，以观察反折、尺度变换后的序列图像。",
94 "Response": "首先，定义时间向量 `n = 0:6`，表示从 0 到 6 的离散时间序列。然后，使用 `stem(n, f)` 绘制原始信号 `f(n)`。接着，通过 `f(-n)` 来进行信号的反折操作，使用 `stem(-n, f)` 绘制反折后的信号。最后，通过 `f(3*n)` 来进行尺度变换，使用 `stem(3*n, f)` 绘制尺度变换后的信号。",
95 "Complex_CoT": "1. 定义时间向量 `n`，设置为 0 到 6；2. 使用 `stem(n, f)` 绘制原始信号 `f(n)`；3. 使用 `stem(-n, f)` 绘制反折信号 `f(-n)`；4. 使用 `stem(3*n, f)` 绘制尺度变换信号 `f(3n)`；5. 使用 `xlabel()`、`ylabel()`、`title()` 添加图形标签，使用 `grid on` 显示网格。"
96 },
97 {
98 "Question": "用 MATLAB 命令求下列函数关于变量 x 的一阶导数。\n(1) y1 = sin(ax^2)\n(2) y2 = x sin(x log x)",
99 "Response": "首先，使用符号工具箱，定义符号变量 `syms x`。对于第一个函数 `y1 = sin(ax^2)`，使用 `diff(sin(a*x^2), x)` 来计算导数。对于第二个函数 `y2 = x sin(x log x)`，使用链式法则和乘积法则，通过 `diff(x*sin(x*log(x)), x)` 来计算导数。最后，使用 `fplot()` 绘制导数图像，使用 `xlabel()`、`ylabel()` 和 `title()` 为图形添加标签。",
100 "Complex_CoT": "1. 使用 `syms x` 创建符号变量；2. 对于 `y1 = sin(ax^2)`，使用 `diff(sin(a*x^2), x)` 计算导数；3. 对于 `y2 = x*sin(x*log(x))`，使用 `diff(x*sin(x*log(x)), x)` 计算导数；4. 使用 `fplot()` 绘制导数图像，使用 `xlabel()`、`ylabel()` 和 `title()` 添加图形标签；5. 使用 `grid on` 显示网格，帮助观察波形。"
101 },
102 {
103 "Question": "用 MATLAB 命令计算定积分 ∫₀¹ (x·e^x)/(1+x)² dx。",
104 "Response": "首先分析被积函数表达式是否适合使用符号积分。由于函数为 (x·e^x)/(1+x)²，是连续且可导的解析函数，可用 `syms` 定义符号变量，利用 `int()` 函数进行符号积分；若追求数值结果，可定义函数句柄，通过 `integral()` 函数计算定积分。",
105 "Complex_CoT": "1. 判断是否符号积分或数值积分；2. 若使用符号方法，使用 `syms x` 定义变量，`int((x*exp(x))/(1+x)^2, 0, 1)` 求积分；3. 若使用数值方法，使用函数句柄 `f = @(x) (x.*exp(x))./(1+x).^2`，再用 `integral(f, 0, 1)` 求值；4. 可选使用 `fplot()` 观察函数图像；5. 通过 `double()` 获取数值结果，使用 `disp()` 输出结果。"
106 },
107 {
108 "Question": "用 MATLAB 命令计算 y(t) = 3u(t-1) 的微分与积分。",
109 "Response": "首先，使用符号工具箱定义符号变量 t，通过 heaviside(t-1) 构造单位阶跃函数，形成 y(t) = 3*heaviside(t-1)。然后，使用 diff(y, t) 求一阶导数，理论结果为冲激函数 3*δ(t-1)，可通过 dirac 表达。接着使用 int(y, t) 求不定积分，得到一段从 t=1 开始线性增长的函数。可选地，使用 fplot 对函数、导数和积分结果进行可视化分析。",
110 "Complex_CoT": "1. 使用 syms 定义符号变量 t；2. 构造表达式 y(t) = 3*heaviside(t-1)；3. 使用 diff(y, t) 计算导数，理论上得到 3*dirac(t-1)；4. 使用 int(y, t) 计算积分，得到分段线性函数；5. 可使用 fplot 画出原函数、导数和积分图像，辅助理解其数学意义。"
111 },
112 {
113 "Question": "已知 f₁(t) = sin(Ωt), f₂(t) = sin(8Ωt)，试用 MATLAB 命令绘出 f₁(t) + f₂(t) 和 f₁(t)f₂(t) 的波形图，其中 f = Ω / 2π = 1Hz。",
114 "Response": "首先将角频率 Ω 设为 2π，对应频率为 1Hz。定义时间向量 t，覆盖多个周期并具有足够采样精度。构造信号 f₁(t) = sin(2πt)，f₂(t) = sin(16πt)，然后计算 f₁(t) + f₂(t) 与 f₁(t)f₂(t)。使用 plot 函数绘制两组信号波形，通过 subplot 分图显示，配合 xlabel、ylabel 和 title 添加图形标签，辅助观察频率叠加和调制特性。",
115 "Complex_CoT": "1. 设置 Ω = 2π，得到 f₁(t) = sin(2πt)，f₂(t) = sin(16πt)；2. 定义时间向量 t，步长小于 1/最大频率的 10 倍，覆盖至少 1 秒；3. 构造组合信号 f₁+f₂ 和 f₁·f₂；4. 使用 subplot 分别绘制组合信号；5. 添加标签并启用 grid on，分析调制特性及频率关系。"
116 },
117 {
118 "Question": "已知信号 f(t) 是幅度为 1、宽度为 4、斜率为 -0.5 的非周期三角波信号，试用 MATLAB 命令画出 f(t) 奇对称分量和偶对称分量的波形图。",
119 "Response": "首先定义时间向量 t，使其在 -2 到 2 范围内对称取值，并构造满足幅度为 1、宽度为 4、斜率为 -0.5 的非周期三角波信号 f(t)。接着利用 f(-t) = fliplr(f(t)) 构造反转信号，再分别用 (f(t) + f(-t))/2 和 (f(t) - f(-t))/2 计算偶对称分量和奇对称分量。使用 plot 命令绘制三者波形，并加上标题与图例标注，辅助观察。",
120 "Complex_CoT": "1. 构造对称时间向量 t，如 t = -2:0.01:2；2. 用分段线性方式构造非周期三角波 f(t)，满足幅度为 1、斜率为 -0.5、宽度为 4；3. 计算 f(-t)，可通过 fliplr 实现；4. 用公式计算偶对称分量 (f(t)+f(-t))/2 和奇对称分量 (f(t)-f(-t))/2；5. 使用 plot 分别绘制原信号、奇分量、偶分量，并用 xlabel、title、legend 等函数标注图形。"
121 },
122 {
123 "Question": "试用 MATLAB 命令画出下列信号的波形图：\n(1) y₁(t) = 3t² + t + 5\n(2) y₂(t) = t·sin(5πt)·u(t)\n(3) f(t) = e⁻ᵗ·sin(10πt) + e⁻ᵗ⧸²·sin(9πt)",
124 "Response": "针对 y₁(t)，可定义时间向量并直接按表达式计算；对 y₂(t)，应先构造单位阶跃函数 u(t)，将其与 t·sin(5πt) 相乘；对 f(t)，为高频衰减正弦叠加信号，应选择足够高的采样精度和合理时间范围，逐项构建后叠加。全部信号绘图均可使用 plot 函数，并加标题注释。",
125 "Complex_CoT": "1. 使用 linspace 或冒号运算定义时间向量；2. y₁(t) 是多项式，直接向量计算；3. y₂(t) 中的 u(t) 可用 heaviside 函数或逻辑索引生成，构造 t·sin(5πt) 后乘以 u(t)；4. f(t) 为两个带不同指数包络的高频正弦信号之和，分别计算再相加；5. 使用 plot 逐个绘图，配合 xlabel、ylabel、title、legend 与 grid on 完善图像展示。"
126 },
127 {
128 "Question": "已知 f(t) = u(t) - u(t-1)，h(t) = u(t-1) - u(t-3)，试求时域卷积积分 y(t) = f(t) * h(t)。",
129 "Response": "首先根据定义构造两个矩形脉冲信号 f(t) 与 h(t)，分别为在 [0,1) 和 [1,3) 区间为 1 的函数。使用逻辑表达式生成这两个信号的数值形式。卷积可用 MATLAB 中的 conv 函数完成，并乘以时间步长 dt 进行归一化。结果长度为原始信号长度之和减一，需构造对应的新时间向量以匹配卷积输出长度。使用 plot 绘制卷积前后信号，观察系统线性时不变响应。",
130 "Complex_CoT": "1. 定义时间向量 t，如 t = 0:dt:4，步长 dt 选取足够小以近似连续时间；2. 使用 (t >= 0) & (t < 1) 构造 f(t)，(t >= 1) & (t < 3) 构造 h(t)；3. 使用 conv(f, h)*dt 计算 y(t)；4. 构造新时间向量 t_conv 对应卷积结果：t_conv = 0:dt:(length(y)-1)*dt；5. 使用 plot 绘制 f(t)、h(t) 和 y(t)，并使用 xlabel、title、legend 等添加注释。"
131 },
132 {
133 "Question": "用 MATLAB 命令求函数 f(t) = sin(t) 与 g(t) = 0.5*(e^{-t} + e^{-3t}) 的卷积。",
134 "Response": "首先构造时间向量 t，如 0 到 10，步长选为 0.01。分别定义函数 f(t) = sin(t) 和 g(t) = 0.5*(exp(-t) + exp(-3t))。使用 conv 函数计算两者的离散卷积，并乘以时间步长 dt 近似连续卷积。卷积结果的时间长度为原始信号之和减 1，因此需构造新的时间轴用于绘图。使用 plot 函数绘制结果并配合 title、xlabel、legend 辅助分析卷积特性。",
135 "Complex_CoT": "1. 设定 t = 0:dt:T，保证指数信号在有效范围内衰减；2. 构造 f = sin(t)，g = 0.5*(exp(-t) + exp(-3*t))；3. 使用 conv(f, g)*dt 得到卷积结果 y(t)；4. 构造卷积结果时间轴 t_conv = 0:dt:(length(y)-1)*dt；5. 使用 plot 绘制 y(t)，并添加标签和网格以便分析卷积后的信号形状。"
136 },
137 {
138 "Question": "如何利用傅里叶变换实现 f(t) = sin(t) 与 g(t) = 0.5·(e^{-t} + e^{-3t}) 的卷积运算？",
139 "Response": "利用傅里叶变换的卷积定理，可将 f(t) = sin(t) 与 g(t) = 0.5*(e^{-t} + e^{-3t}) 的卷积转换为频域乘法形式。首先定义时间向量并构造 f(t) 与 g(t) 的离散表达式，然后使用 fft() 计算两者的频谱。频谱相乘后，再用 ifft() 执行傅里叶反变换，得到时域卷积结果。注意在频域乘法前应统一信号长度，并在绘图时构建合适的时间轴进行结果对比。",
140 "Complex_CoT": "1. 构造时间向量 t，定义 f(t) = sin(t)，g(t) = 0.5*(exp(-t) + exp(-3*t))；2. 使用 fft() 对两个信号计算频谱，必要时使用 zero-padding 补齐长度避免循环卷积效应；3. 将频域结果逐点相乘，得到 Y(w)；4. 使用 ifft() 得到 y(t)，即 f(t)*g(t) 的卷积结果；5. 使用 plot 绘制傅里叶变换实现的卷积结果并与时域 conv() 结果进行比较验证；6. 使用 abs()、real() 视情况提取有效波形。"
141 },
142 {
143 "Question": "用 MATLAB 命令求：离散信号 x = [3, 11, 7, 0, -1.4, 2, 11, 7, 0, -1.4, 2] 和信号 h = [2, 3, 0, -5, 2, 1] 的卷积和。",
144 "Response": "该问题属于典型的离散时间序列卷积操作，MATLAB 中可直接使用内置函数 conv() 实现，无需手动移位与求和。首先将两个序列以向量形式输入 MATLAB，执行 conv(x, h) 得到输出序列 y。卷积结果长度等于 length(x) + length(h) - 1，可据此构建对应时间索引用于绘图分析。若需可视化离散响应，可使用 stem() 绘制 y(n)，并配合 xlabel、ylabel、title 进行注释。",
145 "Complex_CoT": "1. 将 x 和 h 以一维向量形式输入 MATLAB；2. 使用 conv(x, h) 计算离散卷积；3. 输出序列长度为 length(x) + length(h) - 1，需构造新时间索引 n = 0:(length(y)-1)；4. 使用 stem(n, y) 进行离散绘图；5. 添加 xlabel('n')、ylabel('y(n)')、title('x(n)*h(n)') 标注图像，辅助理解卷积的系统响应含义。"
146 },
147 {
148 "Question": "用 MATLAB 命令求：y(n) = sin(2n) * sin(5n) 的卷积。",
149 "Response": "该问题涉及两个离散正弦信号的线性卷积。首先需要限定离散时间 n 的取值范围，如 n = 0 到 50，分别定义 f1(n) = sin(2n) 与 f2(n) = sin(5n)，注意函数参数为弧度。接着使用 conv(f1, f2) 计算离散卷积，结果长度为 length(f1) + length(f2) - 1，需构造新的时间索引用于绘图。最后使用 stem 或 plot 绘制结果，结合 xlabel、ylabel、title 标注信号特性。",
150 "Complex_CoT": "1. 设定 n = 0:N（如 N=50）构造时间轴；2. 使用 sin(2*n) 和 sin(5*n) 分别生成两个离散正弦信号；3. 使用 conv(f1, f2) 得到卷积结果 y；4. 构造新时间索引 n_y = 0:(length(y)-1)；5. 使用 stem(n_y, y) 绘图，并使用 xlabel、title 等命令添加注释，分析频率调制特征。"
151 },
152 {
153 "Question": "用 MATLAB 命令求：两个等宽矩形序列的卷积和，其中矩形序列 x(n) = u(n) - u(n - 10)。",
154 "Response": "该问题涉及单位矩形序列与自身的卷积，即求其自卷积。由于 x(n) 是从 n = 0 到 n = 9 的单位矩形序列，可以使用 ones(1,10) 快速构造该序列。然后通过 conv() 函数进行离散卷积操作，结果将为等腰三角形，长度为 19。为可视化卷积结果，可构造时间索引 n = 0:18，使用 stem 或 plot 进行绘图，并加上图形标注辅助分析。",
155 "Complex_CoT": "1. 使用 ones(1,10) 生成宽度为 10 的矩形序列；2. 使用 conv(x, x) 实现自卷积，得到长度为 19 的结果；3. 构造对应输出索引 n_y = 0:18；4. 使用 stem(n_y, y) 绘制离散卷积结果；5. 添加 title、xlabel 等标注，观察卷积后三角波的对称性与宽度变化。"
156 },
157 {
158 "Question": "用 MATLAB 命令求：某线性时不变离散系统，设激励 x(n) = (0.5)^n·u(n)，系统单位冲激响应 h(n) = u(n + 5) - u(n - 6)，试确定系统的零状态响应 y(n) = x(n) * h(n)。",
159 "Response": "本题目标是求线性时不变系统在特定激励下的零状态响应。激励 x(n) 是从 n = 0 开始的指数衰减序列，可用 (0.5).^n 构造，取有限长度如 0 到 20 近似表示。h(n) 是定义在 n = -5 到 5 的有限对称矩形序列，通过 u(n+5) - u(n-6) 实现，需构造包含负索引的时间轴并用逻辑索引生成序列。使用 conv(x, h) 实现卷积，输出索引范围为 -5 到 N+5，结果可用 stem 绘图并配合 xlabel、title 注释，辅助观察系统响应行为。",
160 "Complex_CoT": "1. 构造激励序列 x(n) = (0.5).^n，取 n = 0:20；2. 构造 h(n) = u(n+5) - u(n-6)，对应 n = -5:5，用逻辑表达式生成值；3. 使用 conv(x, h) 得到 y(n)；4. 构造输出索引 n_y = -5:(length(x)+5) 匹配卷积长度；5. 使用 stem(n_y, y) 绘图，添加 title、xlabel 等标注观察输出的时域结构。"
161 },
162 {
163 "Question": "已知某 LTI 系统的微分方程为：y''(t) + 2y'(t) + 32y(t) = f'(t) + 16f(t)，试用 MATLAB 命令绘出系统的冲激响应和阶跃响应。",
164 "Response": "本题要求通过 MATLAB 绘制线性时不变系统的冲激响应和阶跃响应。首先需将微分方程转换为传递函数形式 H(s)。对两边拉普拉斯变换（零初始条件）得 H(s) = (s + 16) / (s^2 + 2s + 32)。接着用 MATLAB 的 tf() 函数构造该传递函数模型，再分别调用 impulse() 和 step() 绘制冲激响应和阶跃响应曲线。为控制时间范围与分辨率，可设定仿真时间向量 t。",
165 "Complex_CoT": "1. 分析微分方程并通过拉普拉斯变换求得系统传递函数 H(s) = (s + 16)/(s^2 + 2s + 32)；2. 在 MATLAB 中使用 tf([1 16], [1 2 32]) 构造系统模型；3. 使用 impulse(sys, t) 绘制冲激响应，使用 step(sys, t) 绘制阶跃响应；4. 可设定时间向量 t = linspace(0, Tmax, N) 控制响应时间范围；5. 配合 title、xlabel、ylabel 添加图形注释，使用 grid on 提高可读性。"
166 },
167 {
168 "Question": "当初始状态为 y(0⁻) = 1, y′(0⁻) = 2 时，用 MATLAB 命令确定系统 d²y(t)/dt² + 3dy(t)/dt + 2y(t) = df(t)/dt + 3f(t) 的零输入响应，并画出其波形图；如果输入为 f(t) = e^(-3t)·u(t)，试绘制系统输入及其零状态响应的波形图。",
169 "Response": "该问题涉及连续时间 LTI 系统在初始状态不为零时的零输入响应分析，以及在单位阶跃调制指数信号激励下的零状态响应计算。需先通过拉普拉斯变换将系统微分方程转换为传递函数 H(s)，再分别使用 MATLAB 的 initial() 和 lsim() 函数进行响应仿真。初始状态可用 state-space 或传递函数模型结合设置，输入信号 f(t) = e^(-3t)·u(t) 可用向量加函数表达生成。输出结果用 plot 绘图，并配合图形注释清晰展示响应特性。",
170 "Complex_CoT": "1. 将微分方程转换为传递函数 H(s) = (s + 3)/(s² + 3s + 2)；2. 使用 tf() 构造系统传递函数；3. 利用 initial(sys, x0) 绘制零输入响应，初始条件为 y(0) = 1, y′(0) = 2；4. 构造输入信号 f(t) = exp(-3t)·heaviside(t)，并使用 lsim() 计算零状态响应；5. 设置时间向量 t = linspace(0, T, N) 控制仿真时长；6. 使用 plot(t, y) 和 plot(t, f) 分别绘制响应和输入曲线，添加 xlabel、ylabel、title 和 grid 提高可读性。"
171 },
172 {
173 "Question": "线性时不变离散系统 y(n) - y(n-1) + 0.9y(n-2) = x(n)，试绘出系统冲激响应的波形图。",
174 "Response": "该问题给出了一个二阶常系数线性差分方程形式的 LTI 离散系统，要求在单位冲激信号作用下求其冲激响应。首先将该差分方程转化为系统的传递函数形式 H(z)，再用 MATLAB 中的 filter 或 impz 函数求出单位冲激输入下的系统响应。具体地，传递函数的分子为 [1]，分母为 [1, -1, 0.9]。生成单位冲激输入序列 δ(n) 并通过 filter 函数得到输出。最后使用 stem() 绘图，同时添加 xlabel、ylabel、title 以增强可读性。",
175 "Complex_CoT": "1. 从差分方程中提取系统系数，b = [1]，a = [1, -1, 0.9]；2. 构造输入信号 x 为单位冲激序列，如 x(1)=1，其余为 0；3. 使用 filter(b, a, x) 得到系统输出 y，即冲激响应；4. 构造时间轴 n 与响应 y 对应；5. 使用 stem(n, y) 绘制响应波形，并用 title、xlabel、ylabel 和 grid on 进行标注和美化图像。"
176 },
177 {
178 "Question": "已知 x(n) = cos(0.01π·n²) + sin(0.3π·n)，设系统 y(n) = x(n) + 0.8x(n-1)，试用 MATLAB 命令画出输入 x(n) 和输出 y(n) 的波形图。",
179 "Response": "本题涉及构造一个非线性调频和正弦叠加的离散信号，并求其通过非递归线性系统后的输出。输入信号通过 cos 和 sin 函数构造，时间索引 n 可取 0 到 100。系统结构为 y(n) = x(n) + 0.8x(n-1)，属于有限长单位冲激响应系统，延迟项可通过在 x 前插入 0 或使用 circshift 等方式生成。最终使用 plot 或 stem 绘制输入与输出信号，结合 subplot 实现图形对比。",
180 "Complex_CoT": "1. 构造时间索引 n = 0:100；2. 用 cos(0.01π·n²) + sin(0.3π·n) 计算 x(n)，注意使用 .^ 进行元素平方；3. 生成延迟序列 x(n-1)，如 x_shift = [0, x(1:end-1)]；4. 使用 y = x + 0.8 * x_shift 计算输出；5. 使用 subplot 和 plot/stem 同屏绘制 x(n) 和 y(n)，添加 xlabel、title 和 grid 以增强图形表达效果。"
181 },
182 {
183 "Question": "线性时不变离散系统 y(n) - 0.25y(n-1) + 0.5y(n-2) = x(n) + x(n-1)，试绘出系统冲激响应、阶跃响应的波形图，以及激励 x(n) = (0.5)^n·u(n)，初始条件为零时系统响应 y(n) 的波形图。",
184 "Response": "本题要求分析一个二阶 LTI 离散系统在不同输入下的响应，包括冲激响应、阶跃响应和指数信号激励下的系统输出。系统可通过差分方程系数构造传递函数，其中分子系数为 [1 1]，分母系数为 [1 -0.25 0.5]。冲激响应可用 impz 函数直接获取，阶跃和指数输入则可构造输入向量后使用 filter 函数计算输出。绘图部分建议采用 subplot 将三种响应统一展示，配合 xlabel、title 和 grid 提高可读性。",
185 "Complex_CoT": "1. 根据差分方程写出系统的传递函数 H(z) = (1 + z^(-1)) / (1 - 0.25z^(-1) + 0.5z^(-2))；2. 使用 impz(b, a, N) 得到冲激响应；3. 构造阶跃输入 x1 = ones(1,N)，用 y1 = filter(b, a, x1) 计算阶跃响应；4. 构造指数输入 x2 = (0.5).^n，用 y2 = filter(b, a, x2) 计算响应；5. 用 subplot 和 stem 分别绘图，添加标题和标签便于分析对比。"
186 },
187 {
188 "Question": "观察函数 step() 和 impulse() 的调用格式，假设系统的系统函数为 H(s) = (s² + 3s + 7)/(s⁴ + 4s³ + 6s² + 4s + 1)，可以用几种方法绘出系统的阶跃响应曲线？试分别绘制。",
189 "Response": "该问题要求通过多种建模方式绘制连续时间系统的阶跃响应。已知系统传递函数为高阶有理函数，可使用传递函数模型 tf() 和状态空间模型 ss() 两种方式构建系统模型。调用 step() 函数即可得到系统的阶跃响应，impulse() 可用于冲激响应对比。传递函数方式适合系数已知场景，状态空间方式适合矩阵建模分析。建议分别使用 tf() 与 tf2ss() + ss() 构造系统，利用 step() 绘制响应曲线并对比不同建模方式下的等效性。",
190 "Complex_CoT": "1. 将分子系数写为 b = [1 3 7]，分母系数为 a = [1 4 6 4 1]；2. 使用 tf(b, a) 构建传递函数系统 sys1，调用 step(sys1) 绘制阶跃响应；3. 使用 tf2ss(b, a) 将系统转换为状态空间模型 [A,B,C,D]；4. 使用 ss(A,B,C,D) 构造系统 sys2，调用 step(sys2) 绘制响应；5. 使用 subplot 将两种响应曲线同时显示并添加标题说明差异；6. 可调用 impulse(sys) 分析冲激响应以辅助理解系统特性。"
191 },
192 {
193 "Question": "求频率为 1Hz，占空比为 50% 的周期方波信号的傅里叶级数，用 MATLAB 编程实现其各次谐波的叠加，并观察吉布斯现象。",
194 "Response": "本题要求利用傅里叶级数叠加逼近周期方波信号，并通过 MATLAB 动态展示不同谐波阶数对波形重构精度的影响，从而观察吉布斯现象。方波的傅里叶级数包含奇次正弦谐波分量，幅度按 1/n 衰减。需构造时间向量 t，在 t 上对每一个奇次谐波频率计算正弦波并叠加形成逼近波形。通过增加谐波数观察波形在跃变处的超调现象，即为吉布斯现象的表现。",
195 "Complex_CoT": "1. 设置基础频率 f = 1Hz，构造时间向量 t = linspace(0, 3, 1000)；2. 初始化 y = 0，并设定谐波阶数上限 N（如 N=1, 5, 11, 25 等）；3. 使用 for 循环叠加各奇次谐波分量，表达式为 y = y + (4 / (nπ)) * sin(2πnft)，其中 n 为奇数；4. 对不同的 N 分别绘图展示，用 subplot 并添加 title、legend 区分谐波阶数；5. 观察跃变点附近的超调波动以体现吉布斯现象。"
196 },
197 {
198 "Question": "试用 MATLAB 符号运算求解法求下列信号的傅里叶变换：\\\\ \\ f(t) = e^{-2t}u(t)",
199 "Response": "本题要求对因果指数信号 f(t) = e^{-2t}u(t) 求其傅里叶变换。此类信号在数学上具有解析表达，适合使用 MATLAB 的符号工具箱进行符号变换。首先需通过 syms 定义符号变量 t 和 w，并利用 heaviside(t) 表示单位阶跃函数 u(t)。接着构造符号表达式 f = exp(-2*t)*heaviside(t)，调用 fourier(f, t, w) 实现时域到频域的变换。结果为关于 w 的函数，可使用 simplify 对变换结果化简，也可结合 pretty、disp 输出格式化表达式。进一步可使用 fplot 对实部、虚部或幅频响应进行可视化，辅助理解其频域特性。",
200 "Complex_CoT": "1. 定义符号变量 syms t w 表示时域和频域变量；2. 构造信号表达式 f = exp(-2*t)*heaviside(t)；3. 使用 F = fourier(f, t, w) 得到频域函数；4. 用 simplify(F) 简化表达式，或用 pretty(F) 美观显示；5. 若需图示其频谱特性，可用 fplot(real(F)) 或 abs(F) 可视化频谱幅度；6. 通过观察频域结果分析指数衰减在频域的宽带特性。"
201 },
202 {
203 "Question": "试用 MATLAB 符号运算求解法求下列信号的傅里叶变换：\\\\ \\ f(t) = \\sin(t)",
204 "Response": "该问题要求对非因果连续时间信号 f(t) = sin(t) 求傅里叶变换。该信号在全时域定义，不乘以单位阶跃函数，因此不满足常规 L1 条件，傅里叶变换在经典意义下不一定存在，通常需借助广义函数的定义或在符号计算环境中求解。使用 MATLAB 的 Symbolic Math Toolbox，可以对 sin(t) 进行形式傅里叶变换分析。首先使用 syms 定义 t 与 w 两个符号变量，然后用 fourier(sin(t), t, w) 进行变换，观察其返回形式。该函数返回的结果可能以 δ 函数的形式表达频域特性，需结合数学背景判断其合理性。",
205 "Complex_CoT": "1. 用 syms 定义符号变量 t 和 w；2. 构造函数 f = sin(t)；3. 使用 F = fourier(f, t, w) 调用傅里叶变换；4. 可进一步使用 simplify 化简表达式；5. 若结果为广义函数形式（如 δ(w ± 1)），则说明该信号可视为理想频谱在 ±1Hz 处的集中能量分布；6. 建议对结果进行图形解释，并结合 fplot 绘制频域幅度分布，用于理解理想正弦信号的频谱结构。"
206 },
207 {
208 "Question": "试用 MATLAB 符号运算求解法求下列信号的傅里叶变换：\n3δ(t) - e^{-t}u(t)",
209 "Response": "本题要求使用 MATLAB 符号运算工具，对一个包含冲激函数和单位阶跃函数的复合信号进行傅里叶变换的解析求解。首先，应使用 syms 定义时间变量 t 和频率变量 w，作为傅里叶变换的变量基础。然后，构造原始信号表达式，借助 dirac(t) 表示 δ(t)，heaviside(t) 表示 u(t)，指数函数 e^{-t} 可直接用 exp(-t) 表达，因此原信号可构造为 3*dirac(t) - exp(-t)*heaviside(t)。在表达式构建完成后，使用 fourier 函数对该表达式关于变量 t 进行傅里叶变换。fourier(expr, t, w) 函数将返回频率域变量 w 表示的解析解，若有具体傅里叶参数要求（如使用对称形式），可通过 'FourierParameters' 参数设置。计算结果可通过 simplify 化简，或用 pretty 展示更清晰的表达式。整个过程中应理解每个符号函数在时域与频域的对应关系，为后续学习信号变换奠定基础。",
210 "Complex_CoT": "解题的核心是明确符号运算与信号系统中常见函数的对应关系。δ(t) 是冲激函数，其傅里叶变换为常数项；而 e^{-t}u(t) 是一阶衰减函数，其傅里叶变换是有极点的复频域函数。通过 MATLAB 的符号函数 dirac 和 heaviside，可以准确表达这些理想信号。其次，理解 fourier 函数如何处理这些符号函数至关重要，特别是对单位阶跃和冲激函数的频域变换特性要有扎实掌握。通过对 fourier 的使用掌握，学生将能够处理更复杂的复合信号频域分析，并为系统分析与滤波设计提供理论基础。"
211 },
212 {
213 "Question": "用 MATLAB 符号运算求解法求 F(ω) = 1 / (1 + ω²) 的傅里叶反变换 f(t)。",
214 "Response": "本题目标是通过 MATLAB 的符号运算工具，求解频域函数 F(ω) = 1 / (1 + ω²) 对应的时域函数 f(t)，即傅里叶反变换。首先，应使用 syms 定义符号变量 w 和 t，分别代表频率域变量和时间域变量。然后构造符号表达式 F(w) = 1 / (1 + w^2)。接下来，使用 ifourier(expr, w, t) 函数对该表达式进行傅里叶反变换。若题目或教材对傅里叶变换的参数形式有特殊要求，应使用 'FourierParameters' 参数进行设定，确保结果符号与标准一致。变换后结果可通过 simplify 进行化简，或使用 pretty 展示更清晰的解析表达式。",
215 "Complex_CoT": "本题的关键在于理解频域函数的结构及其对应的时域函数特性。F(ω) = 1 / (1 + ω²) 是典型的低通滤波器频率响应形式，其傅里叶反变换结果应与双边衰减信号相关。通过 ifourier 函数可实现频到时的解析转换，过程中需要确保变量和函数表达式都为符号形式。理解傅里叶反变换操作对掌握频域分析在系统建模与信号还原中的应用具有重要意义，进一步可推广至更复杂频域函数的处理。"
216 },
217 {
218 "Question": "用 MATLAB 符号运算求 f(t) = e^{-t^2/2}·u(t) 的傅里叶变换，并用傅里叶反变换算法验证。",
219 "Response": "本题要求用 MATLAB 符号运算工具求解一个单边高斯信号的傅里叶变换，并通过反变换验证其正确性。首先需使用 syms 定义符号变量 t 和 w，构造信号表达式 exp(-t^2/2)*heaviside(t)。该信号非偶对称，傅里叶变换后结果通常较复杂，可能包含误差函数等特殊符号。使用 fourier(expr, t, w) 对其求频域表达式 F(w)。随后，使用 ifourier(F, w, t) 对频域结果进行反变换，返回时域表达式。为验证反变换是否与原始函数一致，应使用 simplify 对反变换结果和原函数差值进行化简，判断是否等价。",
220 "Complex_CoT": "本题中的信号是截断的高斯型函数，即标准高斯函数与单位阶跃函数的乘积，因其不具偶对称性，傅里叶变换后不再是高斯函数，而是复合结构，结果可能含有特殊函数如 erf。通过符号傅里叶变换可深入理解信号定义域对频域结果的影响，特别是单边与双边信号在频域表示上的差异。反变换验证环节要求学生掌握从频域重建原始时域函数的技巧，进一步加深对时频双向变换关系的理解，是符号运算在信号分析中的重要应用。"
221 },
222 {
223 "Question": "使用 MATLAB 命令绘制单边指数信号 f(t) = e^{-a t} u(t) 的幅度谱和相位谱，并通过改变尺度参数 a 的数值，观察尺度变换对频谱图的影响。",
224 "Response": "本题要求通过 MATLAB 数值仿真分析单边指数信号在不同时间尺度参数 a 下的频域特性。首先需构造合适的时间向量 t，并定义信号表达式为 exp(-a*t).*heaviside(t)，其中 a 为可调尺度因子。为获得频谱，使用 fft 对信号进行傅里叶变换，并用 fftshift 将频谱居中，再通过 abs 提取幅度谱，通过 angle 提取相位谱。绘图部分可用 subplot 将多个频谱图展示在同一窗口下，便于观察不同 a 值下的频域变化趋势。通过实验对比不同 a 值时频谱的宽度和形态差异，可清晰观察时间尺度对频率扩展的影响。",
225 "Complex_CoT": "单边指数信号的时间尺度变化直接影响其频域特性。a 越大，信号衰减越快，频谱越宽；a 越小，信号在时域上更“缓慢”，频谱越集中。这是傅里叶变换中时间与频率互为倒数关系的直接体现。通过实际绘图分析，不仅能验证该理论，还能帮助学生建立时频关系的直观理解。在实现过程中，还应注意频率轴的构建与单位一致性，以保证频谱图具有物理意义。"
226 },
227 {
228 "Question": "某系统的频响函数 H(ω) = 1 / (1 + i·10·ω)，试画出其对数幅频特性与相频特性。",
229 "Response": "本题要求利用 MATLAB 绘制系统频响函数 H(ω) = 1 / (1 + i·10·ω) 的对数幅频特性和相频特性。首先应构造频率向量 ω，可使用 linspace 或 logspace 在适当频率范围内（如 0.1 到 100 rad/s）生成频率样本点。然后使用复数表达式计算 H(ω)，即 H = 1 ./ (1 + 1i*10*omega)。对数幅频特性通过 abs 取模后转为 dB 形式，即 20*log10(abs(H))；相频特性通过 angle 获取相位角（弧度），可用 rad2deg 转为角度。最终使用 semilogx 绘制幅频和相频图，必要时用 subplot 分别展示两图。",
230 "Complex_CoT": "该频响函数表示一阶低通系统，其幅度在低频趋于 0 dB，高频逐渐衰减约 -20 dB/dec；相位由 0 向 -90°过渡。通过对数坐标绘图可清晰观察这一变化趋势。频域分析是理解系统动态特性的关键手段，MATLAB 中通过复数表达、幅度计算和相位提取，可以高效完成频率响应分析。掌握该方法对后续滤波器设计、系统建模与控制分析具有重要意义。"
231 },
232 {
233 "Question": "某系统的频响函数 H(ω) = 1 / (1 + i·10·ω)，试画出其对数幅频特性与相频特性。",
234 "Response": "本题要求利用 MATLAB 绘制系统频响函数 H(ω) = 1 / (1 + i·10·ω) 的对数幅频图和相频图。需先构造频率变量向量 ω，常用 logspace 创建对数分布的频率点，使频谱覆盖宽频范围。频响函数 H(ω) 可用复数表达式实现，注意使用逐元素运算。接着计算 H(ω) 的模值并转换为分贝（dB）用于对数幅频图，计算相位角用于相频图。最后使用 semilogx 绘图，并通过 subplot 同时展示幅度与相位特性。",
235 "Complex_CoT": "思考过程如下：首先要明确频响函数 H(ω) 的数学结构是复函数，含有虚数单位 i，因此在 MATLAB 中应使用 1i 表示虚数。由于幅频和相频特性是函数在频率维度上的变化趋势，应在较宽的频率范围上进行采样，logspace 是常用生成对数尺度频率向量的函数，可保证频率覆盖从低到高的多个数量级。接下来计算 H(ω) 时要使用逐元素除法 ./，确保向量运算正确执行。幅度通过 abs(H) 得到，再经 20*log10 转换为 dB，表示信号通过系统后的增益变化。相位通过 angle(H) 得到，为弧度制，可用 rad2deg 转换为角度，便于直观理解系统相位滞后情况。最后绘图阶段，为了清晰对比幅频与相频特性，应使用 subplot 将两图分别显示，横轴使用 semilogx 便于呈现对数尺度频率响应。通过上述过程，能够全面分析该一阶系统的频域行为，包括低通特性、截止频率附近的相位变化等，是理解系统频率响应的重要步骤。"
236 },
237 {
238 "Question": "试画出频响函数 H(ω) = i·ω / (1 + i·ω / ω_c) 的对数幅频特性。",
239 "Response": "本题要求在 MATLAB 中绘制频率响应函数 H(ω) = i·ω / (1 + i·ω / ω_c) 的对数幅频响应。首先需定义频率变量 ω，推荐使用 logspace 生成对数尺度频率向量以便观察宽频特性。设定常数 ω_c，例如 ω_c = 10，随后构造 H(ω) 的复数表达式。使用 abs 函数求复数模，再通过 20*log10(...) 转换为分贝单位，得到对数幅频数据。使用 semilogx 绘制频率（横轴为对数坐标）与幅度（单位为 dB）之间的关系曲线，并添加合适的图例和标签。",
240 "Complex_CoT": "思考过程如下：该频响函数本质为一阶高通滤波器，其数学结构包含虚数单位和频率变量的比例项。首先分析函数形式，分子为 i·ω 表示随频率线性增长，分母为 1 + i·ω / ω_c，表征高通特性。ω 很小时分母趋近于 1，幅度随 ω 增大而线性增长；ω 很大时分母近似为 i·ω / ω_c，整体幅度趋于常数。为了在 MATLAB 中实现，需要构造足够密集的频率点来观察低频与高频区域的变化，logspace 可避免线性分布在低频段过于稀疏。表达式中含有复数除法，必须采用逐元素运算方式构造 H(ω)。完成复频函数定义后，用 abs 提取幅值并转为 dB，使用 semilogx 绘图，这样可以清晰展示滤波器的幅频响应曲线及其频率选择性特征。整个过程中应注重表达式的数值稳定性及图像的可读性，如添加网格线、合适的频率刻度和图例说明。"
241 },
242 {
243 "Question": "已知升余弦脉冲信号为 f(t) = (E/2)·[1 + cos(πt/τ)], |t| ≤ τ。用 MATLAB 编程实现该信号经冲激脉冲抽样后得到的抽样信号 f_s(t) 及其频谱。",
244 "Response": "本题要求利用 MATLAB 实现升余弦脉冲信号 f(t) 的数值构造，并通过冲激脉冲抽样方法得到其抽样信号 f_s(t)，最后计算抽样信号的频谱。首先应生成对称时间向量 t（如 -2τ 到 2τ），构造原始信号 f(t)，在 |t| > τ 区间内赋值为 0，在 |t| ≤ τ 区间内按升余弦表达式定义。随后设定抽样间隔 Ts，通过周期性选取 t 向量中的部分点实现冲激抽样，即构造稀疏序列 f_s(t)。抽样后使用 fft 求其频谱，配合 fftshift 实现频谱居中，abs 得到幅度谱。最后通过 plot 或 stem 绘制原始信号、抽样信号及频谱图。",
245 "Complex_CoT": "思考过程如下：升余弦脉冲在时间上有限支持，必须精确模拟 |t| ≤ τ 区间，因此构造时间向量时需涵盖完整信号形状。信号表达式中包含 cos(πt/τ)，应确保 t/τ 是向量化运算，并通过逻辑索引实现区间限制，即仅在满足条件的 t 上赋值。实现冲激抽样关键是模拟理想冲激序列的作用，可通过在固定采样点保留信号值、其余点置零的方法实现，常用方法包括构造采样掩码或直接赋值。频谱分析需使用 fft 得到复频谱，再通过 fftshift 中心化频率，并以 abs 计算幅度谱。频率轴需依据采样率 Fs 构造，通常为 linspace(-Fs/2, Fs/2, N)。该题不仅涉及时间域信号构造与抽样，还涵盖了频谱计算与离散频率轴构建，整体流程涉及信号处理关键操作，掌握后有助于理解抽样定理及其对频谱的影响。"
246 },
247 {
248 "Question": "已知升余弦脉冲信号为 f(t) = (E/2)[1 + cos(πt/τ)], |t| ≤ τ，设该信号的截止频率为 ωₘ = 2，抽样间隔 Tₛ = 1，采用截止频率为 ω_c = 1.2 × ωₘ 的低通滤波器对抽样信号进行滤波重建，得到重建信号 f̂(t)，并计算其与原始信号 f(t) 之间的绝对误差。",
249 "Response": "本题分为三个步骤：首先构造升余弦脉冲信号 f(t)，在区间 |t| ≤ τ 内用升余弦公式定义，在其他时间置零；其次以抽样间隔 Tₛ = 1 对该信号进行冲激抽样，生成抽样信号 fₛ(t)；最后设计一个截止频率为 ω_c = 2.4 的低通滤波器，对 fₛ(t) 进行卷积滤波，实现重建，得到 f̂(t)。最后将 f̂(t) 与原始信号 f(t) 在对应时间点上作差，取绝对值并统计其误差。",
250 "Complex_CoT": "第一步需构造时间轴 t，覆盖至少 [-τ, τ] 范围，保证完整描述升余弦信号。定义信号 f(t) 时，需在 |t| ≤ τ 条件下应用表达式 (E/2)[1 + cos(πt/τ)]，可用逻辑索引实现。在抽样阶段，设抽样间隔 Tₛ = 1，根据 t 的索引构造抽样掩码，每隔 Tₛ 时间点保留信号值，其余点设为零，得到抽样信号 fₛ(t)。重建时使用低通滤波器滤除高频混叠成分，滤波器截止频率设为 ω_c = 2.4 rad/s。设计滤波器时需注意频率归一化，常用 fir1(N, wc/π) 创建有限长 FIR 滤波器，或直接使用 lowpass 函数。滤波通过 conv 实现抽样信号与滤波器冲激响应的卷积，输出即为重建信号 f̂(t)。最后对 f̂(t) 与 f(t) 逐点比较，计算 abs(f̂(t) - f(t)) 得到误差向量，结合 mean 或 max 得出误差指标。该过程涵盖了时域信号建模、离散抽样、滤波器设计与重建误差分析，是数字信号处理系统仿真的完整实现路径。"
251 },
252 {
253 "Question": "已知调制信号 f(t) = cos(10πt)，载波信号 f_c(t) = cos(80πt)，试编程画出调制与解调过程中的波形图与频谱图。",
254 "Response": "本题目标是利用 MATLAB 实现调制与解调过程的完整仿真，包括时域波形和频域频谱的可视化。首先构造调制信号和载波信号，将二者相乘得到调制信号。然后绘制调制前后的波形图和频谱图。接着在解调环节中将调制信号再次与载波相乘，并设计低通滤波器对结果进行滤波，最后对比解调后的信号与原始调制信号。",
255 "Complex_CoT": "思考过程包括以下几个步骤：首先应设置合理的时间向量 t，考虑调制与载波频率分别为 5Hz 和 40Hz（因为 cos(10πt) 对应频率为 5Hz），采样频率 fs 应至少为 10 倍载波频率以避免混叠，建议 fs ≥ 400Hz。调制信号可直接用 cos(10π*t) 生成，载波信号用 cos(80π*t)。调制过程通过 s(t) = f(t)·f_c(t) 得到，频域上将原始信号频谱搬移到载波频率附近。使用 fft 结合 fftshift 获取频谱，abs 提取幅度信息，并构造频率轴用于绘图。解调过程再次与载波相乘，即 s(t)·f_c(t)，此操作在频域上产生一个低频分量和一个高频分量。为保留低频分量（原信号），必须使用低通滤波器滤除高频部分。滤波器可用 lowpass 函数设计，设定合适截止频率，如 1.5×原信号频率。滤波后即得解调信号，再次使用 fft 和时域图展示解调结果。整个过程中应特别注意时间与频率轴的构造、滤波器设计参数、fft 结果归一化等细节，确保结果准确可视。"
256 },
257 {
258 "Question": "运用 MATLAB 求单位斜变函数和 f(t) = t·e^(–2t)·cos(3t)·u(t) 的拉普拉斯变换。",
259 "Response": "本题需要使用 MATLAB 的符号计算功能分别求解单位斜变函数和复合因果信号 f(t) = t·e^(–2t)·cos(3t)·u(t) 的拉普拉斯变换。需借助 laplace 函数进行解析变换，heaviside 表示单位阶跃函数以保证因果性，最终得到 F(s) 表达式。",
260 "Complex_CoT": "首先，明确单位斜变函数可表示为 r(t) = t·u(t)，它是时间变量 t 与单位阶跃函数 u(t) 的乘积。在 MATLAB 中应使用符号变量 t 表示时间，用 heaviside(t) 表示单位阶跃函数，从而完整定义因果性。定义完 r(t) 后使用 laplace(r(t), t, s) 即可求其拉普拉斯变换。其次，对 f(t) = t·e^(–2t)·cos(3t)·u(t) 进行拉普拉斯变换，其形式是多个函数的乘积，必须使用符号乘法定义清晰的函数表达式。符号表达为：t * exp(-2*t) * cos(3*t) * heaviside(t)。同样调用 laplace 函数进行变换。变换完成后，为使结果更紧凑，可使用 simplify 对表达式进行化简。此类问题考查学生对因果信号、指数衰减、调制余弦等基本函数的频域映射关系理解，掌握 laplace 与 heaviside 的组合用法，是建立微分系统建模与频域分析能力的基础。"
261 },
262 {
263 "Question": "设一像函数 F(s) = (5s - 1)/(s³ - 3s - 2)，试确定其部分分式展开式中的待定系数与极点。",
264 "Response": "本题需使用 MATLAB 对有理函数 F(s) = (5s - 1)/(s³ - 3s - 2) 进行部分分式展开，以求出分式中每一项的待定系数，并确定分母多项式的极点。核心是借助 partfrac 函数进行分解，并配合 factor 或 roots 函数求取极点位置。",
265 "Complex_CoT": "首先将 F(s) 表达为符号形式，需使用 syms 定义变量 s，并构造 F = (5*s - 1)/(s^3 - 3*s - 2)。分母是三次多项式，先使用 factor(F) 或 roots(denominator(F)) 进行因式分解或求根，得到极点位置。然后使用 partfrac(F) 对其进行部分分式展开。该函数会自动匹配出每个极点对应的分式结构，并解出每项前的待定系数（即残差），例如形如 A/(s - p1) + B/(s - p2) + C/(s - p3)。通过查看展开结果可直接获取系数 A, B, C 的符号表达。若要求更清晰展示各项，可使用 simplify 或 pretty 美化表达式。此题过程体现了符号工具在系统函数极点分析与形式分解方面的能力，为后续系统响应计算及拉普拉斯反变换打下基础。"
266 },
267 {
268 "Question": "求下列像函数的拉普拉斯反变换：F₁(s) = (3s + 2)/(s² + 3s + 2)，F₂(s) = (3s² + 2s + 5)/(s³ + 12s² + 44s + 48)。",
269 "Response": "本题需使用 MATLAB 的符号运算功能对两个有理分式函数进行拉普拉斯反变换，得到其对应的时域表达式。使用 ilaplace 函数可以快速获得解析结果，必要时结合 partfrac 展开帮助理解组成部分。适当使用 simplify 可使结果更加清晰。",
270 "Complex_CoT": "首先使用 syms 定义符号变量 s 和 t，然后构造两个待处理的符号表达式 F₁(s) 和 F₂(s)。这两个函数都是有理形式，适用于直接使用 ilaplace 进行反变换。F₁(s) 的分母为二次多项式，F₂(s) 的分母为三次多项式，可能具有实数或复数极点，决定了时间响应中是否包含衰减振荡项。ilaplace(F₁, s, t) 和 ilaplace(F₂, s, t) 可分别得到其在 t 域的解析解，若需要进一步观察解构形式，可先使用 partfrac 将表达式分解成多个简单分式，再逐项进行反变换。此方法有助于理解每个极点对应的时间特性，如指数衰减项、正弦余弦项等。最终可使用 simplify 化简结果，并借助 pretty 以更加直观的格式显示出来。通过这一过程，不仅能够获得反变换结果，还能深入理解频域结构如何映射到时域响应的具体形式。"
271 },
272 {
273 "Question": "求下列像函数的拉普拉斯反变换：F(s) = (s + 3)/(s³ + 5s² + 12s + 8)。",
274 "Response": "本题要求对 F(s) = (s + 3)/(s³ + 5s² + 12s + 8) 进行拉普拉斯反变换，可借助 MATLAB 的符号计算工具，使用 ilaplace 函数直接求解时间域表达式，并结合部分分式展开进一步理解各分量的组成结构。",
275 "Complex_CoT": "首先使用 syms 定义符号变量 s 和 t，并构造出表达式 F(s) = (s + 3)/(s³ + 5s² + 12s + 8)。由于这是一个有理分式，适合使用 ilaplace 函数进行解析求解。该函数会将频域表达式转换为时间域函数，自动识别分母的极点并输出对应的指数或振荡项。为帮助理解结构，也可以先使用 partfrac 对 F(s) 进行部分分式展开，将其分解为多个低阶项如 A/(s + a)、B/(s + b) 等，这些项的反变换在拉普拉斯对照表中都有标准形式。随后对每一分项分别使用 ilaplace 进行求解，再相加得到完整的时间函数。最后建议使用 simplify 化简结果，使表达式更为简洁；也可以使用 pretty 美观显示最终形式。通过这一过程，不仅能获得精确结果，还能深入理解极点与系统时间响应之间的对应关系，有助于掌握从频域模型构造时间响应的基本思维方式。"
276 },
277 {
278 "Question": "为了讨论系统函数 H(s) 的零极点对冲激响应 h(t) 的影响，试用 MATLAB 编程绘制以下两个系统的冲激响应与阶跃响应，并分析差异：H₁(s) = 5/(5s² + s + 5)，H₂(s) = (4s + 5)/(5s² + s + 5)。",
279 "Response": "本题旨在通过对比两个传递函数具有相同极点但不同零点的系统，来分析零点对系统时域响应的影响。可通过 MATLAB 定义两个系统模型 H₁ 和 H₂，分别使用 impulse 和 step 函数计算冲激响应和阶跃响应，并绘图进行对比观察。从响应图中可以看出，虽然两个系统的稳定性一致（因极点相同），但其响应的速率、波形形态和过渡特性会因零点的不同而产生显著差异。",
280 "Complex_CoT": "首先使用 tf 函数在 MATLAB 中定义两个系统的传递函数。H₁ 的分子为常数 5，而 H₂ 的分子为 4s + 5，分母相同，因此两个系统拥有相同的极点，即具有相同的本质动态特性。然后设置统一的仿真时间向量，如 t = 0:0.01:10，确保响应过程充分展开。利用 impulse(sys, t) 分别求出两个系统的冲激响应，再利用 step(sys, t) 求出其阶跃响应。绘图时使用 subplot 将冲激响应和阶跃响应分开展示，使用 legend 和 title 添加图例与标题，便于直观对比。从分析角度看，极点决定系统稳定性和响应形式，而零点会影响系统的响应速度、初始斜率、峰值位置等动态特性。通过两个系统在相同极点下的不同响应差异，可以清晰体现零点对系统时域行为的调制作用，进一步强化对系统结构与动态性能之间关系的理解。"
281 },
282 {
283 "Question": "已知 H(s) = (s² - 1)/(s³ + 2s² + 3s + 2)，画出 H(s) 的零极点图，将其展开为部分分式，并求出 h(t)。",
284 "Response": "本题要求分析系统函数 H(s) 的零极点结构，进行部分分式展开，并求解对应的单位冲激响应 h(t)。在 MATLAB 中，可使用 tf 定义 H(s)，利用 pzmap 绘制极点零点图，通过 partfrac 对其进行部分分式展开，再使用 ilaplace 求出时域响应表达式。最终可用 fplot 将 h(t) 可视化。",
285 "Complex_CoT": "首先将 H(s) 表达为分子 num = [1 0 -1] 和分母 den = [1 2 3 2]，在 MATLAB 中使用 tf(num, den) 构建传递函数模型。调用 pzmap(sys) 可在复平面绘制极点与零点，判断系统稳定性并观察其模态分布。接着，为便于求解拉普拉斯反变换，使用 syms s t 定义符号变量，构造符号表达式 H(s)，调用 partfrac(H) 展开为多个简单分式项。由于每一项都对应于标准的拉普拉斯变换形式，可以分别使用 ilaplace 进行反变换，再将结果求和得到 h(t)。为便于表达和图形展示，应用 simplify 化简结果，使用 fplot(expr, [0, Tmax]) 生成 h(t) 的波形图。通过这个完整流程，可实现频域到时域的映射，深入理解系统结构对响应形式的决定性影响。"
286 },
287 {
288 "Question": "如何用 MATLAB 分析一个已知传递函数 H(s) 的 LTI 系统的冲激响应、阶跃响应以及对特定输入（如 e^{-t}u(t) 和 6cos(t)u(t)）的响应，并绘制波形图？",
289 "Response": "可通过 MATLAB 中的控制系统工具箱函数，对传递函数建模后分别调用 impulse、step 和 lsim 函数计算不同输入下的系统响应，再使用 plot 命令绘图。关键步骤包括建立系统模型、构造时间向量、设置信号输入及调用相应函数得到响应波形。",
290 "Complex_CoT": "首先，需要理解题目所给 H(s) 是一个线性时不变系统的传递函数，形式为有理分式，需提取分子与分母多项式系数，使用 tf(num, den) 函数建模。然后，冲激响应可直接使用 impulse(sys, t) 获得，阶跃响应使用 step(sys, t)。这两个函数均基于系统的零初始状态，自动使用单位输入信号，返回系统的时间响应。接着，对任意输入信号如 e^{-t}·u(t) 和 6·cos(t)·u(t)，需先构造时间向量 t，如 0:0.01:10，然后用元素方式定义输入信号向量 u。调用 lsim(sys, u, t) 计算该信号下系统的零状态响应。绘图时，使用 plot(t, response)，并可配合 title、xlabel、ylabel、legend 做图形标注。为避免误差，所有信号定义需与时间向量一一对应。该流程体现了从系统建模、输入信号设计到响应分析和可视化的完整分析链路。"
291 },
292 {
293 "Question": "如何用 MATLAB 分析系统函数 H(s)=1000(s+2)/[(s+10)(s+50)] 的幅频特性与相频特性，并说明其滤波功能？",
294 "Response": "需建立系统传递函数模型，采用 MATLAB 的 bode 函数绘制幅频与相频特性曲线，通过观察系统在不同频率下的增益与相位变化判断其滤波性质。该系统属于带有一个零点和两个极点的低通滤波器结构，具有抑制高频、保留低频的能力。",
295 "Complex_CoT": "首先提取系统传递函数 H(s) 的分子与分母多项式系数，使用 tf(num, den) 建立连续时间系统模型。然后利用 bode(sys) 函数绘制系统的对数幅频图与相频图，横轴为频率（对数坐标），纵轴分别为增益（dB）和相位（度）。通过图形分析，若系统在低频段幅值较大、高频段幅值逐渐衰减，说明其具有低通滤波特性。为了更清晰理解频域响应，也可使用 bodemag(sys) 单独绘制幅频特性。若需查看特定频点响应，可配合 freqresp 或 bodeplot 函数自定义频率范围。判断系统类型时应结合极零结构分析：本题中，零点在 s=-2，两个极点在 s=-10 和 s=-50，说明系统增益随频率升高而减小，具有明显的低通特征。因此，系统可用于滤除高频干扰，保留低频信号。该分析过程完整体现了频域建模、特性提取与滤波性质判断的综合应用思路。"
296 },
297 {
298 "Question": "已知 X(z) = z / (z - 1)^2, |z| > 1，试求原序列 x(n)，如何使用 MATLAB 得到结果？",
299 "Response": "使用 syms 定义符号变量 z 和 n，构造表达式 X = z / (z - 1)^2；使用 iztrans(X, z, n) 求出时域序列 x(n)；用 ztrans(x(n), n, z) 做正变换，并用 simplify 化简，与原始 X(z) 比较以验证结果。",
300 "Complex_CoT": "先识别 X(z) = z / (z - 1)^2 属于已知标准形式，对应 x(n) = (n + 1)·u(n)。可借助 MATLAB 的符号计算功能求解。首先用 syms 定义变量 z 和 n，便于后续表达式操作；再构造 X(z) 的表达式 z / (z - 1)^2；调用 iztrans 进行反 Z 变换，返回解析表达式 x(n)；为验证结果正确性，可再对得到的 x(n) 使用 ztrans 求其 Z 变换，并通过 simplify 化简与题设的 X(z) 比较，判断是否一致。"
301 },
302 {
303 "Question": "设 X(z) = (z^3 + 1)/(z^3 - z^2 - z - 2)，试确定其部分分式展开式中的待定系数与极点，如何使用 MATLAB 得到结果？",
304 "Response": "将分子和分母分别写成系数向量形式：分子为 [1 0 0 1]，分母为 [1 -1 -1 -2]；使用 residue 函数：[r, p, k] = residue(b, a)，其中 b 和 a 分别为分子和分母系数向量；r 为待定系数，p 为极点，k 为首部多项式（若有）；读取 r 和 p 即可得出所需结果。",
305 "Complex_CoT": "题目要求将 X(z) 展开为部分分式形式，关键是分解为若干一阶因子的线性组合。首先观察 X(z) 为有理函数，分子分母都是三阶多项式，阶数相等。考虑使用 MATLAB 中 residue 函数求部分分式展开，需先将分子 z^3 + 1 和分母 z^3 - z^2 - z - 2 分别转换为多项式系数向量，即 [1 0 0 1] 和 [1 -1 -1 -2]。将其代入 residue 函数，可得到分式展开中的系数向量 r 和对应极点向量 p，若有多项式部分则由 k 给出。r 和 p 对应于展开项 r(i)/(z - p(i))，每个 p(i) 是一个极点。利用这些信息可重建部分分式展开式并进一步分析系统性质。"
306 },
307 {
308 "Question": "设 x₁(n) = [8 + 2(0.5)^n - 9(0.75)^n]·u(n), x₂(n) = n²·u(n)，求其 Z 变换；设 X₁(z) = z / (z - 2), X₂(z) = (z³ + 2z² + 1) / [z(z - 1)(z - 0.5)]，求其反变换，如何使用 MATLAB 实现？",
309 "Response": "使用 syms 定义符号变量 n 和 z；构造 x₁(n) 与 x₂(n) 的表达式后，使用 ztrans 计算其 Z 变换；对于 X₁(z) 与 X₂(z)，构造符号表达式后使用 iztrans 计算其反变换；可用 simplify 化简表达式验证结果。",
310 "Complex_CoT": "x₁(n) 为多个幂函数的线性组合，Z 变换具有线性性，(a^n)·u(n) 的 Z 变换为 z / (z - a)，因此对 x₁(n) 分项后分别使用 ztrans 求变换再相加；x₂(n) 为 n²·u(n)，可直接使用 ztrans 求其变换。对第(2)问，X₁(z) 是基本型，可直接使用 iztrans 求反变换；X₂(z) 为真分式，分母含三个线性因子，iztrans 可自动进行部分分式展开并返回对应时域表达式。使用 syms 定义符号变量 n 和 z 后构造函数表达式，调用 ztrans 或 iztrans 获得目标变换结果；必要时使用 simplify 化简以确认表达式形式与理论一致。"
311 },
312 {
313 "Question": "已知 H(z) = 4z^{-1} / (4 - 9z^{-1} + 2z^{-2})，画出 H(z) 的零极点图，将其展开为部分分式，并求出 h(n)，如何使用 MATLAB 实现？",
314 "Response": "先将 H(z) 中的负幂形式转换为正幂形式，对分子分母分别乘 z² 得到 H(z) = 4z / (2z² - 9z + 4)，从中提取分子和分母的系数向量：b = [0 4 0]，a = [2 -9 4]；使用 zplane(b, a) 绘制零极点图；使用 residue(b, a) 得到部分分式展开的系数和极点；构造符号表达式后使用 iztrans 计算 h(n)。",
315 "Complex_CoT": "原式 H(z) 是负幂形式，使用 MATLAB 前需将其转为常规正幂分式，便于提取系数向量。将 H(z) 的分子分母同乘 z² 得到 H(z) = 4z / (2z² - 9z + 4)，此时可直接读取分子系数 b = [0 4 0]，分母系数 a = [2 -9 4]。为绘制零极点图，调用 zplane(b, a) 观察系统的极点与零点分布。接着使用 residue(b, a) 进行部分分式展开，返回的 r 和 p 分别为每项的系数和极点。最后，若要求 h(n)，可使用 syms 定义符号变量后构造 H(z) 的表达式，用 iztrans 求其反 Z 变换，得到单位脉冲响应 h(n)。此过程可从数值分析与符号运算两个角度完成。"
316 },
317 {
318 "Question": "已知 H(z) = (3z² - 2z) / (z² - 0.9z + 0.25)，试绘出其单位冲激响应和单位阶跃响应的波形图，如何使用 MATLAB 实现？",
319 "Response": "将分子和分母分别写成系数向量 b = [3 -2 0]，a = [1 -0.9 0.25]；使用 impz(b, a, N) 得到前 N 项单位冲激响应；构造长度为 N 的单位阶跃序列 u = ones(1, N)，用 filter(b, a, u) 得到单位阶跃响应；分别使用 stem 绘制两个响应的波形图。",
320 "Complex_CoT": "根据系统函数 H(z)，先将其转化为 MATLAB 可用的系数向量形式，其中分子为 3z² - 2z，可写为 [3 -2 0]，分母为 z² - 0.9z + 0.25，即 [1 -0.9 0.25]。为了求单位冲激响应，可以使用 impz 函数输入 b、a 及所需的长度 N 得到响应序列，也可以构造单位冲激输入 delta = [1 0 0 ...]，再用 filter 函数求系统输出作为 h(n)。单位阶跃响应是对全 1 序列进行滤波后的系统输出，即构造 u = ones(1, N)，用 filter(b, a, u) 得到响应 s(n)。最后使用 stem 函数分别绘制 h(n) 和 s(n)，并加上坐标轴和标题以便分析。"
321 },
322 {
323 "Question": "设 H(z) = 1 / (z - 0.5)，试绘制系统的频率响应函数，如何使用 MATLAB 实现？",
324 "Response": "将 H(z) 表示为系数向量形式，分子 b = [1]，分母 a = [1 -0.5]；使用 freqz(b, a, N) 计算频率响应，其中 N 是频率点数；freqz 会输出复数频率响应，可直接绘制幅度和相位响应图；也可使用 abs 和 angle 提取频率响应的幅度与相位，再用 plot 绘制自定义图像。",
325 "Complex_CoT": "H(z) 是一阶系统函数，频率响应定义为 H(e^{jω})，即将 z 替换为单位圆上的复数 e^{jω}，通常 ω ∈ [0, π]。在 MATLAB 中，需将系统函数表示为系数向量形式，分子为 [1]，分母为 [1 -0.5]。使用 freqz 函数可以在单位圆上对该系统进行频率扫描，freqz(b, a, N) 中的 N 控制频率分辨率，默认输出为幅度和相位响应图。如果需要获取具体的响应数据，可接收 freqz 的输出并分别使用 abs 和 angle 分别处理其模和相角，再绘制自定义图像。该方法适用于分析离散 LTI 系统的频率特性。"
326 },
327 {
328 "Question": "已知 H(z) = (0.2 + 0.3z⁻¹ + z⁻²) / (1 + 0.8z⁻¹ + 0.15z⁻²)，用 residue 函数将 H(z) 展开为部分分式形式，画出零极点分布图，求系统的频率响应和单位脉冲响应，如何使用 MATLAB 实现？",
329 "Response": "将 H(z) 乘 z² 化为正幂形式，得到分子为 [0.2 0.3 1]，分母为 [1 0.8 0.15]；使用 residue(b, a) 展开部分分式；使用 zplane(b, a) 绘制零极点图；使用 freqz(b, a, N) 计算频率响应；使用 impz(b, a, N) 或 filter(b, a, delta) 求单位脉冲响应。",
330 "Complex_CoT": "原始 H(z) 表达式为负幂形式，不便于进行系数提取和函数调用，需通过两边同乘 z² 化为标准的正幂形式，得到分子为 0.2z² + 0.3z + 1，分母为 z² + 0.8z + 0.15。将其转换为系数向量 b = [0.2 0.3 1]，a = [1 0.8 0.15]。调用 residue(b, a) 得到部分分式展开，其中 r 为残数，p 为极点，k 为多项式部分。使用 zplane(b, a) 绘制零极点图以观察系统特性。频率响应可使用 freqz 函数计算，通过指定点数 N 得到在 [0, π] 上的频率响应，可配合 abs 和 angle 分析模和相位。单位脉冲响应可通过 impz(b, a, N) 获取前 N 项响应，或用 filter 对单位冲激序列 [1 0 0 ...] 进行滤波模拟系统响应。"
331 },
332 {
333 "Question": "用 MATLAB 实现信号 y = 5e^(-3t)·sin(4t - 1)·u(t)，应如何操作？",
334 "Response": "首先定义时间变量 t，例如使用 t = 0:dt:T，其中 dt 为时间间隔，T 为最大时间；使用 y = 5 * exp(-3*t) .* sin(4*t - 1) 构造信号表达式；为实现单位阶跃函数 u(t)，可用逻辑表达式 (t >= 0) 对信号 y 进行乘法屏蔽；最后使用 plot(t, y) 绘制波形。",
335 "Complex_CoT": "该信号为连续时间函数，包含指数衰减项、正弦项和单位阶跃函数。首先需建立时间向量 t，范围选在 t ≥ 0，例如 t = 0:0.01:5，确保阶跃函数定义区间正确。y 的表达式为 5 * exp(-3*t) * sin(4*t - 1)，但由于 MATLAB 中数组乘法需显式逐元素操作，需写成点乘形式：y = 5 * exp(-3*t) .* sin(4*t - 1)。为了实现 u(t)，应构造逻辑掩码 (t >= 0)，将其与 y 相乘，实现阶跃函数控制，即 y = y .* (t >= 0)。最后使用 plot 绘图函数显示结果，添加标题与标签用于标识信号含义。"
336 },
337 {
338 "Question": "用 MATLAB 实现信号 y = δ(t - 1)，应如何操作？",
339 "Response": "定义时间变量 t，例如 t = 0:dt:T，其中 dt 是采样间隔；使用 zeros(size(t)) 创建全零向量 y；通过 find(abs(t - 1) < ε) 找到最接近 t = 1 的索引位置，将该位置赋值为 1；使用 stem(t, y) 绘制离散冲激信号。",
340 "Complex_CoT": "δ(t - 1) 表示单位冲激函数在 t = 1 时刻发生，仅在该时刻取值为 1，其余为 0。在 MATLAB 中不能表示真正的冲激函数，因此采用离散近似方式。首先定义一个足够密集的时间向量 t，如 t = 0:0.01:2，以保证能精确接近 t = 1。然后创建与 t 同维度的零向量 y，用于初始化信号。为了定位 t = 1 所在位置，使用逻辑条件 abs(t - 1) < ε，其中 ε 是允许的误差阈值，如 1e-6，通过 find 函数找出对应索引，将该索引位置置为 1。最后使用 stem 函数绘制该离散冲激信号。此方法适用于仿真中对冲激函数的近似建模。"
341 },
342 {
343 "Question": "实现复指数信号 y = e^{(-0.5 + j8)t}，如何使用 MATLAB 实现？",
344 "Response": "定义时间变量 t，例如 t = 0:dt:T；使用 y = exp((-0.5 + 1i*8) * t) 构造复指数信号；使用 real(y) 和 imag(y) 提取实部与虚部，分别使用 plot 进行绘图；如需观察复平面轨迹，可使用 plot3 绘制三维图像。",
345 "Complex_CoT": "该信号为复指数形式，指数中包含实部 -0.5 和虚部 j8，表示一个幅度随时间指数衰减、相位线性变化的复信号。首先定义时间变量 t，需覆盖一定范围，例如从 0 到若干秒，并设置足够小的时间间隔 dt。使用 MATLAB 中 exp 函数和虚数单位 1i 构造信号表达式 y = exp((-0.5 + 1i*8) * t)。该信号是复值，绘图时需将 y 的实部和虚部分别提取，可使用 real(y) 和 imag(y)，再分别绘制。若希望观察其在复平面上的旋转轨迹，可以使用 plot3(t, real(y), imag(y)) 实现三维可视化，展示随时间演化的复信号行为。"
346 },
347 {
348 "Question": "用 MATLAB 实现信号 y = Sa[π(t - 2)]，-4 < t < 8，如何操作？",
349 "Response": "定义时间变量 t，例如使用 t = -4:dt:8；将表达式 y 写成 MATLAB 形式为 y = sinc(t - 2)，因为 MATLAB 中 sinc(x) 实现的是 sin(πx)/(πx)；使用 plot(t, y) 绘制信号波形。",
350 "Complex_CoT": "该信号为抽样函数形式，其中 Sa[π(t - 2)] 可表示为 sinc(t - 2)，原因是 MATLAB 中定义的 sinc 函数为归一化形式 sinc(x) = sin(πx)/(πx)，故将 π(t - 2) 对应为 sinc(t - 2)。首先定义时间变量 t，范围为 -4 到 8，选用适当的采样间隔确保图像光滑，例如 dt = 0.01。根据函数表达式，直接调用 y = sinc(t - 2) 构造信号向量。该函数在 t = 2 附近取最大值，其他区域按抽样函数特性衰减并振荡。最后使用 plot 绘图函数绘制波形，必要时可添加标签和网格提升可读性。"
351 },
352 {
353 "Question": "用 MATLAB 实现一个方波信号，应如何操作？",
354 "Response": "定义时间变量 t，例如 t = 0:dt:T；选择合适频率 f，使用 square(2*pi*f*t) 生成标准方波信号；如需调整占空比，可使用 square(2*pi*f*t, duty) 设置所需比例；使用 plot(t, y) 绘制波形。",
355 "Complex_CoT": "方波是一种周期性二值信号，在一个周期内高电平和低电平各占一定比例。MATLAB 中提供了 square 函数可直接生成该信号。首先定义时间向量 t，设置合适的时间步长 dt 和总时长 T，确保包含多个周期。然后设置频率 f，调用 square(2*pi*f*t) 构造方波；若需控制高电平持续时间，可通过 square 的第二个参数设置占空比，如 square(2*pi*f*t, 25) 表示高电平占 25%。生成信号后，使用 plot 函数将 t 与方波信号对应绘图，观察周期性结构和占空特性。"
356 },
357 {
358 "Question": "用 MATLAB 实现一个非周期型三角波脉冲信号，应如何操作？",
359 "Response": "先定义时间变量 t，例如 t = -1:dt:1；初始化 y 为零向量；设定脉冲中心为 t = 0，宽度为 T，分别在区间 t ∈ [-T, 0] 和 t ∈ [0, T] 内构造上升段和下降段，使用线性表达式实现对称三角波；使用逻辑索引对两个区间赋值；最后使用 plot(t, y) 绘图。",
360 "Complex_CoT": "非周期三角波脉冲信号具有单次上升和下降过程，中心对称且两侧为零。首先定义合适的时间范围，例如 t = -1:0.01:1，覆盖脉冲中心及两侧。初始化 y 为与 t 等长的零向量。设定脉冲中心 t0 = 0，脉宽为 T，如 T = 0.5，则信号在 [-0.5, 0] 内为上升段，y = (t - t1)/(t0 - t1)；在 [0, 0.5] 内为下降段，y = (t2 - t)/(t2 - t0)。用逻辑条件定位 t 所属区间，对应赋值构造信号。边界外值保持为零。最后用 plot 函数将信号绘出，观察脉冲形状是否对称合理。"
361 },
362 {
363 "Question": "用 MATLAB 实现一个周期型三角波脉冲信号，应如何操作？",
364 "Response": "定义时间变量 t，例如 t = 0:dt:T；设定频率 f，使用 sawtooth(2*pi*f*t, 0.5) 构造对称周期三角波；可对输出信号进行幅值缩放或偏移调整；使用 plot(t, y) 绘制信号波形。",
365 "Complex_CoT": "周期三角波是连续重复的线性上升和下降信号，周期性可通过正弦型基础函数构造。MATLAB 中可用 sawtooth 函数生成，形式为 sawtooth(2*pi*f*t, 0.5)，其中 f 为频率，0.5 表示对称三角波，占空比为 50%。首先应定义合适的时间向量 t，保证覆盖多个周期，常设为 t = 0:0.001:N。然后设定频率 f，结合 t 构造 sawtooth 表达式得到周期三角波。若需调整幅度，可对生成结果进行线性变换。最后通过 plot(t, y) 绘图观察周期性变化。该方法结构清晰，计算简洁，适用于周期性波形建模与可视化。"
366 },
367 {
368 "Question": "如何使用 MATLAB 实现信号 f(t) = (3 + e^(-t))·u(t)？",
369 "Response": "定义时间变量 t，例如 t = -2:dt:5；使用 y = (3 + exp(-t)) .* (t >= 0) 实现表达式，其中 (t >= 0) 表示单位阶跃函数；使用 plot(t, y) 绘图。",
370 "Complex_CoT": "该信号为指数衰减函数与常数相加后乘以单位阶跃函数，表示 t ≥ 0 时生效。首先需设置覆盖负值的时间向量以显示 u(t) 的截断效果。指数项通过 exp(-t) 实现，常数直接相加；为控制定义域，需通过逻辑条件 t >= 0 实现 u(t)，并与前面的函数做逐元素乘法。最后通过 plot 显示结果。"
371 },
372 {
373 "Question": "如何使用 MATLAB 实现信号 f(t) = [1 - 2cos(2πt)]·[u(t) - u(t - 2)]？",
374 "Response": "定义时间变量 t，例如 t = -1:dt:3；使用窗函数 (t >= 0) & (t < 2) 控制余弦项有效区间；构造表达式 y = (1 - 2*cos(2*pi*t)) .* ((t >= 0) & (t < 2))；使用 plot(t, y) 绘图。",
375 "Complex_CoT": "该信号在 [0,2) 内是余弦调制信号，其他时间为零，属于加窗信号形式。首先用逻辑表达式构造窗函数 [u(t) - u(t-2)]，即 (t >= 0) & (t < 2)。余弦项通过 cos(2*pi*t) 实现，对其乘上窗函数得到有效信号段。最后用 plot 绘图以确认其在 [0,2) 内振荡，其它处为零。"
376 },
377 {
378 "Question": "如何使用 MATLAB 实现信号 f(t) = Sa[π(t - 1)]？",
379 "Response": "定义时间变量 t，例如 t = -4:dt:6；使用 MATLAB 中的 sinc 函数表示 Sa 函数，表达式为 y = sinc(t - 1)；使用 plot(t, y) 绘图。",
380 "Complex_CoT": "Sa 表示抽样函数，在 MATLAB 中通过 sinc 函数实现，其中 sinc(x) 实际计算的是 sin(πx)/(πx)，因此 Sa[π(t - 1)] 等价于 sinc(t - 1)。需设置较大的时间范围和足够小的步长，以清晰显示主瓣和多个旁瓣。构造表达式 sinc(t - 1) 后直接绘图显示完整波形。"
381 },
382 {
383 "Question": "如何使用 MATLAB 实现信号 f(t) = 2 + 2e^(-πt/4)？",
384 "Response": "定义时间变量 t，例如 t = 0:dt:5；使用表达式 y = 2 + 2 * exp((-pi/4) * t) 构造信号；使用 plot(t, y) 绘制波形。",
385 "Complex_CoT": "该信号为一指数衰减项加上常数项的连续函数，无需分段处理。指数项由 exp((-pi/4)*t) 表示，整体加权并加常数即可。选择正向时间区间 t ≥ 0 足以展示衰减趋势，使用 plot 绘图可清晰看到从 4 衰减至 2 的过程。"
386 },
387 {
388 "Question": "如何使用 MATLAB 生成一个幅值为 1、宽度为 6 的非周期对称三角波信号并绘制波形图？",
389 "Response": "定义时间变量 t，例如 t = -5:dt:5；初始化 y 为零向量；设定有效区间为 t ∈ [-3, 3]，在区间 [-3, 0] 内定义上升段 y = (t + 3)/3，在区间 [0, 3] 内定义下降段 y = (3 - t)/3；使用逻辑索引分别赋值；使用 plot(t, y) 绘制三角波波形图。",
390 "Complex_CoT": "非周期对称三角波表示在一个有限区间内单次上升和下降，中心对称，题中幅值为 1、宽度为 6，因此中心在 t = 0，有效区间为 [-3, 3]。为了实现该波形，首先设定足够密集的时间向量 t 覆盖三角波范围。接着创建与 t 等长的零向量 y，作为信号初始化。然后用逻辑条件 t >= -3 & t <= 0 定义上升段，令 y = (t + 3)/3；用逻辑条件 t > 0 & t <= 3 定义下降段，令 y = (3 - t)/3。两段构造后拼接即得完整三角波。最后使用 plot 函数绘制波形图以观察其形状。"
391 },
392 {
393 "Question": "如何使用 MATLAB 生成信号 x₁ = u(t - 2) 和 x₂ = cos(10t)，并绘制 x₁、x₂ 及 x₁ + x₂ 的波形图？",
394 "Response": "定义时间变量 t，例如 t = 0:dt:T；使用逻辑表达式 (t >= 2) 构造 x₁；使用 cos(10*t) 构造 x₂；令 x3 = x₁ + x₂ 得到叠加信号；使用 subplot 分别绘制 x₁、x₂ 和 x₃ 的波形图。",
395 "Complex_CoT": "x₁ 是延迟2单位的单位阶跃函数，可通过逻辑表达式 (t >= 2) 构造布尔型向量表示；x₂ 是频率为10 rad/s 的余弦信号，通过 cos(10*t) 生成；需设定足够宽的时间区间以完整展现两个信号的行为，如 t = 0:0.01:5；x₁ 与 x₂ 元素一一对应，可直接进行向量加法得到叠加信号 x₃ = x₁ + x₂；用 subplot 分为三幅子图分别绘制 x₁、x₂ 和 x₃，确保观察各信号时域特性及其叠加效果。"
396 },
397 {
398 "Question": "已知连续信号 f(t) = sin(πt)，分别求其延时 0.1s 后的波形、翻转后的波形和频率增加为原来 2 倍的波形，并用 MATLAB 实现。",
399 "Response": "本题涉及信号的三种基本时域变换：延时、翻转和频率扩展。MATLAB 解决流程如下：1）构造时间轴 t，如 t = -1:0.001:1；2）计算原始信号 f = sin(pi * t)；3）延时操作通过替换 t 为 t - 0.1 得到延时信号；4）翻转操作替换 t 为 -t 得到翻转信号；5）频率加倍操作将 t 替换为 2*t 得到频率加倍信号；6）使用 plot 或 subplot 函数绘图展示各波形并加上标签便于比较。",
400 "Complex_CoT": "此问题考查学生对连续时间信号的基本变换理解。在 MATLAB 中，需要先构造一个足够精细的时间向量 t 模拟连续时间。原始信号通过向量计算得到。延时是 f(t - t0) 形式，只需对时间轴减去对应延时时间实现右移；时间翻转为 f(-t)，对应时间轴关于原点对称翻转；频率加倍对应 f(at)，a > 1 表示压缩时间轴，在此为 f(2t)，周期变短，频率加倍。每种变换均通过修改 t 并重新计算 f(t) 实现。绘图部分可用 subplot 实现多图并列，通过对比观察波形变化，加深理解时间变换对信号形态的影响。"
401 },
402 {
403 "Question": "求信号 y = sin(2t) 的微分和积分，并用 MATLAB 实现。",
404 "Response": "本题可通过 MATLAB 的符号计算或数值方法两种途径实现。符号法中，先定义符号变量 t 和函数 y = sin(2*t)，使用 diff(y) 求导，int(y) 求积分；数值法中，先定义时间向量 t 和函数 y = sin(2*t)，使用 diff(y)./diff(t) 实现导数近似，使用 cumtrapz(t, y) 实现积分近似。最后通过 plot 或 fplot 函数绘图展示导数与积分结果。",
405 "Complex_CoT": "首先明确 y = sin(2t) 是连续时间信号，理论导数为 y' = 2cos(2t)，积分为 -0.5cos(2t)。在 MATLAB 中，若使用符号方法，应使用 syms t 定义符号变量，y = sin(2*t) 定义函数，使用 diff(y) 得到导数，int(y) 得到积分，适用于解析解推导；若使用数值方法，应构造时间向量 t = a:dt:b，定义 y = sin(2*t)，用 diff(y)./diff(t) 得到数值导数，此时因长度变化需补点或重构时间轴；积分采用 cumtrapz(t, y) 完成梯形近似积分，适用于数据处理场景。无论哪种方法，最终使用 plot 或 fplot 将原函数、导数、积分结果绘制出来进行直观验证与分析。"
406 },
407 {
408 "Question": "已知信号 f₁(t) = (-t + 4)[u(t) - u(t - 4)], f₂(t) = sin(πt)，试用 MATLAB 绘制出 f₁(-t + 2) 的时域波形图。",
409 "Response": "此为时间反转加右移操作，先对原始信号时间轴 t 取负，再加 2。MATLAB 中令 t1 = -t + 2，再将 f₁ 按照 t1 计算。",
410 "Complex_CoT": "信号 f₁(t) 定义在 [0,4]，其结构为线性乘以单位阶跃差。f₁(-t+2) 等效于先反转，再向右平移 2 单位，即 t → -t+2。画图时应定义足够范围的时间变量 t，例如 -10 到 10，计算变换后的 t1 = -t + 2，然后计算 f₁(t1)。为了保证单位阶跃函数正确运行，应在 t1 的基础上计算 u(t1) 和 u(t1 - 4)。"
411 },
412 {
413 "Question": "已知信号 f₁(t) = (-t + 4)[u(t) - u(t - 4)], f₂(t) = sin(πt)，试用 MATLAB 绘制出 f₁(2t) 的时域波形图。",
414 "Response": "该信号是时间压缩 1/2 倍，即波形加快，宽度缩小。MATLAB 中使用 t1 = 2*t，按此替代变量绘图。",
415 "Complex_CoT": "f₁(2t) 表示时间轴压缩，原本定义在 [0,4] 的信号变为 [0,2]。画图时需设置原始时间向量 t，然后构造 t1 = 2*t，代入原始表达式 f₁(t1)。因为缩放后波形更窄，应适当放宽 t 的定义范围保证完整波形显示。单位阶跃函数应作用在 t1 上。"
416 },
417 {
418 "Question": "已知信号 f₁(t) = (-t + 4)[u(t) - u(t - 4)], f₂(t) = sin(πt)，试用 MATLAB 绘制出 f₁(½t - 2) 的时域波形图。",
419 "Response": "该为先时间拉伸再右移，MATLAB 中定义 t1 = 0.5*t - 2，将其代入原始 f₁ 表达式。",
420 "Complex_CoT": "变换 t → 0.5t - 2 包含两个操作：先放慢（拉伸）时间轴 2 倍，再右移 2 单位，即波形变宽并向右移动。在 MATLAB 中设定 t 后计算 t1 = 0.5*t - 2，将其代入 f₁ 表达式中，并根据 t1 定义单位阶跃函数以限制支持区间。应确保 t 的取值范围足够大以容纳变换后波形。"
421 },
422 {
423 "Question": "已知信号 f₁(t) = (-t + 4)[u(t) - u(t - 4)], f₂(t) = sin(πt)，试用 MATLAB 绘制出 d/dt₁[f₁(½t + 2)] 的时域波形图。",
424 "Response": "先进行时间变换 t1 = 0.5*t + 2，再计算该信号对 t 的导数。MATLAB 可使用 gradient(y, dt) 或 diff(y)/dt 求数值导数。",
425 "Complex_CoT": "首先对 t 进行变换，构造 t1 = 0.5*t + 2，代入原始信号得到 f₁(t1)。由于为离散信号，无法使用解析导数，采用数值方法：计算 f₁(t1) 后，使用 gradient 或 diff 对结果向量 y 进行差分，注意使用合适的时间步长 dt 以确保导数精度。最终使用 plot 绘图即可。"
426 },
427 {
428 "Question": "已知信号 f₁(t) = (-t + 4)[u(t) - u(t - 4)], f₂(t) = sin(πt)，试用 MATLAB 绘制出 ∫_{-∞}^{t} f(2 - τ) dτ 的时域波形图。",
429 "Response": "此为信号时移后积分形式，可理解为变量替换，实际计算为 t 轴上 f(2 - t) 的积分。使用 MATLAB 中的 cumtrapz 函数。",
430 "Complex_CoT": "积分上限为 t，变量替换后为 f(2 - τ)，表示将原始信号反转并右移 2 单位。在 MATLAB 中先构建时间轴 t，再构造反转移位后的 t1 = 2 - t，代入原始信号 f₁(t1)，随后对该向量使用 cumtrapz(t, f₁(2 - t)) 执行累计积分，模拟从 -∞ 到当前时间的积分。由于 t 是有限的，需设置足够宽的时间范围以近似逼近无穷积分。"
431 },
432 {
433 "Question": "已知信号 f₁(t) = (-t + 4)[u(t) - u(t - 4)], f₂(t) = sin(πt)，试用 MATLAB 绘制出 f₁(-t) + f₂(t) 的时域波形图。",
434 "Response": "先对 f₁ 做时间反转，再与 f₂ 相加。MATLAB 中定义 t1 = -t，计算 f₁(t1) 和 f₂(t)，元素对应相加。",
435 "Complex_CoT": "f₁(-t) 表示时间反转，对 t 轴取负构造新时间变量 t1 = -t，带入 f₁ 表达式计算。f₂(t) 为正弦函数无需变换。将两个结果逐点相加，MATLAB 中用向量操作 f1_neg + f2 完成叠加。注意两个信号定义域必须一致，即 t 向量一致，确保对应元素相加无误。"
436 },
437 {
438 "Question": "已知信号 f₁(t) = (-t + 4)[u(t) - u(t - 4)], f₂(t) = sin(πt)，试用 MATLAB 绘制出 -[f₁(t) + f₂(t)] 的时域波形图。",
439 "Response": "先按原始时间轴计算 f₁ 和 f₂，再将两者相加后取反。MATLAB 中使用负号操作即可。",
440 "Complex_CoT": "此为直接信号叠加并取负。f₁(t) 与 f₂(t) 均使用同一时间轴计算，先求和向量 f_sum = f1 + f2，再取其负值 f_neg = -f_sum。该操作属于幅度反转，MATLAB 可直接进行。确保两个信号采样时间一致。"
441 },
442 {
443 "Question": "已知信号 f₁(t) = (-t + 4)[u(t) - u(t - 4)], f₂(t) = sin(πt)，试用 MATLAB 绘制出 f₁(t) × f₂(t) 的时域波形图。",
444 "Response": "逐点相乘，使用 MATLAB 中的点乘操作符 .* 实现。",
445 "Complex_CoT": "此为两个信号的乘积，需对相同时间轴上的 f₁ 和 f₂ 逐点相乘，表示信号调制或加窗操作。定义时间轴后分别计算 f₁(t) 与 f₂(t)，再执行 element-wise 乘法操作 f_prod = f1 .* f2。画图观察相乘后波形特性变化，常用于调幅等信号处理操作。"
446 },
447 {
448 "Question": "求信号 y₁(t) = sin(t) 和 y₁(t) = e^(-2t) 的卷积。",
449 "Response": "在 MATLAB 中处理连续信号卷积，需先定义合理的时间采样区间，使用 linspace 或 0:dt:T 构造时间轴。生成 sin(t) 和 e^(-2t) 的离散信号序列，使用 sin(t) 与 exp(-2*t) 计算。使用 conv 函数进行离散卷积，同时乘以采样间隔 dt 实现对连续卷积的数值逼近。注意卷积后时间轴变化，应通过 t_conv = t1(1)+t2(1):dt:t1(end)+t2(end) 构造新时间轴，并用 plot 函数可视化卷积结果。",
450 "Complex_CoT": "首先分析两个信号：y₁(t) = sin(t) 为正弦波，y₂(t) = e^(-2t) 为指数衰减函数，卷积需满足信号因果性，因此需在 t ≥ 0 上定义。由于 MATLAB 中无法直接对连续信号卷积，需转化为数值计算，即通过高分辨率离散采样逼近积分。第一步定义时间向量 t，选择合适的时间范围 [0, T] 和小的 dt 确保精度，用 0:dt:T 构造。第二步生成信号序列 y1 = sin(t)，y2 = exp(-2*t)。第三步使用 conv(y1, y2)*dt 计算卷积，乘以 dt 是因为数值卷积本质是求和逼近积分。卷积结果长度为 length(y1)+length(y2)-1，因此需通过 t_conv = t(1)+t(1):dt:t(end)+t(end) 构建对应时间轴。最后用 plot(t_conv, y_conv) 可视化卷积结果，观察信号相互作用后的输出形状。"
451 },
452 {
453 "Question": "试用 MATLAB 符号运算法求 y₁(t) = t 和 y₂(t) = sin(t) 的卷积。",
454 "Response": "本题需使用 MATLAB 符号计算工具箱进行解析卷积运算。首先使用 syms 定义符号变量 t 和 tau。将 y₁(t) 表示为 tau，y₂(t) 需变量替换为 sin(t - tau)，用 subs 实现替换。接着利用 int 函数对表达式 tau * sin(t - tau) 在 tau ∈ [0, t] 上积分，得到符号卷积结果。最后可用 simplify 简化表达式，用 fplot 进行可视化。",
455 "Complex_CoT": "卷积的符号定义为 y(t) = ∫₀ᵗ y₁(τ)·y₂(t−τ) dτ，因此需先引入符号变量 t 和积分变量 tau，分别用于输出变量和中间积分变量。y₁(t) = t 在积分中写作 tau，y₂(t) = sin(t) 替换为 sin(t − tau)，通过 subs 函数实现变量替换。然后构造卷积表达式：int(tau * sin(t - tau), tau, 0, t)，得到符号形式的卷积表达式。为便于后续分析，可用 simplify 对表达式进行简化，用 fplot(expr, [a, b]) 对符号函数进行图形展示。符号卷积适合解析处理，尤其当信号具备明确表达式时可得出准确结果。"
456 },
457 {
458 "Question": "用 MATLAB 命令绘制 f₁(t) = u(t), f₂(t) = u(t) - u(t - 3) 的卷积积分 f₁(t) * f₂(t) 的时域波形图。",
459 "Response": "该题为阶跃函数与矩形脉冲的卷积。需定义时间轴 t，构造 f₁ 为 t ≥ 0，f₂ 为 t ≥ 0 减去 t ≥ 3。使用 conv(f1, f2) * dt 进行离散卷积，并构造新时间轴用于绘图。用 plot 进行可视化。",
460 "Complex_CoT": "设定 t ∈ [-1, 10]，选择合适 dt 离散时间。f₁(t) 为单位阶跃函数，可用 f1 = (t >= 0) 表示。f₂(t) 为一个宽度为 3 的矩形脉冲，用 f2 = (t >= 0) - (t >= 3) 表示。使用 conv 函数对两个离散信号进行卷积并乘以 dt 实现近似积分。由于卷积结果长度为 length(f1) + length(f2) - 1，需用 t_conv = t(1)+t(1):dt:t(end)+t(end) 生成新时间轴。最后用 plot(t_conv, y) 绘制卷积结果时域图像。"
461 },
462 {
463 "Question": "用 MATLAB 命令绘制 f₁(t) = u(t), f₂(t) = δ(t) 的卷积积分 f₁(t) * f₂(t) 的时域波形图。",
464 "Response": "该题为单位阶跃函数与冲激函数的卷积。f₁ 为 t ≥ 0，f₂ 近似为在 t=0 处为1，其余为0的向量。使用 conv(f1, f2) * dt 后应得原信号 f₁(t)。构造新时间轴并绘制波形。",
465 "Complex_CoT": "定义时间轴 t ∈ [-1, 5]，dt 足够小。f₁ 为单位阶跃函数 f1 = (t >= 0)，f₂ 为单位冲激函数近似 delta = (t == 0)。由于 δ(t) 是卷积单位元，卷积后应保留原信号形状。执行 conv(f1, delta) * dt，结果为 f₁(t)。构建新时间轴 t_conv 并用 plot 绘图验证结果。此题目旨在验证卷积单位元的性质。"
466 },
467 {
468 "Question": "验证 y₁(t) = sin(t) 与 y₂(t) = u(t) 的卷积满足交换律，即 y₁(t) * y₂(t) = y₂(t) * y₁(t)。",
469 "Response": "在 MATLAB 中使用统一时间轴构造 sin(t) 与单位阶跃函数 u(t)，分别计算两种卷积顺序的结果：conv(y1, y2)*dt 与 conv(y2, y1)*dt。构造统一的新时间轴，对两个结果进行绘图对比或求差异验证其是否完全相等，从而验证交换律成立。",
470 "Complex_CoT": "首先建立时间轴 t = 0:dt:T，dt 取较小值确保精度。构造 y1 = sin(t)，y2 = (t >= 0) 表示单位阶跃函数。使用 conv(y1, y2)*dt 和 conv(y2, y1)*dt 计算两种顺序的卷积结果，卷积后长度为 length(y1)+length(y2)-1。构造新时间轴 t_conv = t(1)+t(1):dt:t(end)+t(end)。将两个卷积结果绘图比较，或通过 max(abs(a - b)) 判断是否全零以验证交换律是否成立。"
471 },
472 {
473 "Question": "验证 y₁(t) = sin(t)、y₂(t) = u(t)、y₃(t) = exp(-t) 满足卷积的分配律：y₁(t) * [y₂(t) + y₃(t)] = y₁(t) * y₂(t) + y₁(t) * y₃(t)。",
474 "Response": "在 MATLAB 中构造 sin(t)、u(t)、exp(-t) 三个信号。先将 u(t) 与 exp(-t) 相加形成组合信号，再与 sin(t) 卷积；另分别对 sin(t) 与 u(t)、exp(-t) 卷积后再相加。对比两种结果的差异或绘图验证分配律是否成立。",
475 "Complex_CoT": "设定统一时间轴 t = 0:dt:T。构造 y1 = sin(t)，y2 = (t >= 0)，y3 = exp(-t)。计算左式卷积：y_sum = y2 + y3，再执行 conv(y1, y_sum)*dt。再分别计算 conv(y1, y2)*dt 与 conv(y1, y3)*dt，并相加得到右式结果。卷积结果长度均为 length(y1)+length(y2)-1，构造统一 t_conv 作为新时间轴。绘图或比较两组卷积值的差异，验证分配律是否成立。"
476 },
477 {
478 "Question": "已知某LTI系统微分方程为 y'''(t) + 3y''(t) + 2y'(t) = 0，初始条件为 y(0⁻) = 2，y'(0⁻) = 1，y''(0⁻) = 2，求其零输入响应。",
479 "Response": "可使用MATLAB符号解函数dsolve解析求解。在命令窗口中定义符号变量y(t)，用diff表示导数，写出方程diff(y,t,3)+3*diff(y,t,2)+2*diff(y,t)==0，并同时设置初始条件y(0)==2，Dy(0)==1，D2y(0)==2。将方程和初始条件作为参数传入dsolve中，即可直接获得解析解。",
480 "Complex_CoT": "dsolve适用于求解具有解析解的常微分方程，特别是线性常系数微分方程。首先定义符号变量 y(t)，用 diff(y,t,n) 分别表示 y 的 n 阶导数。将原方程表示为 diff(y,t,3) + 3*diff(y,t,2) + 2*diff(y,t) == 0。接着写出初始条件 y(0)==2，Dy(0)==1，D2y(0)==2，其中 Dy 和 D2y 是一阶和二阶导数的符号变量，需用 syms y(t) 后自动生成。使用 dsolve([方程], [初始条件]) 调用求解器，返回的即为 y(t) 的显式表达式。该方法优于数值法的优势在于能直接获取解析公式，适用于理论分析或验证。"
481 },
482 {
483 "Question": "已知某LTI系统微分方程为 y''(t) + 2y'(t) + y(t) = 2x(t)，输入为 x(t) = e^(-2t)u(t)，求系统的零状态响应。",
484 "Response": "本题建议使用MATLAB传递函数建模方式，先根据微分方程写出系统的传递函数 H(s) = Y(s)/X(s) = 2 / (s^2 + 2s + 1)。可通过tf函数构造系统模型，其中分子为2，分母为[1 2 1]。再构造输入信号x(t)=e^(–2t)u(t)对应的向量形式，并定义时间轴t，利用lsim函数计算系统对该输入的零状态响应，即 y = lsim(sys, x, t)，其中sys为系统传递函数模型，x为输入序列，t为时间向量。",
485 "Complex_CoT": "题目给出线性常系数非齐次微分方程及具体输入，要求求解系统的零状态响应，表示系统在初始状态为零时对外部输入的响应。由于输入为e^(–2t)u(t)，属于解析形式，适合频域建模。首先对微分方程进行拉普拉斯变换，得到 Y(s)(s^2 + 2s + 1) = 2X(s)，即系统传递函数为 H(s) = 2 / (s+1)^2。该系统可用MATLAB中的tf函数构建传递函数模型。为了模拟系统响应，需要将x(t)=e^(–2t)u(t)在离散时间上表示为向量输入，结合时间向量t构造完整输入序列。利用lsim函数对系统进行线性仿真，输入为sys（系统模型）、x（输入信号向量）和t（时间轴），最终得到y(t)的数值解并可视化分析。若需更精确模拟单位阶跃函数，可使用heaviside或手动限定e^(–2t)在t≥0内定义。该方法适用于处理任意解析输入与线性系统之间的零状态响应分析。"
486 },
487 {
488 "Question": "已知某LTI系统微分方程为 y''(t) + 2y'(t) + y(t) = 2x(t)，输入为 x(t) = e^(-2t)u(t)，试用MATLAB符号法求系统的零状态响应。",
489 "Response": "此题可使用MATLAB符号工具箱中的dsolve结合laplace与ilaplace函数求解。首先定义符号变量y(t)、x(t)，并将微分方程 y'' + 2y' + y = 2x 表示为符号表达式。对该方程进行拉普拉斯变换，设初始条件为零（零状态响应），利用laplace函数分别变换y和x项，转换为代数方程，解出Y(s)。再使用ilaplace对Y(s)进行反变换，得到y(t)的解析表达形式。",
490 "Complex_CoT": "该题求解目标是线性微分方程在零初始条件下的零状态响应，且给定输入x(t)=e^(–2t)u(t)，适合通过符号法使用拉普拉斯变换处理。首先在MATLAB中定义符号变量t、s、y(t)、x(t)，并设定x(t)=exp(–2*t)。然后利用laplace函数分别求出y''、y'、y的拉普拉斯变换，注意设置初始条件为0，得到Y(s)(s^2+2s+1)=2/(s+2)。整理得Y(s)=2/[(s+2)(s+1)^2]，此为频域响应表达式。接着使用ilaplace函数对Y(s)反变换回时域，得到y(t)的解析表达式，完成符号求解。该方法优势在于直接获得闭式解，适合进行符号化分析或理论验证，核心在于精确建立拉普拉斯代数表达及其求解流程。"
491 },
492 {
493 "Question": "已知某LTI系统的微分方程为：y''(t) + 5y'(t) + 3y(t) = x'(t) + 3x(t)，求该系统的冲激响应和阶跃响应。",
494 "Response": "将系统微分方程两边对t作拉普拉斯变换，假设初始条件为零，左边为Y(s)(s² + 5s + 3)，右边为X(s)(s + 3)，得系统传递函数H(s) = Y(s)/X(s) = (s + 3)/(s² + 5s + 3)。用MATLAB中tf函数构造传递函数：分子系数为[1 3]，分母系数为[1 5 3]，用impulse和step函数分别求冲激响应和阶跃响应，配合时间向量绘图即可。",
495 "Complex_CoT": "首先识别该系统为线性时不变系统，可用拉普拉斯变换将微分方程转化为代数方程。在零初始条件下，y''(t) → s²Y(s)，y'(t) → sY(s)，x'(t) → sX(s)，代入得：Y(s)(s² + 5s + 3) = X(s)(s + 3)，从而得传递函数H(s) = (s + 3)/(s² + 5s + 3)。该传递函数形式可直接用MATLAB的tf函数构造，其中num = [1 3]，den = [1 5 3]。定义好时间轴t后，使用impulse(sys, t)可得到冲激响应，step(sys, t)得到阶跃响应。为增强可视化，可使用plot对响应进行标注和分析，辅助判断系统响应特性如稳定性、过渡过程等。"
496 },
497 {
498 "Question": "已知某连续LTI的微分方程为：y''(t) + 3y'(t) + 10y(t) = x'(t) + 4x(t)，其中x(t) = e^{-3t}u(t)。试利用MATLAB卷积积分法求解该系统零状态响应。",
499 "Response": "对微分方程两边进行拉普拉斯变换并假设初始条件为零，得传递函数H(s) = (s + 4)/(s² + 3s + 10)。用MATLAB中tf函数构造传递函数，num为[1 4]，den为[1 3 10]，用impulse函数求出冲激响应h(t)。输入x(t) = e^{-3t}u(t)，可用exp和逻辑索引生成离散表达。使用conv函数对h(t)和x(t)进行卷积，乘以时间间隔dt近似连续卷积结果，得到零状态响应y(t)。也可直接用lsim函数求系统对x(t)的响应。",
500 "Complex_CoT": "首先确认题设系统为线性时不变系统，适合用卷积积分法求零状态响应。将微分方程两边拉普拉斯变换，y'' → s²Y(s)，y' → sY(s)，x' → sX(s)，得Y(s)(s² + 3s + 10) = X(s)(s + 4)，从而H(s) = (s + 4)/(s² + 3s + 10)。用MATLAB中tf函数输入分子[1 4]、分母[1 3 10]生成系统sys。使用impulse(sys, t)求出系统单位冲激响应h(t)。构造时间向量t，并生成x(t) = exp(-3*t) .* (t >= 0)。用conv(h, x)*dt实现卷积运算，补齐时间轴后绘图显示y(t)。若使用lsim(sys, x, t)可直接得到响应而无需手动卷积。"
501 },
502 {
503 "Question": "已知某连续LTI的微分方程为：y''(t) + 4y'(t) + 3y(t) = 0，试用MATLAB求起始条件为y(0₊) = 4，y'(0₊) = 2的零输入响应。",
504 "Response": "将二阶微分方程转化为状态空间形式，设x₁ = y，x₂ = y'，得到一阶状态方程dx/dt = A·x，其中A为系统矩阵。构建状态空间模型后，利用MATLAB的ss函数建立系统对象，再使用initial函数结合初始状态向量x0 = [4; 2]和时间向量t，求系统在零输入下的响应。",
505 "Complex_CoT": "首先将微分方程y'' + 4y' + 3y = 0重写为状态空间形式。设x₁ = y，x₂ = y'，则有dx₁/dt = x₂，dx₂/dt = -4x₂ - 3x₁，写成矩阵形式为dx/dt = A·x，其中A = [[0 1]; [-3 -4]]。输出y(t) = x₁，对应C = [1 0]，D = 0。构造状态空间系统sys = ss(A,B,C,D)，由于为零输入响应，输入为零，即B为零向量。然后设置初始状态x0 = [4; 2]，定义时间向量t = 0:dt:T，使用initial(sys, x0, t)计算y(t)。最后用plot绘图分析响应特性，确认系统是否收敛、是否存在振荡等行为。"
506 },
507 {
508 "Question": "利用MATLAB求解系统 y''(t) + y'(t) + 7y(t) = 2x'(t) + x(t) 在冲激、阶跃、正弦输入下的零状态响应。",
509 "Response": "将微分方程两边进行拉普拉斯变换，初始条件为零，得到系统传递函数 H(s) = (2s + 1)/(s² + s + 7)。用 tf 函数构造该传递函数，分子为 [2 1]，分母为 [1 1 7]。利用 impulse 求冲激响应，step 求阶跃响应，构造 sin(ωt) 用 lsim 求正弦响应，响应结果用 plot 绘图观察。",
510 "Complex_CoT": "将 y''(t) + y'(t) + 7y(t) = 2x'(t) + x(t) 进行拉普拉斯变换，得 Y(s)(s² + s + 7) = X(s)(2s + 1)，从而得到传递函数 H(s) = (2s + 1)/(s² + s + 7)。在MATLAB中用 tf([2 1], [1 1 7]) 构造传递函数 sys。调用 impulse(sys, t) 得到冲激响应，step(sys, t) 得到阶跃响应，正弦输入如 x(t) = sin(3t)，构造输入向量后用 lsim(sys, x, t) 求响应。时间向量用 linspace 或 0:dt:T 构造，绘图时用 plot(t, y)，结合图像观察系统动态特性。"
511 },
512 {
513 "Question": "利用MATLAB求解系统 y'''(t) + 7y''(t) + 6y'(t) + 32y(t) = x(t) 在冲激、阶跃、正弦输入下的零状态响应。",
514 "Response": "将微分方程两边进行拉普拉斯变换，初始条件为零，得到系统传递函数 H(s) = 1/(s³ + 7s² + 6s + 32)。用 tf 函数构造传递函数，分子为 1，分母为 [1 7 6 32]。利用 impulse 求冲激响应，step 求阶跃响应，构造正弦信号用 lsim 求正弦响应，配合 plot 绘图分析。",
515 "Complex_CoT": "对方程 y'''(t) + 7y''(t) + 6y'(t) + 32y(t) = x(t) 进行拉普拉斯变换，在零初始条件下有 Y(s)(s³ + 7s² + 6s + 32) = X(s)，得 H(s) = 1/(s³ + 7s² + 6s + 32)。在MATLAB中使用 tf([1], [1 7 6 32]) 构造传递函数 sys，分别用 impulse(sys, t)、step(sys, t) 计算冲激和阶跃响应。正弦输入如 x(t) = sin(5t)，构造时间向量和输入信号后使用 lsim(sys, x, t) 求正弦响应。所有响应用 plot(t, y) 绘图以分析系统对不同输入的动态行为。"
516 },
517 {
518 "Question": "试用 MATLAB 求信号 x = e^{-t^2} 的傅里叶变换。",
519 "Response": "定义符号变量 t 和 w，使用 syms t w。定义信号表达式 x = exp(-t^2)。调用 fourier(x, t, w) 计算傅里叶变换，结果以 w 为频率变量的解析表达式输出。可使用 simplify 对结果化简，使用 fplot 绘制频域表达式随 w 变化的图像。",
520 "Complex_CoT": "傅里叶变换要求基于积分定义对时间函数进行变换，信号 x = exp(-t^2) 是连续实函数，傅里叶变换定义为 ∫x(t)·exp(-j·w·t) dt。使用 syms 定义符号变量 t 和频率变量 w，建立表达式 x = exp(-t^2)。使用 fourier(x, t, w) 调用 MATLAB 内建的傅里叶变换符号运算，返回解析表达式。变换结果通常包含指数形式，可进一步用 simplify 精简结构。若需查看频域分布，可用 fplot 对变换结果在频率轴上进行绘图，分析其谱集中性与衰减特性。"
521 },
522 {
523 "Question": "试用 MATLAB 求 X(jω) = 3 / (5 + jω) 的傅里叶逆变换。",
524 "Response": "定义符号变量 w 和 t，使用 syms w t。用 1i 表示虚数单位，将频域表达式定义为 X = 3 / (5 + 1i*w)。调用 ifourier(X, w, t) 计算逆变换，结果为时域表达式 x(t)。可用 simplify 对结果化简，用 fplot 绘图观察时域特性。",
525 "Complex_CoT": "目标是从频域表达式 X(jω) = 3 / (5 + jω) 还原出其对应的时域信号 x(t)。首先在 MATLAB 中用 syms 定义符号变量 w 和 t，确保后续表达式为符号形式处理。定义频域函数 X = 3 / (5 + 1i*w)，注意 MATLAB 中复数单位必须为 1i 而非 j。调用 ifourier(X, w, t) 使用傅里叶逆变换积分公式 ∫X(jω)·exp(jωt) dω，输出对应的 x(t)。该结果为解析表达式，可进一步用 simplify 精简其符号结构。若需观察该时域信号的行为，可使用 fplot 对 x(t) 在 t 的指定区间进行绘图。整个过程中需注意频域表达式的形式规范性与变量定义的准确性，避免数值与符号混用导致变换失败。"
526 },
527 {
528 "Question": "试用 MATLAB 绘制单边指数信号 x(t) = e^{-3t}·u(t) 的幅度谱和相位谱。",
529 "Response": "定义符号变量 t 和 w，构造信号 x = exp(-3*t) * heaviside(t)。使用 fourier(x, t, w) 计算频域表达式 X。用 abs(X) 得到幅度谱，用 angle(X) 得到相位谱。用 matlabFunction 将符号表达式转为函数句柄，生成频率向量 w = linspace(...)，用 fplot 或 plot 绘制幅度谱和相位谱。",
530 "Complex_CoT": "信号为实数单边指数信号，先定义符号变量 t 和 w。表达式构造为 x = exp(-3*t) * heaviside(t)，使用 fourier(x, t, w) 得到解析频谱 X(w)。X 为复数函数，用 abs 计算幅度谱，用 angle 计算相位谱。将 X 转换为数值函数句柄，使用 matlabFunction。设定频率范围，例如 w = linspace(-20, 20, 1000)，分别绘制 abs(X(w)) 和 angle(X(w))。如相位跳变较大，可用 unwrap 处理相位连续性，最终分别绘制幅度谱和相位谱。"
531 },
532 {
533 "Question": "试用 MATLAB 数值法求单边指数信号 x(t) = e^{-3t}·u(t) 的幅度谱。",
534 "Response": "设定时间范围 t = 0:dt:T，T 取足够大以确保信号充分衰减。构造信号 x = exp(-3*t)。使用 fft(x) 计算离散傅里叶变换，用 fftshift 调整频谱顺序。计算采样频率 Fs = 1/dt，采样点数 N = length(x)，构造频率轴 f = (-N/2:N/2-1)*(Fs/N)。使用 abs 提取幅度谱，乘以 dt 进行归一化，用 plot(f, abs(X)) 绘制幅度谱图像。",
535 "Complex_CoT": "信号为单边实指数信号，需在有限时间区间内数值计算其频谱。首先设定采样时间范围 t = 0:dt:T，T 应足够大以保证 exp(-3t) 衰减至接近零。定义信号 x = exp(-3*t)，此时无需单独乘以阶跃函数。使用 fft(x) 计算其离散傅里叶变换，结果为复数向量。调用 fftshift 对频谱进行居中处理，使零频率位于中间位置。计算采样频率 Fs = 1/dt，采样点数 N = length(x)，用公式 f = (-N/2:N/2-1)*(Fs/N) 构造对应频率轴。用 abs(X) 得到幅度谱，为近似连续傅里叶变换，还需乘以 dt 归一化。最后使用 plot 绘制频率与幅度的关系图。"
536 },
537 {
538 "Question": "假设 x(t) = 2e^{-2t}·u(t)，y(t) = 2e^{-2(t-1)}·u(t-1)，试用 MATLAB 绘制出 x(t) 和 y(t) 的幅度谱和相位谱，并比较时域平移在频域中相位的变化。",
539 "Response": "定义符号变量 t 和 w，构造信号 x = 2*exp(-2*t)*heaviside(t)，y = 2*exp(-2*(t-1))*heaviside(t-1)。分别对 x 和 y 使用 fourier 计算频域表达式 X 和 Y。使用 abs 和 angle 分别计算两者的幅度谱与相位谱。使用 fplot 绘制 X 与 Y 的幅度谱和相位谱，对比观察相位变化。",
540 "Complex_CoT": "x(t) 和 y(t) 是一对具有时间平移关系的单边指数信号，傅里叶变换中时移对应频域乘以 e^{-jωτ}，引起线性相位变化。使用 syms 定义符号变量 t 和 w，构造 x = 2*exp(-2*t)*heaviside(t)，y = 2*exp(-2*(t-1))*heaviside(t-1)。调用 fourier(x, t, w) 和 fourier(y, t, w) 分别求得 X(w) 和 Y(w)。使用 abs(X)、angle(X)、abs(Y)、angle(Y) 分别求得两信号的幅度谱与相位谱。使用 fplot 绘制频率范围内的谱图，并将 X 与 Y 的相位谱进行比较，可观察到 Y(w) 比 X(w) 多出一个线性相位项 –ω，验证傅里叶变换的时移性质。"
541 },
542 {
543 "Question": "设目信号为 x(t) = u(t+1) – u(t–1)，用它对载频信号 cos(ω₀t) 进行调幅，得到 y(t) = cos(ω₀t)·x(t)。试用 MATLAB 计算 x(t) 和 y(t) 的频谱函数，并绘制它们的频谱图，观察傅里叶变换的频移特性。",
544 "Response": "定义符号变量 t、w、w0，构造信号 x = heaviside(t+1) – heaviside(t–1)，y = cos(w0*t) * x。使用 fourier(x, t, w) 和 fourier(y, t, w) 分别计算频域表达式 X 和 Y，利用 abs 提取幅度谱，用 fplot 绘制 X 和 Y 的幅度谱，对比观察频移效果。",
545 "Complex_CoT": "信号 x(t) 是宽度为 2 的矩形函数，傅里叶变换结果为 sinc 型函数，其频谱主瓣宽度与时间宽度成反比。y(t) 是调幅信号，调制等效于时域乘积，对应频域卷积，根据傅里叶变换性质，cos(ω₀t) 的频谱包含两个冲击函数，调制会将 x(t) 的频谱移动到 ±ω₀。在 MATLAB 中，首先定义符号变量 t、w、w0，然后构造 x = heaviside(t+1) – heaviside(t–1)，y = cos(w0*t)*x，使用 fourier 计算 X = fourier(x, t, w)，Y = fourier(y, t, w)。用 abs(X)、abs(Y) 提取频谱幅度，使用 fplot 在指定频率范围内绘制图像。比较 X 和 Y 的频谱形状，X 为以零频为中心的低通型，Y 为对称分布在 ±ω₀ 附近的带通型，验证傅里叶变换中的频移定理。"
546 },
547 {
548 "Question": "试用 MATLAB 命令求信号 x₁(t) = sin(2πt)/(2πt) 的傅里叶变换，并绘制其幅度谱和相位谱。",
549 "Response": "定义符号变量 t 和 w，构造表达式 x = sin(2*pi*t)/(2*pi*t)。使用 fourier(x, t, w) 计算傅里叶变换，结果记为 X。使用 abs(X) 和 angle(X) 分别提取幅度谱与相位谱，使用 fplot 绘图展示。",
550 "Complex_CoT": "该信号为标准 sinc 函数形式，其傅里叶变换为矩形函数。在 MATLAB 中使用 syms 定义符号变量 t 和 w，表达式写为 x = sin(2*pi*t)/(2*pi*t)。使用 fourier(x, t, w) 进行符号傅里叶变换，返回频域表达式 X(w)。再分别用 abs(X) 提取幅度谱，angle(X) 提取相位谱。使用 fplot 绘制两者的频率响应图像，观察频域中的带限特性和恒定相位特性。"
551 },
552 {
553 "Question": "试用 MATLAB 命令求信号 x₁(t) = sin[2π(t–2)] / [2π(t–2)] 的傅里叶变换，并绘制其幅度谱和相位谱。",
554 "Response": "定义符号变量 t 和 w，构造表达式 x = sin(2*pi*(t-2))/(2*pi*(t-2))。使用 fourier(x, t, w) 计算傅里叶变换，结果为 X。用 abs(X) 和 angle(X) 提取幅度谱与相位谱，并使用 fplot 绘制。",
555 "Complex_CoT": "该信号为第一个 sinc 函数延迟 2 个单位后的形式，时域平移对应频域乘以 e^{-j2w}。在 MATLAB 中定义符号变量 t 和 w，设 x = sin(2*pi*(t-2))/(2*pi*(t-2))，调用 fourier(x, t, w) 得到频域表达式 X(w)。傅里叶变换结果幅度与第一个信号相同，唯一区别在于相位出现线性变化项 -2w。用 abs 提取幅度谱，angle 提取相位谱，并用 fplot 绘制，以验证时移引起的相位线性变化。"
556 },
557 {
558 "Question": "试用 MATLAB 命令求信号 x₁(t) = [sin(πt)/(πt)]² 的傅里叶变换，并绘制其幅度谱和相位谱。",
559 "Response": "定义符号变量 t 和 w，构造表达式 x = (sin(pi*t)/(pi*t))^2。使用 fourier(x, t, w) 计算傅里叶变换，结果为 X。用 abs(X) 和 angle(X) 提取幅度谱与相位谱，用 fplot 绘图展示。",
560 "Complex_CoT": "该信号为 sinc 函数平方形式，等效于矩形函数自卷积，其傅里叶变换为带宽加倍的三角函数形式。在 MATLAB 中定义符号变量 t 和 w，构造表达式 x = (sin(pi*t)/(pi*t))^2，使用 fourier(x, t, w) 计算频域表达式 X(w)。使用 abs(X) 和 angle(X) 分别获取幅度谱与相位谱，使用 fplot 对其在指定频率范围内绘图。通过频谱观察可确认其主瓣更宽，频谱中心对称，且相位为偶函数。"
561 },
562 {
563 "Question": "试用 MATLAB 求 X(jω) = e^{-9ω²} 的傅里叶逆变换，并绘制其时域信号波形图。",
564 "Response": "定义符号变量 w 和 t，构造频域表达式 X = exp(-9*w^2)。使用 ifourier(X, w, t) 计算傅里叶逆变换，得到时域函数 x(t)。使用 simplify 对结果化简，用 fplot 绘制时域信号波形。",
565 "Complex_CoT": "X(jω) 为高斯型频谱，根据傅里叶变换性质，频域为高斯函数时，对应时域也是高斯函数。在 MATLAB 中，使用 syms 定义符号变量 w 和 t，表达频域函数为 X = exp(-9*w^2)。调用 ifourier(X, w, t) 计算解析形式的逆变换，得到时域函数 x(t)。用 simplify 化简表达式结构，用 fplot 绘制 x(t) 在合适区间内的图像，观察其中心对称的高斯脉冲形状。"
566 },
567 {
568 "Question": "试用 MATLAB 求 X(jω) = cos(πω)/(1 - 4ω²) 的傅里叶逆变换，并绘制其时域信号波形图。",
569 "Response": "定义符号变量 w 和 t，构造频域表达式 X = cos(pi*w)/(1 - 4*w^2)。使用 ifourier(X, w, t) 计算傅里叶逆变换，结果记为 x(t)。使用 simplify 简化表达式，并用 fplot 绘制波形图。",
570 "Complex_CoT": "频谱 X(jω) 为分式与余弦函数乘积，具有偶函数特性。使用 syms 定义符号变量 w 和 t，将频域表达式表示为 X = cos(pi*w)/(1 - 4*w^2)。使用 ifourier(X, w, t) 求解逆变换，得到时域函数 x(t)。为清晰表达，用 simplify 对结果进行代数化简，再使用 fplot 绘制 x(t) 在一定范围内的图像，分析其振荡结构与衰减行为，理解频谱特征对时域波形的影响。"
571 },
572 {
573 "Question": "试用 MATLAB 求 X(jω) = Sa²(ω/4) 的傅里叶逆变换，并绘制其时域信号波形图。",
574 "Response": "定义符号变量 w 和 t，将 Sa²(ω/4) 表达为 sinc²(w/4/pi)。构造表达式 X = sinc(w/(4*pi))^2。使用 ifourier(X, w, t) 得到时域函数 x(t)，用 simplify 化简结果，并用 fplot 绘图。",
575 "Complex_CoT": "Sa²(ω/4) 是带限频谱的平方形式，等价于 sinc² 频谱，对应时域信号为有限支撑的卷积型函数。在 MATLAB 中，用 syms 定义 w 和 t，将频域表达式重写为 X = sinc(w/(4*pi))^2，便于 MATLAB 识别 sinc 函数。使用 ifourier(X, w, t) 得到时域函数 x(t)，用 simplify 进行代数简化。再使用 fplot 绘制其时域波形，观察主瓣宽度与sinc平方结构带来的平滑边缘过渡。"
576 },
577 {
578 "Question": "已知两个门信号的卷积为三角波信号，试用 MATLAB 验证傅里叶变换的时域卷积定理。",
579 "Response": "定义符号变量 t 和 w，构造门函数 x(t) = heaviside(t+1) – heaviside(t–1)，令 h(t) 与 x(t) 相同。使用 conv(x, h, 'same') 或符号积分 int(x(τ)*h(t–τ), τ, –inf, inf) 得到时域卷积。再分别对 x(t) 和 h(t) 使用 fourier 求 X(w)、H(w)，将两者相乘后用 ifourier 得到频域乘积对应的时域信号，最后比较两个结果，使用 simplify 和 fplot 观察其是否相同。",
580 "Complex_CoT": "矩形信号 x(t) 和 h(t) 可通过 heaviside(t+1) – heaviside(t–1) 构造，定义符号变量 t、w、τ。时域卷积通过 int(x(τ)*h(t–τ), τ, –inf, inf) 计算，理论上应得到一个对称的三角波函数。使用 fourier 分别对 x(t)、h(t) 求傅里叶变换 X(w)、H(w)，再计算 Y(w) = X(w)*H(w)，表示频域相乘。调用 ifourier(Y, w, t) 得到频域乘积对应的时域函数，设为 y₂(t)。用 simplify 分别化简卷积结果 y₁(t) 与 y₂(t)，使用 fplot 绘图比较二者波形是否一致，从而验证傅里叶变换中的卷积定理。"
581 },
582 {
583 "Question": "已知系统的微分方程为 y‴(t) + 3y″(t) + 6y′(t) + 3y(t) = 4x′(t) + 2x(t)，求该系统的频率响应，并用 MATLAB 绘制其幅频特性图和相频特性图。",
584 "Response": "定义符号变量 w，应用傅里叶变换性质，将 y‴、y″、y′、x′ 依次表示为 (j*w)^n·Y(w) 或 X(w) 的形式。代入原方程后整理得到 H(w) = Y(w)/X(w)。用 abs(H) 得到幅频特性，用 angle(H) 得到相频特性。使用 matlabFunction 将 H 转为可绘图的函数句柄，定义频率范围后分别用 fplot 绘制幅度和相位曲线。",
585 "Complex_CoT": "该系统为线性常系数微分系统，频率响应可通过傅里叶变换导出。设 Y(w) 为 y(t) 的频域表示，X(w) 为 x(t) 的频域表示，利用傅里叶变换性质，y‴(t) → (j*w)^3·Y，y″(t) → (j*w)^2·Y，y′(t) → j*w·Y，x′(t) → j*w·X，x(t) → X。将这些表达式代入原方程得到：[(j*w)^3 + 3(j*w)^2 + 6(j*w) + 3]·Y = [4(j*w) + 2]·X，进而得到 H(w) = Y/X = [4j*w + 2] / [(j*w)^3 + 3(j*w)^2 + 6(j*w) + 3]。化简后分别计算 abs(H) 和 angle(H)。使用 matlabFunction(H) 生成数值函数句柄，定义频率变量 w = linspace(...)，使用 fplot 绘制幅度和相位曲线，从图中可分析该系统的频率响应特性。"
586 },
587 {
588 "Question": "已知某低通滤波器的频率响应为 H(ω) = 1 / (–ω² + 5jω + 3)，若外加激励信号为 x(t) = 4cos(t) + 2cos(20t)，试用 MATLAB 求其稳态响应。",
589 "Response": "定义符号变量 ω，构造频率响应 H = 1 / (–ω^2 + 5j*ω + 3)。分别令 ω = 1 和 ω = 20，使用 subs 代入得到 H1 和 H2，分别用 abs 和 angle 得到两个频率分量的幅值与相位响应。构造稳态输出为 4*abs(H1)*cos(t + angle(H1)) + 2*abs(H2)*cos(20*t + angle(H2))，使用 fplot 绘制该表达式观察系统响应。",
590 "Complex_CoT": "系统为线性时不变系统，激励为两个不同频率的余弦信号，其稳态响应可用频率响应法分项求解并叠加。输入信号频率分别为 ω1 = 1 和 ω2 = 20，分别代入系统的频率响应函数 H(ω) = 1 / (–ω² + 5jω + 3)，计算 H1 = H(1) 和 H2 = H(20)。再计算 abs(H1)、angle(H1)、abs(H2)、angle(H2)，得到对应频率分量的幅度缩放与相位偏移。根据稳态响应公式，输出为 4·|H1|·cos(t + ∠H1) + 2·|H2|·cos(20t + ∠H2)。最后在 MATLAB 中使用 fplot 绘制该表达式，观察滤波器对两个频率成分的衰减与相位变化特性。"
591 },
592 {
593 "Question": "已知系统的微分方程为 y‴(t) + y″(t) + 5y′(t) + 3y(t) = –5x′(t) + 3x(t)，当外加激励信号为 x(t) = cos(2t) 时，试用 MATLAB 命令求系统的稳态响应。",
594 "Response": "定义符号变量 w，利用导数的傅里叶变换性质将原微分方程变为频域代数方程，求出频率响应 H(w) = Y(w)/X(w)。将 w = 2 代入 H(w) 得到 H(2)，用 abs 和 angle 分别求出幅度与相位。根据频率响应法，稳态响应为 abs(H(2))·cos(2t + angle(H(2)))，用 fplot 绘制该响应。",
595 "Complex_CoT": "将系统微分方程变换到频域，利用傅里叶变换中导数对应 jω 的乘积关系，得到 (jω)^3·Y + (jω)^2·Y + 5jω·Y + 3Y = –5jω·X + 3X。整理后得频率响应 H(ω) = Y/X = (–5jω + 3)/[(jω)^3 + (jω)^2 + 5jω + 3]。在 MATLAB 中，定义符号变量 w，构造 H(w)，再代入 ω = 2，使用 subs 计算 H(2)，通过 abs 和 angle 得到幅值与相位。输入信号为 cos(2t)，稳态响应为 |H(2)|·cos(2t + ∠H(2))。用 fplot 在合适时间范围绘制该表达式，得到系统对余弦激励的稳态输出。"
596 },
597 {
598 "Question": "已知三角脉冲信号为 x(t) = 4 - |t|，(0 ≤ t ≤ 4)，试用 MATLAB 实现该信号经冲脉冲抽样后得到的抽样信号 xₛ(t) 及其频谱。",
599 "Response": "使用 t = 0:dt:4 构造时间向量，定义三角信号 x(t) = 4 - abs(t)。设定抽样周期 Ts，利用 mod(t, Ts) == 0 得到抽样时刻位置，构造单位冲激序列或乘法掩码，将 x(t) 乘以抽样序列得到 xₛ(t)。使用 fft 计算频谱，配合 fftshift 整理频谱对称性，使用 abs 计算幅度谱并用 plot 或 stem 绘图观察频域特性。",
600 "Complex_CoT": "首先设置时间范围 t = 0:dt:4，定义连续时间三角信号为 x(t) = 4 - abs(t)。设定抽样周期 Ts，例如 0.5、1 等，构造抽样脉冲序列 δ_T(t)，可用逻辑条件 mod(t, Ts) == 0 构建抽样位置掩码。将 x(t) 与该掩码逐点相乘得到冲激抽样信号 xₛ(t)。对 xₛ(t) 使用 fft 进行频谱分析，得到离散频域信号，再用 fftshift 对频谱进行中心化处理。通过 abs 取模，获得幅度谱，用 linspace 或自定义方式构造频率轴 f，最后用 plot 或 stem 绘制频谱图，观察主频谱及其镜像频谱重复情况。此过程体现冲激抽样在频域引起的周期扩展特性。"
601 },
602 {
603 "Question": "试用三角脉冲信号来验证抽样定理。",
604 "Response": "构造三角脉冲信号 x(t) = 4 - |t|，定义时间范围 t = 0:dt:T，使用多个不同的抽样周期 Ts 分别对信号进行冲激抽样，得到对应的抽样信号 xₛ(t)。使用 fft 和 fftshift 分析每个抽样信号的频谱，用 abs 提取幅度谱并绘图，对比在不同抽样周期下是否出现频谱混叠。可选用 interp1 对抽样信号进行重构，观察是否能还原原始信号，从而验证抽样定理。",
605 "Complex_CoT": "设置连续时间变量 t = 0:dt:T 构造三角信号 x(t) = 4 - |t|，定义两个抽样周期 Ts₁ 和 Ts₂，其中 Ts₁ 满足抽样定理（Fs₁ > 2B），Ts₂ 不满足（Fs₂ < 2B）。分别使用 mod(t, Ts)==0 构造抽样掩码，将其与 x(t) 相乘得到两个抽样信号 xₛ₁(t) 和 xₛ₂(t)。对三角信号及两个抽样信号使用 fft 和 fftshift 得到频谱，观察 xₛ₁(t) 是否保留原频谱，xₛ₂(t) 是否出现频谱重叠混叠。使用 abs 获取幅度谱并使用 plot 或 stem 进行频谱图绘制。如需进一步验证抽样定理，可用 interp1 对 xₛ₁(t) 和 xₛ₂(t) 进行重构，比较重建结果是否与原信号一致，从而确认抽样率是否足够。"
606 },
607 {
608 "Question": "对于三角脉冲信号 x(t) = 4 - |t|，(0 ≤ t ≤ 4)，假设截止频率为 ωₘ = π/2，抽样间隔为 Tₛ = 1，试用 MATLAB 恢复抽样信号，并计算恢复后的信号与原信号的绝对误差。",
609 "Response": "构造高分辨率时间轴 t = 0:dt:4，定义原始信号 x(t) = 4 - abs(t)。用 Ts = 1 在 t 上提取抽样点和对应时间，记为 xₛ(n) 和 nTs。用 sinc 插值公式恢复信号，对每个 t 计算所有抽样点对其的加权和，即 x̂(t) = ∑ xₛ(n)·sinc((t–nTs)/Ts)。重构后用 abs(x̂ – x) 计算每点绝对误差，可进一步用 max 或 mean 统计误差指标，用 plot 绘制原信号、重构信号和误差。",
610 "Complex_CoT": "设定 t = 0:dt:4 构造原始时间轴，定义三角信号 x(t) = 4 – |t|。设 Ts = 1，构造抽样时刻 tn = 0:Ts:4，记录对应抽样值 xₛ(n) = x(tn)。为了重构信号，应用理想带限重构公式：x̂(t) = ∑ xₛ(n)·sinc((t–nTs)/Ts)，可使用 meshgrid 或循环完成对每个时间点的加权累加。构造重构信号 x̂(t) 后，与原信号 x(t) 按点相减取绝对值得到误差向量，用 max 计算最大误差，mean 计算平均误差。绘图展示原信号与重构信号的重合程度，并绘制误差曲线辅助分析恢复精度。"
611 },
612 {
613 "Question": "已知信号 x(t) = Sa(100t)，计算该信号的奈奎斯特频率，并以不同采样频率对该信号进行采样，画出采样前后信号的频谱，观察随着采样频率的变化，信号频谱有何变化。",
614 "Response": "构造高分辨率时间轴 t = –T:dt:T，并用 sinc(100*t) 定义连续信号 x(t)，使用 fft 和 fftshift 计算并绘制其原始频谱。设定多组采样频率 Fs，分别对信号进行采样，采样方式为从 x(t) 中按采样周期提取样点，抽样后用 zero-padding 扩展为原长度或使用 interp1 插值补全，再使用 fft 分析各采样频率下的频谱。比较不同 Fs 条件下频谱形态，观察是否满足奈奎斯特准则，判断频谱是否出现混叠。",
615 "Complex_CoT": "信号 x(t) = Sa(100t) 实质为 sinc(100t)，其频谱为矩形函数，频率带宽为 100Hz，对应奈奎斯特频率为 200Hz。先设置高分辨率时间轴 t = –T:dt:T，并构造 x = sinc(100*t)，使用 fft 计算频谱并用 fftshift 将直流分量移至中心，绘制频谱观察带限特性。然后选取多个采样频率 Fs，包含 Fs < 200、Fs = 200 和 Fs > 200 三种情况。对每个 Fs 构造抽样时间点 t_sample = 0:1/Fs:T，提取样本 x_sample = x(t_sample)，使用 zero-padding 或 interp1 生成与原长度一致的信号 x_pad。分别对各 x_pad 使用 fft 分析采样后的频谱结构。最后比较不同 Fs 条件下频谱图，验证 Fs ≥ 2B 时频谱保持，Fs < 2B 时频谱发生重叠混叠，直观体现抽样频率变化对频谱的影响。"
616 },
617 {
618 "Question": "已知余弦脉冲信号为 x(t) = E·cos(πt/τ), |t| ≤ τ/2，作该余弦脉冲信号经过冲激抽样得到的信号 xₛ(t) 及其频谱，再利用 xₛ(t) 恢复信号 xᵣ(t)，并求其与原信号的绝对误差。",
619 "Response": "设定时间轴 t = –τ/2:dt:τ/2，构造原信号 x(t) = E·cos(πt/τ)。选择抽样周期 Ts，利用 mod(t, Ts)==0 或下标索引提取抽样点，构造稀疏抽样信号 xₛ(t)。分别对 x(t) 和 xₛ(t) 使用 fft 与 fftshift 得到频谱，用 abs 提取幅度谱并绘图比较。对 xₛ(t) 使用 sinc 插值方法进行重建，得到 xᵣ(t)，计算 abs(xᵣ – x) 得到逐点绝对误差，用 max 或 mean 获取误差统计值，评估重建精度。",
620 "Complex_CoT": "先设定时间轴 t = –τ/2:dt:τ/2，构造原始余弦脉冲信号 x(t) = E·cos(πt/τ)，并限制其定义域范围。选择抽样周期 Ts，构造抽样时刻向量 tₛ，通过逻辑索引或 find 提取对应样值，形成抽样序列 xₛ(t)。将抽样值插入原时间轴对应位置，其他位置设为零，得到 xₛ(t)。对原信号和抽样信号分别进行 fft 和 fftshift 得频谱，对比其主瓣与频谱重复情况，分析抽样带来的频域影响。使用 sinc 插值核构造理想重建公式 xᵣ(t) = ∑ xₛ(n)·sinc((t – nTs)/Ts)，可使用 meshgrid 与矩阵乘法实现高效计算。最后与原信号进行点对点比较，用 abs(xᵣ – x) 得到误差序列，评估重建准确性并绘制误差曲线。"
621 },
622 {
623 "Question": "试用 MATLAB 求函数 x(t) = sin(2t)·u(t) 的拉普拉斯变换。",
624 "Response": "定义符号变量 t 和 s，构造表达式 x = sin(2*t)*heaviside(t)。使用 laplace(x, t, s) 求解其拉普拉斯变换，结果为频域函数 X(s)。如需化简形式可用 simplify(X)。",
625 "Complex_CoT": "首先明确函数定义域为 t ≥ 0，因此需引入单位阶跃函数构造因果性信号。在 MATLAB 中使用 syms 定义符号变量 t 和 s，然后构造表达式 x(t) = sin(2*t)*heaviside(t)。使用 laplace(x, t, s) 调用符号变换命令即可得到解析形式的拉普拉斯变换 X(s)。为便于后续运算或表达整洁，可使用 simplify 对结果进行符号化简。最终结果为标准形式 2 / (s^2 + 4)，可作为频域分析基础。"
626 },
627 {
628 "Question": "试用 MATLAB 求函数 x(t) = e^(–3t)·u(t) 的拉普拉斯变换。",
629 "Response": "使用 syms 定义符号变量 t 和 s，构造表达式 x = exp(–3*t)*heaviside(t)，使用 laplace(x, t, s) 得到其拉普拉斯变换 X(s)。如有需要可使用 simplify 进一步化简结果。",
630 "Complex_CoT": "该函数为指数衰减型乘以单位阶跃函数，定义域为 t ≥ 0。使用 syms 定义变量后，在 MATLAB 中写出 x(t) = exp(–3*t)*heaviside(t)，输入 laplace(x, t, s) 可直接得到其拉普拉斯变换。根据定义，拉普拉斯变换结果应为 1 / (s + 3)，MATLAB 返回的表达式与该形式等价。如表达式含指数项或复杂分母，可使用 simplify 规整化形式，以便后续符号运算或逆变换。"
631 },
632 {
633 "Question": "试用 MATLAB 求函数 x(t) = e^(–2t)·sin(t)·u(t) 的拉普拉斯变换。",
634 "Response": "用 syms 定义变量 t 和 s，构造表达式 x = exp(–2*t)*sin(t)*heaviside(t)，使用 laplace(x, t, s) 计算其拉普拉斯变换，输出为 X(s)。可使用 simplify 对结果整理表达形式。",
635 "Complex_CoT": "该函数为阻尼正弦信号乘以单位阶跃函数，属于常见时域构型。使用 syms 定义变量后，在 MATLAB 中定义表达式 x(t) = exp(–2*t)*sin(t)*heaviside(t)，再调用 laplace(x, t, s) 得到其拉普拉斯变换。按拉普拉斯变换线性性及表格公式，其变换结果应为 1 / ((s + 2)^2 + 1)。若 MATLAB 返回的是展开式或非标准形式，使用 simplify 可以进一步化简为常用表达。该类型频域函数常用于系统响应、阻尼特性等分析。"
636 },
637 {
638 "Question": "利用 MATLAB 求 X₁(s) = (2s + 4) / (s² + 4s + 3) 的拉普拉斯逆变换。",
639 "Response": "使用 syms 定义 s 和 t，构造表达式 X1 = (2*s + 4)/(s^2 + 4*s + 3)，用 ilaplace(X1, s, t) 求其拉普拉斯逆变换。必要时使用 simplify 化简，确保结果以标准形式输出。",
640 "Complex_CoT": "该函数是一个严格真分式，分母为二次可因式化项。首先使用 syms 定义符号变量 s 和 t，在 MATLAB 中写出分式表达式 X1。因为分母 s² + 4s + 3 可因式为 (s+1)(s+3)，若 ilaplace 无法直接识别表达式结构，可以先使用 apart(X1) 进行部分分式分解，再传入 ilaplace 求解。最后用 simplify 化简结果，可得到对应的时域响应形式，通常为指数函数的线性组合。"
641 },
642 {
643 "Question": "利用 MATLAB 求 X₂(s) = (s² + 8) / (s² + 5s + 6) 的拉普拉斯逆变换。",
644 "Response": "定义符号变量 s 和 t，构造 X2 = (s^2 + 8)/(s^2 + 5*s + 6)，再用 ilaplace(X2, s, t) 求解逆变换。可选用 simplify 整理结果。",
645 "Complex_CoT": "X₂(s) 是一个假分式，分子与分母同阶，先使用 long division（即 expand 或 polynomial division）将其表示为多项式加真分式的组合。在 MATLAB 中可以直接使用 ilaplace 求解，系统会自动完成这种形式的转换。但为保证表达式处理清晰，可先使用 quorem 或 expand 将其变形为 A + B/(s² + 5s + 6) 的形式，再分别求解常数项和分式部分的逆变换。最后对整体结果用 simplify 化简为清晰的时域表达式。"
646 },
647 {
648 "Question": "利用 MATLAB 求 X₃(s) = 1 / (s(s + 1)²) 的拉普拉斯逆变换。",
649 "Response": "使用 syms 定义 s 和 t，构造表达式 X3 = 1 / (s*(s + 1)^2)，调用 ilaplace(X3, s, t) 求逆变换，必要时使用 simplify 处理结果。",
650 "Complex_CoT": "该函数为典型的三阶部分分式结构，分母含有一阶项和二阶重根。可使用 apart(X3) 对其做部分分式分解，将其转化为标准形式如 A/s + B/(s+1) + C/(s+1)²，然后再逐项使用 ilaplace 求解。也可将整体表达式直接输入 ilaplace，MATLAB 会自动识别并分项处理。最终结果通常为常数项、指数项及 t·e^(–t) 型项的组合，反映拉普拉斯域中高阶极点对时域行为的影响。"
651 },
652 {
653 "Question": "利用 MATLAB 求 X₄(s) = (s² + 2s + 5) / [(s + 3)(s + 5)²] 的拉普拉斯逆变换。",
654 "Response": "定义符号变量 s 和 t，构造表达式 X4 = (s^2 + 2*s + 5)/((s + 3)*(s + 5)^2)，使用 ilaplace(X4, s, t) 求解，结合 simplify 或 expand 理清输出结构。",
655 "Complex_CoT": "此函数的分母包含一个一阶项和一个二阶重根，整体是一个真分式。在 MATLAB 中构造表达式后，可直接使用 ilaplace 求解。若变换结果不清晰，可用 apart(X4) 将其拆解为三项标准部分分式，例如 A/(s+3) + B/(s+5) + C/(s+5)^2，再逐项进行逆变换求解。最终结果将包含多个指数函数与 t·e^(–at) 型项，可体现复频域中不同极点对时域响应的贡献。"
656 },
657 {
658 "Question": "试用 MATLAB 实现部分展开式法求拉普拉斯变换式 X(s) = (3s² + s – 1) / [s(s – 1)(s + 2)] 的反变换。",
659 "Response": "使用 syms 定义符号变量 s 和 t，构造表达式 X = (3*s^2 + s - 1)/(s*(s - 1)*(s + 2))，用 apart(X, s) 展开为标准部分分式形式，再对每一项分别使用 ilaplace 求逆变换，最后将所有项相加，并用 simplify 化简结果。",
660 "Complex_CoT": "首先分析表达式 X(s) 为真分式，分母由三个互异一次因子构成，适合使用部分分式展开法。使用 syms 定义变量后，构造完整表达式 X(s)。调用 apart(X, s) 将其拆分为三项之和形式 A/s + B/(s - 1) + C/(s + 2)，MATLAB 会返回各项对应系数表达式。接着对每一部分项分别使用 ilaplace 进行逆变换，由于各项均为标准一阶形式，结果应为 a·u(t)、b·e^(t)·u(t)、c·e^(–2t)·u(t) 等形式。最后将所有项累加得到完整的时域函数表达式，并使用 simplify 对最终结果做统一整理。该过程清晰体现了由频域函数到时域信号的结构还原与分项解析过程。"
661 },
662 {
663 "Question": "试用 MATLAB 实现 x(t) = (2 + t)·e^(–3t) 的拉普拉斯变换。",
664 "Response": "使用 syms 定义变量 t 和 s，构造表达式 x = (2 + t)*exp(–3*t)，调用 laplace(x, t, s) 得到 X(s)，可使用 simplify 对结果进行化简。",
665 "Complex_CoT": "x(t) 为 (2 + t)·e^(–3t)，包含一个 t·e^(–at) 型项，其拉普拉斯变换需要利用线性叠加和已知变换规则。在 MATLAB 中，先定义符号变量 t 和 s，构造表达式为 (2 + t)*exp(–3*t)。调用 laplace 函数即可自动处理各项变换，并完成频域相加。最终结果为 X(s) 的解析表达式，可以用 simplify 化简，使其呈现标准分式结构，方便观察极点与系统特性。"
666 },
667 {
668 "Question": "试用 MATLAB 实现 x(t) = e^(–3t)·cos(2t)·ε(t) 的拉普拉斯变换。",
669 "Response": "使用 syms 定义变量 t 和 s，构造表达式 x = exp(–3*t)*cos(2*t)*heaviside(t)，使用 laplace(x, t, s) 求解，结果用 simplify 整理输出。",
670 "Complex_CoT": "该信号为指数衰减余弦调制信号 e^(–3t)·cos(2t)·u(t)，其拉普拉斯变换形式标准，可直接通过 MATLAB 中的 laplace 函数求得。在实现上，先定义符号变量 t 和 s，再写出表达式 x = exp(–3*t)*cos(2*t)*heaviside(t)。调用 laplace(x, t, s) 后，MATLAB 将自动匹配形式并给出 X(s)，其结果一般为二次形式分母的真分式，可用 simplify 精化输出，得到以 s 为变量的标准频域表达式。"
671 },
672 {
673 "Question": "求 X(s) = 1 / [(s – 1)(s – 2)(s – 3)] 的拉普拉斯反变换。",
674 "Response": "使用 syms 定义 s 和 t，构造表达式 X = 1 / ((s – 1)*(s – 2)*(s – 3))，使用 ilaplace(X, s, t) 求解。若想查看中间部分分式展开结果，可用 apart(X, s) 后再分别求解每项的逆变换，最后相加。",
675 "Complex_CoT": "该函数为三个互异一阶实极点构成的真分式，可直接用 ilaplace 求逆变换。若需手动推导结构，可先用 apart 展开为 A/(s – 1) + B/(s – 2) + C/(s – 3)，再对每项分别使用 ilaplace 得到对应的 e^(t)、e^(2t)、e^(3t) 形式，最终结果为三个指数项的线性组合。该过程体现了一阶实极点在时域中的基本响应结构。"
676 },
677 {
678 "Question": "求 X(s) = (5s + 3) / [(s – 1)(s² + 2s + 5)] 的拉普拉斯反变换。",
679 "Response": "用 syms 定义 s 和 t，构造表达式 X = (5*s + 3)/((s – 1)*(s^2 + 2*s + 5))，直接调用 ilaplace(X, s, t) 可获得结果。如需解析结构可先用 apart 展开，再对每项单独求逆变换。",
680 "Complex_CoT": "该函数由一个实极点与一个不可约复二次极点组成，可直接使用 ilaplace 求逆变换。若需分项解析，可先用 apart 展开为 A/(s – 1) + (Bs + C)/(s² + 2s + 5)，再分别对一阶项和复极点项进行逆变换，前者为指数函数，后者变换结果为 e^(–t)·[cos(2t) 或 sin(2t)] 的形式。最终结果是实数指数项与振荡项的组合。"
681 },
682 {
683 "Question": "求 X(s) = e^(–s) / [s(s² + 1)] 的拉普拉斯反变换。",
684 "Response": "定义 s 和 t 后构造表达式 X = exp(–s)/(s*(s^2 + 1))，直接使用 ilaplace(X, s, t) 可自动返回包含单位延迟的时域表达式，或先对无 e^(–s) 部分求逆变换再乘以 heaviside(t – 1) 并将 t 替换为 t – 1。",
685 "Complex_CoT": "该函数包含 e^(–s)，表示单位时间延迟。处理方法为先计算 F(s) = 1 / [s(s² + 1)] 的逆变换，得到非延迟响应 f(t)，再根据拉普拉斯时移性质将结果替换为 f(t – 1)·heaviside(t – 1)。在 MATLAB 中直接调用 ilaplace 对整体表达式求解，系统会自动考虑时移并给出包含单位阶跃函数的正确表达式。"
686 },
687 {
688 "Question": "求 X(s) = e^(–s) / [(s + 2)((s + 2)² + 1)] 的拉普拉斯反变换。",
689 "Response": "定义 s 和 t，构造 X = exp(–s)/((s + 2)*((s + 2)^2 + 1))，直接使用 ilaplace(X, s, t) 即可得到延迟型复频点的时域表达式，含 heaviside(t – 1)。",
690 "Complex_CoT": "该函数包含延迟因子 e^(–s) 与一阶平移复二次项。先构造非延迟部分 F(s) = 1 / [(s + 2)((s + 2)^2 + 1)]，使用 ilaplace 得到对应 e^(–2t)·[cos(t) + sin(t)] 型表达式。结合时移性质，最终形式为上述结果延迟一单位时间，即 f(t – 1)·heaviside(t – 1)。在 MATLAB 中直接使用 ilaplace 处理带有 e^(–s) 的表达式即可自动得到正确形式。"
691 },
692 {
693 "Question": "已知系统函数为 H(s) = (s² + 5) / (s² + 2s + 5)，系统初始状态为 y(0⁻) = 0，y'(0⁻) = –2，输入 x(t) = ε(t)，求系统完全响应 y(t)、零输入响应 yzi(t)、零状态响应 yzs(t)。",
694 "Response": "先构造系统微分方程 y''(t) + 2y'(t) + 5y(t) = x''(t) + 5x(t)。零输入响应由齐次方程 y'' + 2y' + 5y = 0 和初始条件 y(0) = 0, y'(0) = –2 通过 dsolve 求得；零状态响应由 Y(s) = H(s)·1/s 得到，使用 ilaplace 求出 yzs(t)；最终完全响应为两者之和。",
695 "Complex_CoT": "根据 H(s) = (s² + 5)/(s² + 2s + 5)，分母为系统特征多项式，得到系统对应微分方程为 y'' + 2y' + 5y = x'' + 5x。输入为阶跃函数 ε(t)，对应拉普拉斯变换为 1/s，因此 X(s) = 1/s，Yzs(s) = H(s)·X(s) = (s² + 5) / (s(s² + 2s + 5))，使用 ilaplace 求得零状态响应 yzs(t)。零输入响应部分令右侧为 0，即 y'' + 2y' + 5y = 0，结合初始条件 y(0) = 0，y'(0) = –2，使用 dsolve 求得 yzi(t)。将 yzi(t) 与 yzs(t) 相加得到系统完全响应 y(t)。"
696 },
697 {
698 "Question": "已知系统函数为 H(s) = (s + 1) / (s² + 2s + 5)，试用 MATLAB 求出系统函数的零、极点位置，并绘制零极点分布图。",
699 "Response": "使用 tf 函数定义系统的分子和分母多项式，构造传递函数对象 H。使用 zero(H) 获取零点，使用 pole(H) 获取极点。最后调用 pzmap(H) 绘制零极点图。",
700 "Complex_CoT": "系统函数为有理分式，其零点来自分子多项式 s + 1，极点来自分母多项式 s² + 2s + 5。首先在 MATLAB 中使用 tf 定义系统的分子系数 [1 1] 和分母系数 [1 2 5]，形成传递函数对象 H。通过 zero(H) 命令提取零点位置，pole(H) 命令提取极点位置，这些命令返回的是零点和极点的复数坐标。接着使用 pzmap(H) 可视化零极点在复平面上的分布，其中 × 表示极点，○ 表示零点。结果可用于进一步分析系统的稳定性和动态响应特性。"
701 },
702 {
703 "Question": "H(s) = 1 / (s + 2)，绘制系统零极点图并求单位冲激响应 h(t)。",
704 "Response": "使用 tf 函数定义系统 H = tf([1], [1 2])，使用 zero(H) 和 pole(H) 分别求出零点和极点，使用 pzmap(H) 绘制零极点图，使用 impulse(H, t) 计算单位冲激响应并绘图。",
705 "Complex_CoT": "H(s) 是单极点系统，无零点，极点为 s = –2。用 tf 定义传递函数后，zero(H) 返回空集，pole(H) 返回 –2。用 pzmap 可在复平面标出极点位置。impulse(H, t) 返回 h(t)，其理论形式为 e^(–2t)·u(t)，为单指数衰减型响应，用 plot(t, h) 展示其波形。"
706 },
707 {
708 "Question": "H(s) = 1 / (s – 3)，绘制系统零极点图并求单位冲激响应 h(t)。",
709 "Response": "定义 H = tf([1], [1 –3])，用 zero(H)、pole(H) 查看零极点，用 pzmap(H) 绘图，用 impulse(H, t) 生成冲激响应曲线。",
710 "Complex_CoT": "系统无零点，极点为 s = 3。极点在右半平面，表示系统不稳定。用 pzmap 显示极点位置后，再用 impulse 得到 h(t) = e^(3t)·u(t)，为指数发散。应注意数值绘图时可能出现快速增长，时间范围应适当控制以避免溢出。"
711 },
712 {
713 "Question": "H(s) = 2 / (s² + 4)，绘制系统零极点图并求单位冲激响应 h(t)。",
714 "Response": "定义 H = tf([2], [1 0 4])，用 pole 和 zero 提取极点零点，用 pzmap 画图，用 impulse(H, t) 得 h(t) 并绘图。",
715 "Complex_CoT": "分子为常数，说明系统无零点。分母为 s² + 4，对应共轭虚数极点 s = ±2j，系统无阻尼、具有纯振荡特性。h(t) = sin(2t)·u(t) 是单位幅度的正弦波响应。用 pzmap 显示虚轴对称极点，impulse 得出的波形为振荡曲线。"
716 },
717 {
718 "Question": "H(s) = 2s / (s² + 1)²，绘制系统零极点图并求单位冲激响应 h(t)。",
719 "Response": "使用 H = tf([2 0], [1 0 1 0]) 定义系统，zero 和 pole 计算极点零点，用 pzmap 查看分布，用 impulse 得 h(t)。",
720 "Complex_CoT": "分子为线性项，说明零点在 s = 0，分母为二次项平方，极点为 s = ±j，且为二重极点。系统是谐振系统，冲激响应 h(t) 包含 t·cos(t)·u(t) 形式的振荡增长特征。用 pzmap 可见虚轴上双极点和零点重合，impulse 显示的波形具有振幅随时间增长的振荡形态。"
721 },
722 {
723 "Question": "H(s) = 1 / s，绘制系统零极点图并求单位冲激响应 h(t)。",
724 "Response": "用 tf([1], [1 0]) 定义系统，pole(H) 得到极点在原点，使用 pzmap 绘图，impulse(H, t) 得单位阶跃函数响应。",
725 "Complex_CoT": "系统为纯积分器，极点在 s = 0，表示边界稳定但不可控。无零点，冲激响应 h(t) = u(t)，为单位阶跃函数。pzmap 显示极点位于原点，impulse 得到单位阶跃波形。注意积分型系统频率响应慢，图形应在较长时间区间绘制以观察稳定性趋势。"
726 },
727 {
728 "Question": "H(s) = 1 / s²，绘制系统零极点图并求单位冲激响应 h(t)。",
729 "Response": "定义 H = tf([1], [1 0 0])，用 pole 查看双重极点位置，用 pzmap 绘图，用 impulse 生成响应波形。",
730 "Complex_CoT": "系统为二阶积分器，极点为 s = 0 且为二重极点，表示系统边界不稳定且存在持续增长趋势。冲激响应为 h(t) = t·u(t)，即单位斜坡函数。用 pzmap 可见双极点重合于原点，impulse 得到线性增长波形，体现积分器的能量累积特性。"
731 },
732 {
733 "Question": "H(s) = (s + 1) / [(s + 1)^2 + 4]，绘制系统零极点图并求单位冲激响应 h(t)。",
734 "Response": "定义传递函数 H = tf([1 1], [1 2 5])，使用 zero(H) 和 pole(H) 提取零极点，用 pzmap(H) 绘图，再用 impulse(H, t) 生成 h(t) 并绘制波形图。",
735 "Complex_CoT": "将分母展开得 s² + 2s + 5，对应复数共轭极点 s = –1 ± 2j，分子为 s + 1，表示零点在 s = –1。传递函数具有阻尼振荡特性。使用 tf 构造系统，zero 得 –1，pole 得复极点对。用 pzmap 查看分布，impulse 得到的响应应为衰减振荡波形，理论形式类似 e^(–t)·cos(2t)。"
736 },
737 {
738 "Question": "H(s) = s / [(s + 1)^2 + 4]，绘制系统零极点图并求单位冲激响应 h(t)。",
739 "Response": "使用 tf([1 0], [1 2 5]) 定义系统，使用 pole、zero 和 pzmap 分析极点零点，用 impulse 求 h(t) 并绘图。",
740 "Complex_CoT": "分母与前题相同，仍为 s² + 2s + 5，极点为 –1 ± 2j，表示阻尼系统。分子为 s，零点在原点。该结构代表一个原点零点作用于阻尼共轭极点对系统。用 tf 构建模型后，pole 显示极点位置，zero 返回零点为 0。响应应为衰减振荡但有初始上升趋势，h(t) 理论上类似于 e^(–t)·sin(2t) 形式。"
741 },
742 {
743 "Question": "已知系统微分方程为 y''(t) + 5y'(t) + 6y(t) = x'(t) – x(t)，y(0⁻) = 3，y'(0⁻) = 1，x(t) = δ(t) – δ(t–1)，求系统的零输入响应 yzi(t)、零状态响应 yzs(t)、全响应 y(t)。",
744 "Response": "使用 dsolve 求解齐次方程 y'' + 5y' + 6y = 0，结合初始条件 y(0)=3, y'(0)=1 得到 yzi(t)。对非齐次项 x(t)=δ(t)–δ(t–1)，在拉普拉斯域构造系统函数 H(s) = (s–1)/(s²+5s+6)，使用 Y(s)=H(s)·X(s)，X(s)=1–e^(–s)，再用 ilaplace 求得 yzs(t)。最终全响应为 y(t) = yzi(t) + yzs(t)。",
745 "Complex_CoT": "将系统方程转换为 y'' + 5y' + 6y = x' – x，表示右侧为非齐次输入项。首先设 x(t)=0，构造齐次方程，用 dsolve 指定 y(0)=3、y'(0)=1，求得 yzi(t)。接着构造拉普拉斯表达式，X(s)=1–e^(–s)，X'(s)=s·(1–e^(–s))，代入得到 Y(s)=(s–1)(1–e^(–s))/(s²+5s+6)。使用 ilaplace 求得 yzs(t)。两部分相加得到完整响应。也可用 tf 定义系统函数 H(s)，用 impulse 或 lsim 分别对 δ(t) 和 δ(t–1) 输入进行系统响应模拟后相减，得到 yzs(t)。"
746 },
747 {
748 "Question": "已知系统函数 H(s) = 1 / [s(s + 1)(s + 2)]，试用 MATLAB 命令画出其零、极点分布图与时域波形，并判断系统的稳定性。",
749 "Response": "使用 tf([1], conv([1 0], conv([1 1], [1 2]))) 构造系统函数，使用 zero(H)、pole(H) 查看零点与极点，用 pzmap(H) 绘图，用 impulse(H, t) 得到时域响应。根据极点实部是否全小于零判断系统稳定性。",
750 "Complex_CoT": "该系统为三阶系统，无零点，极点位于 s = 0、s = –1、s = –2。零极点分布显示系统存在一个原点极点，因此系统为边界稳定，即响应不发散但不收敛于零。定义传递函数后使用 pole 函数提取极点位置，使用 pzmap 查看极点分布，使用 impulse 可得单位冲激响应 h(t)，其响应为非发散但非零稳态。"
751 },
752 {
753 "Question": "已知系统函数 H(s) = (s + 2) / [(s + 2)(s² + 3)]，试用 MATLAB 命令画出其零、极点分布图与时域波形，并判断系统的稳定性。",
754 "Response": "将 H(s) 约简后为 1 / (s² + 3)，使用 tf([1], [1 0 3]) 定义系统，使用 zero(H)、pole(H) 提取零极点，调用 pzmap(H) 绘图，使用 impulse(H, t) 得冲激响应，判断极点是否全部在左半平面。",
755 "Complex_CoT": "约去公共因子后，系统变为 H(s) = 1 / (s² + 3)，无零点，极点为 s = ±j√3，为虚轴上的共轭复数对。系统为理想振荡系统，冲激响应为持续振荡的正弦波。使用 tf 构造系统函数后，zero 返回空集，pole 返回 ±j√3，使用 pzmap 显示极点位于虚轴上，使用 impulse 得到 h(t) = sin(√3·t)/√3 形式，系统为边界稳定，响应不衰减也不发散。"
756 },
757 {
758 "Question": "利用 MATLAB 实现单位样值序列。",
759 "Response": "定义离散时间变量 n，通常取对称区间如 n = -10:10。使用逻辑表达式 (n == 0) 生成单位样值序列 δ[n]，其中仅当 n=0 时为 1，其余为 0。使用 stem(n, δ) 绘图显示离散序列形状。",
760 "Complex_CoT": "单位样值序列 δ[n] 定义为在 n=0 处取值为 1，其余为 0，是离散时间系统分析的基础信号。在 MATLAB 中可先构造时间索引变量 n，如 n = -10:10，然后用逻辑判断生成 δ[n]：令 x = (n == 0)，该表达式返回一个与 n 等长的数组，n=0 对应位置为 1，其余为 0。最后用 stem(n, x) 绘制离散序列图形，展示 δ[n] 的冲激特性。也可使用 impz 函数生成系统冲激响应作为 δ[n] 的替代表示形式。"
761 },
762 {
763 "Question": "利用 MATLAB 实现单位阶跃序列。",
764 "Response": "先定义离散时间变量 n，取对称区间如 n = -10:10。利用逻辑条件 (n >= 0) 生成单位阶跃序列 u[n]，当 n ≥ 0 时为 1，其余为 0。使用 stem(n, u) 可绘制该离散序列的阶跃特性图。",
765 "Complex_CoT": "单位阶跃序列 u[n] 是最常见的离散时间基本信号之一，定义为 n < 0 时取 0，n ≥ 0 时取 1。MATLAB 中应首先创建离散时间索引向量，如 n = -10:10。然后使用逻辑表达式 x = (n >= 0) 构造对应值向量，逻辑运算自动将满足条件位置赋值为 1，其余为 0，即实现了单位阶跃函数的定义。通过 stem(n, x) 绘图可清晰显示阶跃结构，用 xlabel、ylabel 和 grid 提高可读性。该序列常用于离散系统响应分析中模拟系统在单位输入作用下的累积反应。"
766 },
767 {
768 "Question": "利用 MATLAB 实现矩形序列 R₃(n)。",
769 "Response": "定义时间变量 n，例如取 n = -10:10，利用逻辑条件 (n >= -3) & (n <= 3) 构造矩形序列 R₃(n)，表示在 n ∈ [–3, 3] 区间上值为 1，其余为 0。使用 stem(n, R) 可绘制离散序列图。",
770 "Complex_CoT": "矩形序列 R₃(n) 是长度为 7、中心对称的有限序列，定义域为 n ∈ [–3, 3]，在该区间内取值为 1，其余为 0。在 MATLAB 中，应先构造足够覆盖的离散时间索引向量 n，如 n = -10:10，然后使用逻辑条件构造值向量：x = (n >= -3) & (n <= 3)，该表达式将在指定区间生成值为 1 的布尔向量。为便于后续处理可用 double(x) 转换为数值类型。使用 stem(n, x) 绘图可直观展示该有限宽度矩形信号的时域结构。该构造方法通用于任意宽度和中心的矩形序列建模。"
771 },
772 {
773 "Question": "利用 MATLAB 分别绘制实指数序列：x₁(n) = 1·2ⁿ，x₂(n) = 2·1ⁿ，x₃(n) = 4·0.6ⁿ，x₄(n) = (–0.5)ⁿ，x₅(n) = 3·(–1)ⁿ，x₆(n) = 2·(–1.3)ⁿ。",
774 "Response": "首先定义时间变量 n = 0:20，然后分别使用 A·B.^n 的形式构造六个指数序列。利用 stem(n, x) 绘制每个序列的时域波形，必要时用 subplot 分别展示。注意使用 .^ 运算符处理向量幂运算，负底数会造成交替震荡，底数大于 1 导致发散，0 < 底数 < 1 导致衰减。",
775 "Complex_CoT": "设定统一的离散时间变量 n = 0:20，以便对比多个序列。依次构造各指数信号：x₁ = 2ⁿ 表现为指数增长，x₂ = 常数序列为恒定值，x₃ = 0.6ⁿ 表现为指数衰减，x₄ = (–0.5)ⁿ 在衰减的同时奇偶交替符号，x₅ = (–1)ⁿ 为单位幅度交替震荡，x₆ = (–1.3)ⁿ 是发散同时交替震荡的非稳定信号。每个序列可通过 x = A * B.^n 构造，绘图时使用 stem 命令表示离散信号。可以用 subplot(3,2,i) 布局展示六个图形，对比观察不同底数对指数信号增长、衰减、震荡行为的影响。该类题目旨在训练学生理解指数信号的数学性质及 MATLAB 中的向量化计算方法。"
776 },
777 {
778 "Question": "利用 MATLAB 绘制正弦序列 x(n) = 2·sin(nπ/6 + π/3)。",
779 "Response": "定义时间变量 n，例如取 n = 0:30，利用向量表达式 x = 2*sin(pi*n/6 + pi/3) 构造该正弦序列。使用 stem(n, x) 绘制离散图像，确保图像反映幅值、频率、初相位的变化特性。",
780 "Complex_CoT": "观察该正弦序列的表达式可得：幅值为 2，角频率为 π/6，初相位为 π/3，属标准形式的离散正弦信号。在 MATLAB 中，应先设定足够长的时间范围，例如 n = 0:30，用于捕捉多个周期。然后通过向量方式构造序列：x = 2*sin(pi*n/6 + pi/3)，注意必须使用 pi 表示圆周率，保持角度单位一致性。绘图时使用 stem 函数清晰展示离散信号的点状结构，可辅以 xlabel, ylabel, title 添加说明。通过观察图像周期性、起始相位位置及最大值，可以验证正弦信号构造是否正确。"
781 },
782 {
783 "Question": "利用 MATLAB 绘制虚指数序列 x(n) = e^{j·nπ/6} 的实部、虚部、模和相位随时间变化的曲线，并分析其时域特性。",
784 "Response": "定义时间变量 n = 0:30，利用向量表达式构造复指数序列 x = exp(1j * pi * n / 6)。分别使用 real(x)、imag(x)、abs(x)、angle(x) 提取四个分量，用 subplot 和 stem 绘制实部、虚部、模和相位随 n 的变化图。",
785 "Complex_CoT": "复指数信号 x(n) = e^{j·nπ/6} 是一个单位模旋转序列。其模为恒定 1，实部为 cos(nπ/6)，虚部为 sin(nπ/6)，相位为 nπ/6，呈线性增长。MATLAB 中应先设定时间索引 n = 0:30，然后构造 x = exp(1j * pi * n / 6)，再分别提取：实部 real(x)、虚部 imag(x)、模 abs(x)、相位 angle(x)。通过 subplot(2,2,i) 将四个子图统一展示，使用 stem 或 plot 函数绘图。观察得：实部虚部为正余弦波，周期为 12，模恒为 1，相位为一条斜率为 π/6 的直线，体现复信号在单位圆上匀速旋转的时域特性。"
786 },
787 {
788 "Question": "利用 MATLAB 绘制复指数序列 x(n) = (0.6)^n * e^{jπn/8} 的时域波形，并分析其时域特性。",
789 "Response": "首先定义离散时间序列 n，例如使用 n = 0:50。然后构造复指数信号 x = (0.6).^n .* exp(1j*pi*n/8)。接着分别计算并绘制实部 real(x)、虚部 imag(x)、幅度 abs(x)，推荐使用 stem 函数表示离散信号。为便于比较，可使用 subplot 将三种图像并列显示。观察并分析幅度随时间的衰减趋势，以及实部、虚部的周期性正弦变化。",
790 "Complex_CoT": "1. 识别信号形式 x(n) = (0.6)^n * e^{jπn/8}，是一个复指数序列，包含衰减与相位旋转；2. MATLAB 中复指数表达式可直接使用 `exp(1j*pi*n/8)`，注意使用 `1j` 表示虚数单位；3. 使用 `n = 0:50` 生成时间序列，适当范围便于观察波形演化；4. 使用向量表达式构造完整序列 x；5. 利用 `real(x)`、`imag(x)`、`abs(x)` 分别提取实部、虚部与幅度；6. 用 `stem` 绘图函数绘制离散波形，配合 `subplot` 实现多个图并列显示；7. 分析波形特性：幅度是否按 (0.6)^n 衰减，实部和虚部是否为周期正余弦信号，其周期是否为 16，是否存在随时间衰减的调制效应。"
791 },
792 {
793 "Question": "利用 MATLAB 绘制离散序列 x(n) = 3·2^n 的波形图。",
794 "Response": "定义时间序列 n，例如 n = 0:10，然后构造序列 x = 3*(2.^n)。使用 stem(n, x) 绘制离散波形，添加标题与坐标轴标签辅助分析。",
795 "Complex_CoT": "1. 明确信号是一个实数指数序列，指数底数为 2，系数为 3；2. MATLAB 中指数运算需用点乘形式 `.^` 保证向量化；3. 定义离散时间变量 n 为非负整数向量，如 0 到 10；4. 构造表达式 x = 3*(2.^n)，实现指数信号；5. 使用 stem 函数进行绘图，适合表示离散采样点；6. 观察波形随时间快速增长的趋势，理解其发散特性。"
796 },
797 {
798 "Question": "利用 MATLAB 绘制离散序列 x(n) = 1 + cos(nπ/4 - π/3) 的波形图。",
799 "Response": "定义时间序列 n，例如 n = 0:20，然后构造序列 x = 1 + cos(pi*n/4 - pi/3)。使用 stem(n, x) 绘图，添加标题便于辨识。",
800 "Complex_CoT": "1. 此序列为余弦信号叠加常数项，具有周期性；2. MATLAB 中余弦函数 cos 支持向量运算，π 用 pi 表示；3. 使用向量 n 构造角度表达式 pi*n/4 - pi/3；4. 构造表达式 x = 1 + cos(...) 实现整体序列；5. 采用 stem 函数绘图表达离散性；6. 观察波形的周期与相位偏移，理解 cos(·) 的相位影响。"
801 },
802 {
803 "Question": "利用 MATLAB 绘制离散序列 x(n) = cos(nπ/3) + cos(3n) 的波形图。",
804 "Response": "定义时间序列 n，例如 n = 0:20，然后构造两个部分 x1 = cos(pi*n/3)，x2 = cos(3*n)，再相加得 x = x1 + x2。使用 stem(n, x) 绘制总信号波形。",
805 "Complex_CoT": "1. 此序列为两个余弦信号之和，一个为实频率，另一个为数字频率；2. 构造 pi*n/3 和 3*n 的角度表达式；3. 使用 cos 函数分别计算两个分量；4. 将两个向量相加得到最终信号 x；5. 通过 stem 函数绘图观察合成信号的周期与调制效应；6. 分析不同频率叠加后对波形形状的影响。"
806 },
807 {
808 "Question": "利用 MATLAB 绘制离散序列 x(n) = (–1/3)^n·ε(n) 的波形图。",
809 "Response": "定义时间序列 n，例如 n = -5:10，构造单位阶跃函数 u = double(n>=0)，再构造序列 x = ((-1/3).^n).*u，使用 stem(n, x) 绘图。",
810 "Complex_CoT": "1. 识别信号为条件指数函数，仅在 n≥0 时存在；2. MATLAB 中可用逻辑表达式 n>=0 生成布尔数组，再转为 double 得到 ε(n)；3. 使用点乘 `.*` 限制指数函数作用域；4. 构造 x = ((-1/3).^n).*double(n>=0)，完成整体表达式；5. stem 函数绘制波形，观察其交替正负且逐渐衰减的趋势；6. 分析 ε(n) 对信号截断的影响，以及 (–1/3)^n 的振荡性和收敛性。"
811 },
812 {
813 "Question": "利用 MATLAB 绘制复指数序列 x(n) = 2·e^{(2/5 + j·π/3)·n} 的实部、虚部、模和相角随时间变化的波形图。",
814 "Response": "定义时间序列 n，例如 n = 0:30，构造复指数表达式 x = 2*exp((2/5 + 1j*pi/3)*n)。利用 real(x)、imag(x)、abs(x)、angle(x) 分别计算实部、虚部、模和相角。用 subplot 配合 stem(n, ·) 绘制四个波形图，并分别标注标题与坐标轴。",
815 "Complex_CoT": "1. 识别为复指数信号，包含实指数增长项和旋转项；2. 定义 n 为向量以支持时变绘图，例如 n = 0:30；3. 构造指数部分 (2/5 + 1j*pi/3)*n，确保虚部用 1j 表示；4. 使用 exp 执行复指数运算，生成完整序列 x；5. 通过 real(x)、imag(x)、abs(x)、angle(x) 提取四个分量；6. 用 stem 绘制每个分量的时域波形，使用 subplot 分图展示，保持图像整洁；7. 观察模随 n 增大呈指数增长，实部与虚部呈振荡性增长，相角线性增加。"
816 },
817 {
818 "Question": "利用 MATLAB 绘制复数表达式 x(n) = 1.3·e^{j·π/3} + e^{j·π/6} 的实部、虚部、模和相角的波形图。",
819 "Response": "由于表达式中不含变量 n，x 是复常数。可直接计算 x = 1.3*exp(1j*pi/3) + exp(1j*pi/6)，再用 real(x)、imag(x)、abs(x)、angle(x) 得到四个分量，分别用 stem 绘制为常量序列图，时间序列可设为 n = 0:10。",
820 "Complex_CoT": "1. 分析表达式发现无时间变量 n，x 为一复常数；2. 直接用 exp 计算 x = 1.3*exp(1j*pi/3) + exp(1j*pi/6)；3. 用 real(x)、imag(x)、abs(x)、angle(x) 提取四个分量；4. 为了作图，构造恒定向量，如 real(x)*ones(size(n))；5. 定义时间序列 n = 0:10 作为绘图横轴；6. 用 stem 函数绘制四条恒定波形，使用 subplot 排列展示；7. 强调该信号为常量，波形无时间变化，用于展示复数特性。"
821 },
822 {
823 "Question": "已知 x₁(n) = R₄(n)，x₂(n) = 0.8ⁿ[ε(n) – ε(n–4)]，绘制 x₃(n) = x₁(n)·x₂(n) 的波形图。",
824 "Response": "设定时间变量 n = 0:7。构造 x₁(n) 为宽度为4的单位矩形序列，可使用逻辑条件 (n >= 0) & (n < 4)，并用 double() 转换为数值序列。构造 x₂(n) 为窗限指数序列，使用表达式 0.8.^n .* ((n >= 0) & (n < 4))，实现窗函数乘法。将 x₁(n) 和 x₂(n) 点乘，使用 `.*` 运算得到 x₃(n)。调用 stem(n, x₃) 绘制离散波形图，添加标题和坐标轴标签。",
825 "Complex_CoT": "1. 明确 x₁(n) 是有限支持的单位矩形序列 R₄(n)，在 n ∈ [0,3] 上为 1，其余为 0，可用逻辑判断实现；2. x₂(n) 是在相同区间内有效的指数信号，指数形式为 0.8 的 n 次方，乘以窗函数实现限长；3. 使用向量表达式定义 n = 0:7，使得信号序列有完整显示空间；4. 构造 x₁(n) 时用逻辑判断生成布尔序列，再转为数值类型；5. 构造 x₂(n) 时先对 n 执行元素指数运算，再与窗函数点乘；6. 对 x₁ 和 x₂ 执行元素乘积操作，生成 x₃(n)；7. 用 stem 函数将 x₃(n) 绘图显示，观察其波形随时间变化。"
826 },
827 {
828 "Question": "已知 x₁(n) = R₄(n)，x₂(n) = 0.8ⁿ[ε(n) – ε(n–4)]，绘制 x₄(n) = [x₁(n) + x₂(n)]·[x₁(n) – x₂(n)] 的波形图。",
829 "Response": "定义时间变量 n = 0:7。构造 x₁(n) 为单位矩形序列，方法为 (n >= 0) & (n < 4)，并用 double() 转换为数值向量。构造 x₂(n) 为窗限指数序列，使用 0.8.^n .* ((n >= 0) & (n < 4)) 表达式。先对 x₁(n) 与 x₂(n) 做加法与减法，得到 s1 = x₁ + x₂ 和 s2 = x₁ - x₂，再对 s1 和 s2 做逐点乘法，得到 x₄(n) = s1 .* s2。使用 stem(n, x₄) 绘图并标注图像。",
830 "Complex_CoT": "1. x₁(n) 是一个持续4个采样点的矩形窗函数，在 n ∈ [0,3] 上为 1，其余为 0，通过逻辑条件生成；2. x₂(n) 是窗限的右侧指数序列，通过对指数信号 0.8^n 与时间窗函数相乘构造；3. 选取统一时间范围 n = 0:7，确保后续加减和点乘操作维度一致；4. 使用向量操作分别计算 x₁ + x₂ 和 x₁ - x₂，生成两个中间序列；5. 对上述两个中间结果执行元素乘法，构造目标信号 x₄(n)；6. 用 stem 绘制结果波形，便于观察结构变化；7. 本题要求掌握多个向量信号运算的组织顺序与构造方式，理解组合表达式在时域上的等效展开。"
831 },
832 {
833 "Question": "已知 x(n) = R₄(n)，绘制 -x(–n) 的离散波形图。",
834 "Response": "定义时间序列 n = -5:5，用逻辑表达式 (n >= 0) & (n < 4) 构造原始矩形序列 x(n)，并用 double() 转换为数值型。将时间反转通过对原始序列按中心对称翻转实现，得到 x(–n)。再对结果整体取负号，得到 –x(–n)。使用 stem(n, y) 绘制结果图形，添加标题说明变换类型。",
835 "Complex_CoT": "1. 原始信号 x(n) 为 R₄(n)，其非零区间为 n ∈ [0,3]，构造方式为 (n >= 0) & (n < 4)，再转为 double 类型；2. –x(–n) 包含两个变换：时间反转与幅度反转；3. 时间反转要求将 x(n) 变换为关于原点对称的序列，可通过构造 n 与 -n 的映射或对原信号翻转实现；4. 翻转后再乘以 –1 得到幅度反转结果；5. 使用统一的 n = -5:5 显示反转后的完整区间；6. 使用 stem(n, –x(–n)) 绘图，并设置适当标签展示反转效果。"
836 },
837 {
838 "Question": "已知 x(n) = R₄(n)，绘制 x(n + 3) 的离散波形图。",
839 "Response": "设置时间轴 n = -5:5，构造原始信号 x(n) = double((n >= 0) & (n < 4))。要实现 x(n + 3)，相当于将原信号左移 3 个单位。通过将原时间变量 n 替换为 n_shifted = n + 3 来构造新索引，并重新计算信号值。最后用 stem(n, y) 绘制平移后的波形。",
840 "Complex_CoT": "1. 原始矩形序列 R₄(n) 在 n ∈ [0,3] 处为 1，其它为 0，可用逻辑表达式构造；2. x(n + 3) 是一个左移操作，即信号向时间轴负方向移动，内容不变但位置提前；3. 构造索引变换 n_shifted = n + 3，然后在该新时间点上计算原始 x(n)；4. 或直接构造表达式 double((n + 3 >= 0) & (n + 3 < 4)) 实现平移后序列；5. 统一使用 n = -5:5 范围绘图，观察平移效果；6. 利用 stem 函数绘制，并设置标题说明左移性质。"
841 },
842 {
843 "Question": "已知 x(n) = R₄(n)，绘制 x(–n + 2) 的离散波形图。",
844 "Response": "定义时间变量 n = -5:5，构造表达式 x(–n + 2) 等价于 x(–(n – 2))，即先右移 2，再进行时间反转。可通过构造 n1 = –n + 2，将此映射作为索引代入 x 的构造条件中，用 double((n1 >= 0) & (n1 < 4)) 得到变换后的信号值。最终使用 stem(n, y) 绘图。",
845 "Complex_CoT": "1. x(n) 是单位矩形窗函数 R₄(n)，定义在 n ∈ [0,3]，其构造方法为逻辑判断；2. 表达式 x(–n + 2) 属于复合变换，包括时间反转与时移；3. 可先解析为 x(–(n – 2))，即先平移再反转；4. 将 n 映射到 n1 = –n + 2，对 n1 使用条件判断构造窗函数序列；5. 统一使用 n = -5:5 显示完整变化区间，确保变换前后长度一致；6. 用 stem 绘制变换后波形图，直观展示复合变换效果；7. 本题强调复合变换中操作顺序的重要性，需精准控制映射关系。"
846 },
847 {
848 "Question": "差分方程为 5y(n) + 3y(n–1) + 2y(n–2) = x(n) – 2x(n–1)，其中 x(n) = 2ⁿ·ε(n)，试用 MATLAB 绘制系统的零状态响应。",
849 "Response": "在 MATLAB 中，首先要定义时间范围，例如设置 n = 0:20，表示我们要计算和绘图的时间区间。然后构造输入信号 x(n)，因为 x(n) = 2^n·ε(n)，所以可以用向量表达式实现：x = 2.^n .* (n >= 0)，其中 (n >= 0) 是阶跃函数 ε(n) 的作用。接着初始化输出 y(n)，由于系统是零状态响应，初始所有输出都设为 0，可以用 y = zeros(size(n)) 创建相同大小的零向量。因为差分方程中 y(n) 与 y(n–1)、y(n–2) 有递归关系，所以不能用向量一次性计算，必须用 for 循环来逐步求解。差分方程可以改写成 y(n) = (1/5)·[x(n) – 2·x(n–1) – 3·y(n–1) – 2·y(n–2)]。在 MATLAB 中，从 k = 1（对应 n=1）开始循环，用索引 k–1、k–2 分别表示前一项和前两项。注意：当 k=1 或 k=2 时，x(k–1)、y(k–1) 等前向项可能越界，因此需提前设置这些初始条件为 0。用 for 循环迭代填充 y 向量之后，使用 stem(n, y) 绘制离散时间响应图，并添加合适的标题、坐标轴标签，完成整个过程。",
850 "Complex_CoT": "1. 识别系统为一阶常系数线性差分方程，含递归项 y(n–1)、y(n–2)，因此不能直接用向量运算求解，需要递推计算；2. 输入信号为 x(n) = 2^n·ε(n)，意味着从 n=0 开始是指数上升序列，负时间处为 0；3. 设置统一时间序列 n = 0:20，确保能观察足够长的系统响应过程；4. 利用 MATLAB 向量表达式生成输入信号 x(n) = 2.^n .* (n >= 0)；5. 使用 zeros(size(n)) 创建零向量 y(n) 并初始化，表示系统零状态响应；6. 将差分方程变形为显式求解格式 y(n) = (1/5)[x(n) – 2x(n–1) – 3y(n–1) – 2y(n–2)]，便于编程实现；7. 使用 for 循环对每个时间点 k 计算 y(k)，其中 k 对应时间 n，使用索引操作访问前两项 y(k–1)、y(k–2)，并在边界处对越界索引赋零值；8. 最后用 stem 函数将 y(n) 离散波形绘制出来，观察系统对输入的零状态响应；9. 此题训练学生理解递归系统的数值求解方法与程序结构。"
851 },
852 {
853 "Question": "差分方程为 5y(n) + 3y(n–1) + 2y(n–2) = x(n) – 2x(n–1)，输入为单位冲激信号 δ(n)，试用 MATLAB 绘制该系统的单位样值响应。",
854 "Response": "首先在 MATLAB 中定义时间变量 n = 0:20，表示计算的时间范围。构造输入信号 x(n)，因为是单位冲激信号，x 只在 n = 0 处为 1，其余为 0，因此可以创建一个与 n 等长的零向量，然后将 x(1) 赋值为 1（对应 n=0 的位置）。输出 y(n) 初始化为零向量，长度与 n 相同。因为差分方程中 y(n) 依赖 y(n–1) 和 y(n–2)，所以需要使用 for 循环从 k = 1 开始进行递推计算。将差分方程整理为显式形式 y(k) = (1/5) * [x(k) – 2*x(k–1) – 3*y(k–1) – 2*y(k–2)]，在循环中逐项代入，注意当 k=1 或 k=2 时涉及的前向项可能越界，需要用条件判断处理边界，如超出索引时赋值为 0。迭代结束后，使用 stem(n, y) 绘制 y(n) 的离散波形图，展示系统的单位样值响应。",
855 "Complex_CoT": "1. 明确单位样值响应是系统对 δ(n) 输入的零状态响应，需使用差分方程数值递推求解；2. 差分方程中包含 y(n)、y(n–1)、y(n–2)，为递归结构，因此不能用简单向量表达式，而需依赖迭代计算；3. 输入为单位冲激信号 δ(n)，在 n = 0 时取 1，其余为 0，可用长度为 N 的零向量并设 x(1) = 1 构造；4. 设置时间范围 n = 0:20，覆盖足够响应过程；5. 初始化输出向量 y 为零，用 for 循环从 k = 1 到 length(n) 递推计算 y(k)；6. 在每次迭代中代入公式 y(k) = (1/5)[x(k) – 2x(k–1) – 3y(k–1) – 2y(k–2)]，其中前一项不足时设为 0；7. 使用 stem 函数绘图，将时间序列 n 与输出 y 对应绘制，得到系统在单位冲激作用下的输出响应；8. 本题强调差分方程边界处理、冲激信号构造与递推方法的结合使用。"
856 },
857 {
858 "Question": "系统单位样值响应为 h(n) = (1/2)^n·[ε(n) – ε(n–10)]，输入为 x(n) = ε(n) – ε(n–5)，试用 MATLAB 求其零状态响应 yₙ(n)。",
859 "Response": "首先在 MATLAB 中构造输入信号 x(n) 和单位样值响应 h(n)。分别定义两个时间轴 n1 = 0:4 和 n2 = 0:9，x(n) 在 n1 上为全 1，其他为 0；h(n) 在 n2 上为 (1/2).^n，有效长度为 10。使用逻辑表达式构造窗口并与指数项相乘。接着调用 conv(x, h) 执行线性卷积，得到零状态响应 y(n)，其总长度为 14，即 n = 0:13。最后使用 stem(n, y) 绘制输出响应，添加标题与坐标轴标签。",
860 "Complex_CoT": "1. 识别题目本质为线性时不变系统的卷积求解问题，零状态响应可通过 y(n) = x(n) * h(n) 得出；2. 输入 x(n) 为宽度为 5 的单位序列，可通过 n = 0:4 构造并赋值为 1，其余为 0；3. 系统单位响应 h(n) 为有限长度指数序列，有效范围为 n ∈ [0,9]，构造方法为 (1/2).^n .* ((n >= 0) & (n < 10))；4. 卷积后序列长度为 5 + 10 – 1 = 14，对应时间轴 n = 0:13；5. 使用 MATLAB 内置函数 conv(x, h) 进行线性卷积，自动实现滑动求和；6. 利用 stem(n, y) 将结果绘制为离散时间波形；7. 本题关键在于掌握单位阶跃差表达窗函数的方法，以及利用 conv 函数进行时域卷积求解零状态响应。"
861 },
862 {
863 "Question": "已知线性时不变系统的输入输出分别为 x(n) = R₄(n)，y(n) = 0.8ⁿ·[ε(n) – ε(n–5)]，试用 MATLAB 求系统的单位冲激响应。",
864 "Response": "在 MATLAB 中，首先定义时间变量 n = 0:10，构造输入信号 x(n) 为长度为 4 的单位矩形序列，可用逻辑表达式 (n >= 0) & (n < 4)，再用 double() 转换为数值向量。输出信号 y(n) 为窗限指数序列，可写作 y = 0.8.^n .* ((n >= 0) & (n < 5))。构造好 x 和 y 后，调用 deconv 函数：h = deconv(y, x)，得到系统单位冲激响应 h(n)。用 stem(0:length(h)-1, h) 绘制响应波形，并通过 title、xlabel、ylabel 添加注释信息。",
865 "Complex_CoT": "1. 本题目要求通过输入输出反推系统单位冲激响应，即通过 y(n) = x(n) * h(n) 反求 h(n)，这属于系统的辨识问题；2. 输入 x(n) 是单位矩形序列 R₄(n)，非零区间在 n ∈ [0,3]，输出 y(n) 是窗限指数序列，在 n ∈ [0,4] 为 0.8ⁿ，其余为 0；3. 利用 MATLAB 构造输入输出序列，x(n) 用逻辑表达式构造为全 1 的长度 4 序列，y(n) 用指数函数乘以窗函数生成；4. 使用 deconv(y, x) 执行反卷积操作，该函数自动返回满足 y = conv(x, h) 的冲激响应 h；5. 反卷积结果 h 的长度为 length(y) – length(x) + 1，应结合时间轴 n = 0:length(h)-1 定义横坐标；6. 利用 stem 函数绘图观察单位冲激响应波形；7. 本题重点在于将卷积逆过程建模为可求解问题，并熟练使用 MATLAB 的 deconv 函数完成系统分析。"
866 },
867 {
868 "Question": "差分方程为 4y(n) + y(n–1) – 3y(n–2) = x(n)，输入为 x(n) = cos(nπ/6)·ε(n)，试用 MATLAB 求出系统在 n = 0~30 范围内的零状态响应。",
869 "Response": "首先在 MATLAB 中定义时间序列 n = 0:30。构造输入信号 x(n) 为 x = cos(pi*n/6) .* (n >= 0)。初始化输出向量 y，与 n 同长度，全部设为 0。将差分方程重写为递推形式：y(n) = (1/4) * [x(n) – y(n–1) + 3*y(n–2)]。使用 for 循环从 k = 1 开始计算 y(k)，每一步都代入前一项和前两项的 y 值以及当前的 x(k)。边界条件 y(–1)、y(–2) 设置为 0。循环完成后用 stem(n, y) 绘图，并添加标题和坐标轴标签。",
870 "Complex_CoT": "1. 差分方程含有 y(n)、y(n–1)、y(n–2) 等项，属于线性常系数递归系统，无法向量化直接求解；2. 输入信号为余弦函数调制的单位阶跃信号，可用向量表达式构造；3. 将方程整理为显式形式，便于在 for 循环中逐步计算；4. 初始化 y 为全零向量，准备递推求解；5. 设定时间索引 k = 1 到 31（对应 n = 0 到 30），在每一步中判断边界处理，确保当 k<2 或 k<3 时使用的 y(k–1)、y(k–2) 为 0；6. 每步迭代使用已知值代入公式计算 y(k)；7. 用 stem 函数绘制离散波形图，并标明横轴为时间，纵轴为响应值；8. 强调零状态响应意味着所有初始条件为零，仅输入引起系统响应。"
871 },
872 {
873 "Question": "差分方程为 y(n) + (1/3)y(n–1) – (1/2)y(n–2) = x(n) + 3x(n–1)，输入为 x(n) = (1/4)^n·ε(n)，试用 MATLAB 求出系统在 n = 0~30 范围内的零状态响应。",
874 "Response": "定义时间变量 n = 0:30，在 MATLAB 中构造输入信号 x = (1/4).^n .* (n >= 0)。初始化输出向量 y 为全零。将差分方程整理为递推形式：y(n) = –(1/3)*y(n–1) + (1/2)*y(n–2) + x(n) + 3*x(n–1)。使用 for 循环从 k = 1 开始，逐项代入计算，每一步使用前两项的 y 和当前、前一项的 x 值。对边界 k=1 和 k=2 时的前向项进行判断并补零。递推完毕后使用 stem(n, y) 绘图展示响应。",
875 "Complex_CoT": "1. 该系统是带输入延迟项和输出延迟项的二阶线性差分系统，求解需依赖递推；2. 输入信号为单位阶跃函数调制的指数信号，可用 (1/4).^n 与逻辑条件构造；3. 为确保系统初始状态为零，输出向量 y 用零初始化；4. 重写差分方程为显式形式 y(n) = f(x, y) 形式，方便循环中代入计算；5. 设定索引循环范围与时间一致，即从 k = 1 到 length(n)，每一步计算 y(k) 时需使用前一项和前两项的值，k=1 或 k=2 时手动处理边界；6. 用向量索引方式处理 x(k–1)，边界补零；7. 完成递推后用 stem 绘图，观察系统响应在当前输入激励下的变化；"
876 },
877 {
878 "Question": "差分方程为 3y(n) + y(n–1) – (1/2)y(n–2) = x(n) + 5x(n–1)，求该系统的单位样值响应，并绘制序列波形。",
879 "Response": "在 MATLAB 中，首先定义时间范围 n = 0:30，构造输入信号 x，其中 x(1) = 1，其他为 0，表示单位冲激。初始化输出向量 y 全为零，长度与 n 一致。将差分方程整理为递推形式：y(n) = (1/3) * [x(n) + 5x(n–1) – y(n–1) + 0.5y(n–2)]，从 n = 0 开始用 for 循环逐步计算每个 y(n)。循环中用索引 k 表示当前时间，注意处理 k–1、k–2 越界时将相应项设为 0。计算完成后使用 stem(n, y) 绘制单位样值响应波形，并用 title、xlabel、ylabel 添加图像注释。",
880 "Complex_CoT": "1. 识别差分方程为二阶线性非齐次方程，含前两项输出与前一项输入；2. 单位样值响应对应输入 x(n) = δ(n)，可用 x = zeros(1, N); x(1) = 1 构造；3. 将方程改写为显式形式，目标是将 y(n) 表达为已知的 x(n)、x(n–1)、y(n–1)、y(n–2) 的函数；4. 初始化 y 为全零向量，为递推做好准备；5. 用 for 循环从 k = 1 开始递推计算 y(k)，处理边界 k=1、k=2 时的越界情况；6. 每步迭代都代入公式 y(k) = (1/3) * [x(k) + 5x(k–1) – y(k–1) + 0.5y(k–2)]；7. 最终用 stem 绘图展示离散响应序列；8. 本题重点在于从差分方程推导递推公式，并实现边界控制与序列可视化。"
881 },
882 {
883 "Question": "差分方程为 10y(n) + 6y(n–1) – 3y(n–2) = 3x(n) + x(n–1)，求该系统的单位样值响应，并绘制序列波形。",
884 "Response": "在 MATLAB 中定义时间变量 n = 0:30，创建输入信号 x，x(1) = 1，其余为 0，表示单位冲激。初始化 y 为全零向量。将差分方程重写为 y(n) = (1/10) * [3x(n) + x(n–1) – 6y(n–1) + 3y(n–2)]，利用 for 循环从 k = 1 到 31 依次计算每个 y(k)，遇到索引越界的项如 x(0)、y(–1)、y(–2) 统一设为 0。递推结束后，使用 stem(n, y) 绘图并添加标签说明。",
885 "Complex_CoT": "1. 系统为二阶递归型差分系统，右边含有当前与滞后输入，左边含有当前与滞后输出；2. 构造单位冲激信号 x，方法为创建全零向量后将 x(1) 设为 1；3. 将方程整理为显式形式便于编程求解；4. 初始化输出向量 y 为全零，准备递推计算；5. 从 k = 1 到 31 使用 for 循环，k 表示当前离散时刻，逐步计算 y(k)；6. 递推表达式为 y(k) = (1/10) * [3x(k) + x(k–1) – 6y(k–1) + 3y(k–2)]，在边界 k=1、k=2 时要设置 x(k–1)、y(k–1)、y(k–2) 为 0；7. 计算完成后用 stem 绘制 y(n)，观察单位样值响应的动态特性；"
886 },
887 {
888 "Question": "单位样值响应为 h(n) = (1/3)^n · [ε(n) – ε(n–8)]，输入为 x(n) = R₅(n)，求系统的零状态响应。",
889 "Response": "在 MATLAB 中，定义时间范围 n = 0:20。构造单位样值响应 h(n) 为 (1/3).^n .* ((n >= 0) & (n < 8))，构造输入 x(n) 为 double((n >= 0) & (n < 5))，即宽度为 5 的单位矩形序列。使用 conv(x, h) 计算系统零状态响应 y(n)。卷积后得到的 y(n) 长度为 length(x) + length(h) – 1，可用 stem(0:length(y)–1, y) 绘图，添加 title、xlabel、ylabel 进行标注。",
890 "Complex_CoT": "1. 识别系统为线性时不变系统，零状态响应可通过线性卷积 y(n) = x(n) * h(n) 求得；2. h(n) 是窗限指数序列，非零区间为 n ∈ [0,7]，可用向量化方式构造；3. 输入 x(n) 是单位矩形序列 R₅(n)，在 n ∈ [0,4] 上为 1，其余为 0；4. 使用 MATLAB 中 conv 函数对 x(n) 与 h(n) 做线性卷积，自动得到完整响应；5. 响应序列的时间轴从 0 开始，总长度为 5+8–1=12；6. 用 stem 绘制响应波形，结合标题等增强图像表达效果；"
891 },
892 {
893 "Question": "单位样值响应为 h(n) = (1/3)^n · [ε(n) – ε(n–8)]，输入为 x(n) = (1/4)^n · [ε(n) – ε(n–3)]，求系统的零状态响应。",
894 "Response": "在 MATLAB 中，构造时间变量 n = 0:20。构造 h(n) 为 (1/3).^n .* ((n >= 0) & (n < 8))，输入信号 x(n) 构造为 (1/4).^n .* ((n >= 0) & (n < 3))，即窗限指数序列。用 conv(x, h) 计算输出 y(n)，其有效长度为 length(x) + length(h) – 1。用 stem(0:length(y)–1, y) 绘图，观察系统输出波形。",
895 "Complex_CoT": "1. 该系统卷积响应为有限长度指数信号与另一个有限长度指数信号的卷积；2. x(n) 是三段窗限的右侧指数序列，在 n ∈ [0,2] 上为 (1/4)^n，其余为 0；3. h(n) 是在 n ∈ [0,7] 上定义的指数序列 (1/3)^n；4. 利用向量点乘和逻辑表达式生成 x 与 h；5. 使用 conv 函数计算系统响应 y(n) = conv(x, h)，MATLAB 自动进行线性卷积；6. 输出长度为 3+8–1=10，对应时间轴从 0 到 9；7. 用 stem 绘制响应，验证指数输入在有限支持系统下的衰减卷积效应。"
896 },
897 {
898 "Question": "单位样值响应为 h(n) = (1/3)^n · [ε(n) – ε(n–8)]，输入为 x(n) = cos(nπ/3) · [ε(n) – ε(n–6)]，求系统的零状态响应。",
899 "Response": "在 MATLAB 中设定时间轴 n = 0:20。构造单位样值响应 h(n) 为 (1/3).^n .* ((n >= 0) & (n < 8))，输入信号 x(n) 为 cos(pi*n/3) .* ((n >= 0) & (n < 6))，表示窗限余弦信号。使用 conv(x, h) 得到零状态响应 y(n)，并用 stem(0:length(y)–1, y) 绘图表示。",
900 "Complex_CoT": "1. 输入为带有有限支持窗函数的余弦信号，即在 n ∈ [0,5] 上为 cos(nπ/3)，其余为 0；2. h(n) 是指数衰减窗限信号，在 n ∈ [0,7] 有效；3. 构造 x(n) 与 h(n) 时使用逻辑表达式乘以余弦或指数序列完成窗限；4. 使用 conv(x, h) 求系统输出，等价于 y(n) = x(n) * h(n)；5. 卷积结果长度为 6+8–1=13，输出时间轴为 n = 0:12；6. 用 stem 绘制离散响应波形，清晰展示余弦输入对系统的调制响应特性；"
901 },
902 {
903 "Question": "已知周期序列 x(n) = 5 (2 ≤ n ≤ 6)，其余为 0，周期 N = 10，求其离散傅里叶级数 Xₙ(k) = DFS[x(n)]。",
904 "Response": "在 MATLAB 中，首先构造一周期内的长度为 10 的序列 x(n)。x(n) 在 n = 2 到 6 的位置取值为 5，其他位置为 0，定义方式为：x = zeros(1, 10); x(3:7) = 5，其中下标偏移因 MATLAB 从 1 开始。然后使用 fft(x, 10) 计算该序列的离散傅里叶级数（等价于长度为 10 的 DFT），结果为包含 10 个复数的频谱向量 X(k)。可使用 abs(X) 和 angle(X) 分别查看幅度谱和相位谱。最后使用 stem(0:9, abs(X)) 或 angle(X) 绘制频谱图像。",
905 "Complex_CoT": "1. 题中要求计算周期为 N=10 的离散傅里叶级数 Xₙ(k)，其本质等价于对一周期长度为 10 的序列执行 DFT；2. 根据 x(n) 的定义，在 n = 2~6 上取值为 5，其他为 0，因此应创建长度为 10 的向量，并对对应位置赋值；3. MATLAB 下标从 1 开始，因此 n = 2 对应索引 3，n = 6 对应索引 7，因此设置 x(3:7) = 5 即可；4. 使用 fft 函数计算 DFT，即 X = fft(x, 10)，返回 10 点频谱；5. X(k) 为复数序列，若需进一步分析，可使用 abs() 提取模值，angle() 提取相位；6. 可调用 stem(0:9, abs(X)) 绘制幅度图，辅助理解频域能量分布；"
906 },
907 {
908 "Question": "已知 x₁(n) = sin(πn/4)，x₂(n) = sin[π(n–1)/4]，试用 MATLAB 验证傅里叶级数的时移特性。",
909 "Response": "在 MATLAB 中，首先设定时间范围 n = 0:15，构造原始序列 x1 = sin(pi*n/4)。构造时移序列 x2，可通过将 x1 向右移一位得到，例如使用 circshift(x1, 1)。接着使用 fft(x1) 和 fft(x2) 分别计算两者的离散傅里叶变换，得到频域表达 X1 和 X2。根据傅里叶级数的时移特性，x(n–1) 的频域响应应为 X1(k)·exp(–j·2π·k/N)，其中 N 为序列长度。用 k = 0:N–1 构造旋转因子，并与 X1 相乘得到理论值 X2_theory。对比实际计算得到的 X2 和理论值 X2_theory 的幅度与相位，可通过 abs() 和 angle() 提取，并使用 stem() 绘图验证，若二者吻合则说明时移特性成立。",
910 "Complex_CoT": "1. 傅里叶级数的时移特性指出，x(n–n₀) 的频域响应等于原频谱乘以旋转因子 exp(–j·2π·k·n₀/N)；2. 为验证该特性，需构造 x1(n) 和其右移一位的 x2(n)；3. 在 MATLAB 中设定周期 N=16，并用 sin(pi*n/4) 构造 x1，x2 可通过 circshift(x1, 1) 得到；4. 使用 fft(x1) 和 fft(x2) 分别求得频域表达 X1 和 X2；5. 构造理论旋转因子为 exp(–1j*2*pi*k/N)，k = 0:N–1，并乘以 X1 得到 X2_theory；6. 计算实际 X2 与理论 X2_theory 的幅度和相位并作图对比，检查 abs(X2) ≈ abs(X2_theory)、angle(X2) ≈ angle(X2_theory) 是否成立；7. 若幅度相同、相位相差旋转量，则时移特性验证通过，说明傅里叶级数对时移具有严格的复指数旋转特性。"
911 },
912 {
913 "Question": "试用 MATLAB 求周期序列 x(n) = sin(πn/2) + cos(πn/3) 的傅里叶级数。",
914 "Response": "首先在 MATLAB 中确定周期。x(n) 由两个正弦信号叠加组成，分别为周期为 4 的 sin(πn/2) 和周期为 6 的 cos(πn/3)。因此总周期为两者最小公倍数，即 12。定义时间向量 n = 0:11，用向量表达式生成周期序列 x = sin(pi*n/2) + cos(pi*n/3)。然后使用 fft(x) 对其进行离散傅里叶变换，结果即为该周期信号在一周期内的傅里叶级数系数，共 12 项复数。可通过 abs(X) 查看频谱幅度，angle(X) 查看相位信息，也可用 stem(0:11, abs(X)) 绘制频谱图表示。",
915 "Complex_CoT": "1. 明确傅里叶级数计算可通过对周期信号一周期序列执行 DFT 实现；2. 输入信号为 x(n) = sin(πn/2) + cos(πn/3)，各分量分别周期为 4 和 6；3. 两分量的最小公倍数为 12，因此取 n = 0:11 构造周期序列；4. 用向量表达式构造信号 x(n) 为 sin(pi*n/2) + cos(pi*n/3)；5. 使用 fft(x) 得到长度为 12 的复数序列，即傅里叶级数系数；6. 可选操作包括 abs(X) 得到幅度谱，angle(X) 得到相位谱，用 stem(0:11, ·) 可视化结果；7. 本题关键在于判断周期、构造一周期数据、使用 FFT 提取频域特征，是理解傅里叶级数本质的典型过程。"
916 },
917 {
918 "Question": "设 x₁(n) = R₃(n)，x₂(n) = ∑_{r=–∞}^{∞} x(n + 7r)，求 X₂(k)，并画图表示 x₂(n)、X₂(k)。",
919 "Response": "在 MATLAB 中，首先定义有限序列 x₁(n) = R₃(n)，即在 n = 0,1,2 取值为 1，其余为 0，可用 x1 = double((n >= 0) & (n < 3)) 表示。由于 x₂(n) 是将 x₁(n) 以周期 7 扩展得到的周期序列，只需取一个周期 n = 0:6，在该区间内构造 x₂(n)，满足 x2(n) = x1(mod(n, 3))。接着使用 fft(x2, 7) 计算其 7 点傅里叶级数系数 X₂(k)。最后分别绘制 x₂(n)、X₂(k) 的幅度谱 abs(X2) 和相位谱 angle(X2)，可使用 stem 函数展示结果。",
920 "Complex_CoT": "1. 理解 x₁(n) 为有限长单位矩形序列 R₃(n)，仅在 n = 0 到 2 上为 1，其他为 0；2. x₂(n) 是 x₁(n) 以周期 7 进行无限周期扩展形成的周期序列，因此只需关注一周期内的值，即 n = 0 到 6；3. 利用 MATLAB 的 mod 函数可实现对任意 n 做周期映射：x2(n) = x1(mod(n, 3))，得到一个周期为 7 的周期信号；4. 该周期序列的傅里叶级数系数可通过 DFT 获得，即 X₂(k) = fft(x2, 7)，长度为 7 的复数频谱向量；5. 提取幅度与相位信息可用 abs(X2)、angle(X2) 表示；6. 使用 stem 绘制 x₂(n)、abs(X2)、angle(X2) 随 k 变化的波形，完整表达时域与频域对应关系；"
921 },
922 {
923 "Question": "已知有限长度序列 x(n) = {3, –1, 2, 4, 3, –2, 0, 1, –4}，设 x(n) 右移三位得到 x₁(n)，试分别利用圆周移位的定义和圆周时移特性求出并绘制 x₁(n) 的 N 点 DFT X₁(k) 和 X₂(k)，并将 X₁(k) 和 X₂(k) 进行比较。",
924 "Response": "在 MATLAB 中，首先定义原始序列 x = [3, -1, 2, 4, 3, -2, 0, 1, -4]，长度为 N = 9。圆周右移三位得到 x₁(n)，可通过 circshift(x, 3) 实现。用 fft(x1) 得到 X₁(k)。接着用 fft(x) 计算原始序列的 DFT 得到 X(k)，利用圆周时移特性构造旋转因子 exp(-1j*2*pi*k*n0/N)，其中 n0 = 3，k = 0:N-1，通过 X₂(k) = X(k) .* exp(...) 得到 X₂(k)。分别使用 stem 绘制 |X₁(k)|、∠X₁(k) 与 |X₂(k)|、∠X₂(k)，进行幅度和相位对比验证两者是否相等。",
925 "Complex_CoT": "1. 理解题意为验证 DFT 的圆周时移特性；2. 原始序列长度 N = 9，右移三位即 n₀ = 3，采用 circshift(x, 3) 生成圆周右移序列 x₁(n)；3. 利用 fft(x1) 求得 X₁(k)，代表直接对右移后序列求得的频谱；4. 利用时移特性公式 X₂(k) = X(k)·exp(–j·2π·k·n₀/N)，首先用 fft(x) 得到 X(k)，再计算旋转因子；5. 构造 k = 0:N–1 的索引向量，构建复指数因子与 X(k) 点乘得到 X₂(k)；6. 分别绘制 X₁(k)、X₂(k) 的幅度谱与相位谱，通过 abs() 和 angle() 提取，使用 stem(0:N-1, ...) 展示；7. 比较两者在幅值与相位上的一致性，从图像中验证 DFT 的时移定理是否成立；8. 本题综合考查 DFT 的计算、移位操作以及理论与数值的对应性验证。"
926 },
927 {
928 "Question": "求序列 x₁(n) = ε(n) – ε(n – n₀)，0 < n₀ < N 的 N 点 DFT。",
929 "Response": "在 MATLAB 中，首先定义时间序列 n = 0:N–1，然后用逻辑表达式构造长度为 n₀ 的单位矩形序列：x1 = double((n >= 0) & (n < n0))。该序列在 0 到 n₀–1 上为 1，其余为 0。使用 fft(x1, N) 计算其 N 点 DFT 得到频谱 X1。最后使用 stem(0:N–1, abs(X1)) 和 stem(0:N–1, angle(X1)) 分别绘制幅度谱与相位谱。",
930 "Complex_CoT": "1. 识别 x₁(n) 为一个长度为 n₀ 的单位矩形窗函数，即有限阶跃函数；2. 在 MATLAB 中通过双逻辑条件生成该序列：x(n) = 1 在 n ∈ [0, n₀–1]，其余为 0；3. DFT 长度指定为 N，因此使用 fft(x, N) 计算频谱；4. 频谱为复数向量，包含幅度和相位信息；5. 使用 abs() 和 angle() 分别获取幅度谱与相位谱，用 stem() 绘图分析频率成分；6. 本题重点是窗函数在频域中对应的谱展宽特性。"
931 },
932 {
933 "Question": "求序列 x₂(n) = 4 + cos²(2πn/N)，n = 0, 1, ..., N–1 的 N 点 DFT。",
934 "Response": "在 MATLAB 中，先设定 n = 0:N–1，利用三角恒等式将 cos² 表达为：cos²(θ) = 0.5 + 0.5·cos(2θ)，从而构造序列：x2 = 4.5 + 0.5*cos(4*pi*n/N)。再使用 fft(x2, N) 计算其 N 点 DFT。最后使用 stem 绘制 abs(X2) 与 angle(X2) 的图像。",
935 "Complex_CoT": "1. 识别 x₂(n) 为一个直流项加余弦平方项；2. 使用恒等变换 cos²(θ) = 0.5 + 0.5·cos(2θ) 将其转为基本余弦函数形式，便于编码实现；3. 构造时域序列后，直接使用 fft(x2, N) 得到频谱；4. 幅度谱应呈现出直流分量和双边频率成分，反映调制频率特征；5. 使用 abs() 和 angle() 获取频域特性，通过 stem(0:N–1, ·) 可视化观察其谱结构；6. 本题体现周期实值信号中多项式分量与余弦调制在频谱中的表现形式。"
936 },
937 {
938 "Question": "求序列 x₃(n) = δ(n – n₀)，0 < n₀ < N 的 N 点 DFT。",
939 "Response": "在 MATLAB 中，先用 x3 = zeros(1, N) 初始化为全零向量，再将第 n₀+1 个位置赋值为 1，表示单位冲激在 n₀ 处。使用 fft(x3, N) 计算其 N 点 DFT，频谱结果应为 N 点复指数：X3(k) = exp(–j·2π·k·n₀/N)。可用 abs() 和 angle() 提取其幅度与相位，并用 stem 绘制。",
940 "Complex_CoT": "1. 冲激函数 δ(n – n₀) 在 n = n₀ 处为 1，其余为 0；2. DFT 理论结果应为复指数序列 e^(–j·2π·k·n₀/N)，具有等模、线性相位特性；3. 构造方式是用全零向量并将 x(n₀+1) 设为 1（注意 MATLAB 从 1 开始计数）；4. 使用 fft(x3, N) 得到频谱，理论上模值全为 1，相位为线性下降；5. abs() 和 angle() 分析幅度和相位结构，使用 stem 显示离散谱特征；6. 本题考查 DFT 对离散时间冲激信号的解析与实现能力，是频谱分析的基本案例。"
941 },
942 {
943 "Question": "已知一个 16 点的时域非周期 δ(n – n₀) 信号和阶跃信号 u(n – n₀)，n₀ = 4，用 N = 64 点进行 FFT 变换，作其时域信号图及信号频谱图。",
944 "Response": "在 MATLAB 中，设置时间索引 n = 0:15。构造 δ(n – 4) 的方法是在长度为 16 的全零向量中第 5 项赋值为 1；构造 u(n – 4) 的方法是对 n ≥ 4 的位置赋值为 1。将两信号用 zeros(1, 64) 扩展为 64 点，并将前 16 项填入原信号，剩余部分为 0。分别对这两个 64 点序列使用 fft(x, 64) 计算频域表达。计算得到的频谱结果为复数序列，可通过 abs() 提取幅度谱，通过 angle() 提取相位谱。定义频率索引 k = 0:63，分别使用 stem(n, x) 绘制时域图，stem(k, abs(X)) 绘制频域幅度图。",
945 "Complex_CoT": "1. δ(n – 4) 是在第 5 个采样点为 1 的冲激信号，其余为 0；u(n – 4) 是从第 5 个点开始为 1 的阶跃信号；2. 为满足 FFT 要求，需要将原始信号补零扩展到 64 点长度；3. 使用 fft() 计算扩展后序列的频谱，得到长度为 64 的复数结果；4. 使用 abs() 获取频谱幅度，angle() 获取频谱相位；5. 使用 stem 绘制信号图与频谱图，观察冲激与阶跃信号在频域上的响应特征；6. 频率轴为 k = 0:63，可视化时建议分别绘制 δ(n – 4)、u(n – 4) 及其对应频谱图。"
946 },
947 {
948 "Question": "已知 x(n) 为长度 N = 5 的矩形序列，试用 MATLAB 分析当 FFT 分别取 8、32、64 点时 x(n) 的频谱变化。",
949 "Response": "在 MATLAB 中，首先定义矩形序列 x = ones(1,5)，表示长度为 5 的单位矩形序列。分别对该序列进行 8 点、32 点、64 点的 FFT 变换，可使用 fft(x, N_fft) 实现，其中 N_fft 依次取 8、32、64。此操作会自动对原始长度为 5 的信号进行零填充。对每种变换结果，使用 abs() 提取频谱幅度，定义频率坐标 k = 0:N_fft–1 或 fs = (0:N_fft–1)/N_fft，并用 plot 或 stem 绘图观察。通过对比不同 N_fft 下的幅度谱，可看到频谱主瓣形状相同但分布更细密，说明 FFT 点数越大，频率分辨率越高。",
950 "Complex_CoT": "1. 构造矩形序列 x(n) = 1（n = 0 到 4），其他为 0，可用 ones(1,5) 表示；2. 为实现不同频域分辨率分析，对 x(n) 分别进行 8、32、64 点 FFT 变换；3. 使用 fft(x, N_fft) 可在 FFT 长度大于原信号长度时自动零填充，得到更高分辨率的频谱；4. 提取幅度谱使用 abs() 函数，频率坐标使用 k = 0:N_fft–1 或归一化为 fs = (0:N_fft–1)/N_fft；5. 用 stem 或 plot 分别绘制三种变换下的频谱幅度曲线进行比较；6. 频谱主瓣保持一致，旁瓣位置和细节因点数增多而展现得更清晰；7. 通过该对比可以观察频域插值效果，理解 FFT 点数增加对频谱细节分析的重要作用。"
951 },
952 {
953 "Question": "试用 MATLAB 求函数 x(n) = (1/2)^n 的 z 变换。",
954 "Response": "在 MATLAB 中，首先使用 syms 定义符号变量 n 和 z。然后表达该序列为 f = (1/2)^n，使用 ztrans(f, n, z) 进行 Z 变换。可用 pretty() 查看变换结果的格式化输出，必要时可用 simplify() 简化表达式。最终得到的是 f(n) 的 Z 变换表达式，关于 z 的分式函数。",
955 "Complex_CoT": "1. 明确 x(n) = (1/2)^n 是基本指数型序列；2. 使用 syms n z 定义变量；3. 表达该序列为 f = (1/2)^n；4. 使用 ztrans(f, n, z) 求解其 Z 变换；5. 得到关于 z 的有理函数形式结果；6. 可用 pretty() 展示、simplify() 精简表达式；7. 本质为识别基本类型信号并调用 ztrans 完成解析变换。"
956 },
957 {
958 "Question": "试用 MATLAB 求函数 x(n) = sin(nπ) 的 z 变换。",
959 "Response": "首先在 MATLAB 中定义符号变量 n 和 z。由于 sin(nπ) 的值恒为 0，可直接验证该序列为零序列。构造表达式 f = sin(pi*n)，然后用 ztrans(f, n, z) 执行变换，MATLAB 返回结果为 0。可使用 isAlways(f == 0) 进一步验证函数恒为零。",
960 "Complex_CoT": "1. sin(nπ) 在所有整数 n 处为 0，说明该序列为零序列；2. 使用 syms n z 定义符号变量；3. 构造表达式 f = sin(pi*n)；4. 执行 ztrans(f, n, z)，应返回 0；5. 若有疑问可用 isAlways(f == 0) 检查其恒等为零；6. 此类恒等为零的函数变换结果仍为零，无需继续展开计算。"
961 },
962 {
963 "Question": "试用 MATLAB 求函数 x(n) = n(n–1)/3 的 z 变换。",
964 "Response": "先用 syms 定义符号变量 n 和 z，然后构造表达式 f = n*(n–1)/3。调用 ztrans(f, n, z) 计算其 Z 变换，得到关于 z 的分式表达式。可使用 pretty() 查看整洁形式，也可用 collect() 按 z 展开项进一步观察函数结构。",
965 "Complex_CoT": "1. x(n) = n(n–1)/3 为二次多项式型序列；2. 定义符号变量 syms n z；3. 表达为 f = n*(n–1)/3；4. 用 ztrans(f, n, z) 计算其 Z 变换；5. MATLAB 输出关于 z 的闭式表达；6. 可用 pretty() 查看形式，或使用 collect() 查看按幂次展开结果；7. 本质为对多项式类信号执行符号 Z 变换，得到 z 域系统函数表达式。"
966 },
967 {
968 "Question": "试用 MATLAB 求函数 X₁(z) = (z² + 5z + 1)/(z² – 4) 的 z 反变换。",
969 "Response": "在 MATLAB 中，首先使用 syms 定义符号变量 z 和 n。然后将函数表达为 X1 = (z^2 + 5*z + 1)/(z^2 - 4)。使用 iztrans(X1, z, n) 执行反 Z 变换，结果为 x(n) 的解析表达式。可使用 simplify() 对结果进行代数简化，pretty() 可用于格式化输出结果。",
970 "Complex_CoT": "1. 明确题意是已知 Z 域函数，求其对应的时域序列；2. 定义符号变量 syms z n；3. 将分式表达式 X1 = (z^2 + 5*z + 1)/(z^2 - 4) 构造为符号函数；4. 使用 iztrans(X1, z, n) 求得对应的时间序列 x(n)；5. 使用 simplify() 优化形式，pretty() 展示结果；6. 因分母为 z² – 4 = (z – 2)(z + 2)，属于简单有理分式，iztrans 可自动完成部分分式展开和反变换；7. 最终输出为关于 n 的显式表达形式。"
971 },
972 {
973 "Question": "试用 MATLAB 求函数 X₁(z) = (z² + 3z)/[(z – 1)(z – 2)(z – 3)] 的 z 反变换。",
974 "Response": "在 MATLAB 中，使用 syms 定义符号变量 z 和 n。构造函数 X1 = (z^2 + 3*z)/((z - 1)*(z - 2)*(z - 3))。调用 iztrans(X1, z, n) 可得到对应的时间序列 x(n)。可用 simplify() 对结果进行化简，并通过 pretty() 输出简洁表达式。",
975 "Complex_CoT": "1. 目标为对有理分式执行 Z 反变换；2. 使用 syms z n 定义符号变量；3. 构造表达式 X1 为三阶分母、二阶分子的比值；4. iztrans(X1, z, n) 自动完成部分分式展开与 Z 域反变换；5. 可用 simplify() 简化，pretty() 结构化输出；6. 因包含 3 个一阶极点，反变换结果为 3 项幂函数加权组合；7. 结果为关于 n 的显式解，表示原时域序列 x(n)。"
976 },
977 {
978 "Question": "已知 X(z) = z(z³ + 2z² – 4z + 8) / [(z – 2)²(z² + 4)]，试用 MATLAB 实现部分展开式法求 X(z) 的 z 反变换 x(n)。",
979 "Response": "在 MATLAB 中，使用 syms 定义符号变量 z 和 n，并构造有理表达式 X(z)。使用 partfrac(X, z) 对表达式进行部分分式展开，将其拆解为多个可直接反变换的标准分式。然后使用 iztrans(expr, z, n) 对展开结果的每一项分别执行 Z 反变换，得到对应的时域表达式 x(n)。可使用 simplify() 对结果进行代数化简，并通过 pretty() 查看结构清晰的输出结果。",
980 "Complex_CoT": "1. 定义符号变量 syms z n；2. 构造函数 X(z) = z*(z³ + 2z² – 4z + 8)/[(z – 2)²(z² + 4)]；3. 使用 partfrac(X, z) 对 X(z) 进行部分分式展开，分母包含一个重实极点和一个复数共轭对极点，对应展开结果将含有形如 A/(z–2) + B/(z–2)² + (Cz+D)/(z²+4) 的项；4. 对每个分式项单独使用 iztrans() 求反变换；5. 使用 simplify() 统一化简得到最终的时域表达式 x(n)；6. 若需验证结果正确性，可使用 ztrans(x(n), n, z) 检查是否能还原 X(z)。"
981 }
982]