简介：面向西电通信工程学院“信号与系统”课程大作业的资料包，聚焦多普勒效应在通信系统中的应用分析与建模。内容围绕多普勒效应原理、信号处理模拟及实验验证展开，可帮助学习者完成从理论推导、算法设计到结果分析的全过程。包内共有6个文件，以4个MATLAB脚本（.m）为主，按题目模块实现多普勒信号的生成、调制与频谱分析；另含1份PDF版作业要求和1份Word版报告，分别对应题目解析与结果整理。压缩包仅1.03MB，结构紧凑，便于快速下载。目前已有1002人学习该资源。通过这套资料，读者可同时获得可直接运行的仿真代码、报告写作框架和数据分析方法，既能加深对信号与系统中采样、频谱分析、系统响应等知识点的理解，也能为雷达、声纳等实际场景中的多普勒处理提供可复用的解决思路。对西电通院学生而言，是一份贴近课程要求、便于对照参考的大作业完整方案。 既然点进《通院指南针》这个系列，说明你多半正在为信号与系统这门课的大作业掉头发。我当初接到“多普勒效应”这个题目的时候，第一反应是：这不就是中学物理课讲过的火车鸣笛变调吗？于是找了个公式一套，用MATLAB画出一个频谱搬移的图，然后写了一堆物理原理，觉得自己稳了。
结果老师看了一眼就问我：“你这个作业，跟傅里叶变换有什么关系？跟系统分析有什么关系？”我当场愣住。后来我才意识到，多普勒效应这个题目，看起来像是物理题，其实是信号与系统的应用题。你不仅要会算频移，更要用这门课的语言——时变时延、相位变化、频谱搬移、短时傅里叶变换——把现象解释清楚。这篇文章就是把我重做、答辩、包括后来帮学弟学妹改这份作业的经验完整写下来，给正要写大作业的同学一个可以直接照着走的路子。
1. 先把最容易被带偏的一件事想明白：这门课要的不是多普勒公式
1.1 为什么那么多同学把作业写成了物理实验报告
多普勒效应太“有名”了，所以大家上手极快：查一下频移公式，然后造一个带频移的余弦信号，做完FFT，指着图上移动的峰写一句“可以看到频移”。这种作业老师见得多了，基本不会给高分，因为它本质上是一篇物理实验报告，不是信号与系统作业。
信号与系统这门课的核心词汇是：系统、时域、频域、傅里叶变换、卷积、滤波器。多普勒效应恰好可以完美地装进这套框架里：发射信号是激励，目标反射信道是系统，接收信号是响应。你要做的是用傅里叶变换、频率分辨率、瞬时频率这些工具去分析这个“时变系统”，而不是单纯复述物理现象。这个思路一旦扭转过来，题目的难度其实不高，但高度会上去不少。
1.2 用信号与系统的语言重写“多普勒效应”
从信号与系统的角度看，多普勒效应的本质是一个“时变时延系统”。这句话是整份作业的理论基石。目标运动的时候，回波时延是随时间变化的，时延变化导致信号相位持续变化，相位持续变化在频域里就体现为频率搬移。反过来，如果目标静止，回波只是一个固定时延的副本，不会产生新的频率分量。
在报告里，我建议把发射机、目标信道、接收机画成一个完整的系统框图，然后强调一句话：“只有时变系统才能产生新的频率分量。”这是时不变系统做不到的。物理上叫多普勒效应，信号与系统里就是频谱搬移。写清楚这一点，老师才会觉得你用的是本课程的思维，而不是在百度百科上抄了一段物理科普。
2. 从物理场景到信号模型的完整推导
2.1 静止目标：回波就是带时延的副本
先看最简单的场景。发射信号是单频余弦 x(t) = cos(2π f0 t)，目标在距离 R0 处静止不动，那么回波是：
r(t) = x(t - τ0)，其中 τ0 = 2R0/c
这里的 c 是介质中的波速。用傅里叶变换看一眼：回波的幅度谱没有任何变化，频谱上只是多了一个线性相位 e^{-j2πfτ0}。这说明一个非常重要的事：静止目标的全部信息都在相位谱的斜率里，而幅度谱看不出所以然。
很多同学写到这里就停了，觉得太简单。其实这一步是后面一切推导的地基——你要通过这个最简单的例子向老师证明，你理解“相位携带信息”。测距本质上就是从线性相位斜率里估计时延，这是信号与系统课程里“相位谱”概念的实际应用。
2.2 运动目标：时变时延如何换算成频谱搬移
现在让目标以速度 v 朝接收机运动，距离变成 R(t) = R0 - vt。回波时延也跟着变：
τ(t) = 2R(t)/c = τ0 - (2v/c)t
把时变时延代回去，展开相位：
r(t) = cos[2π f0 (t - τ(t))] = cos[2π f0 (t - τ0) + 2π f0 (2v/c)t]
注意第二项，这分明就是一个新的频率分量：f0 + (2v/c)f0。也就是说，回波的瞬时频率从 f0 变成了 f0 + fd，这里的 fd = 2vf0/c 正是教材上的多普勒频移公式。
这里最值钱的地方在于：我没有背公式，只用了“时延是时间的函数”这一个假设。整个推导链条是物理到数学的自然过渡，中间每一步都用到了信号与系统的时空观：时移、相位求导、瞬时频率。把这个推导过程完整写进报告，属于一眼就能看出课程功底的加分操作。
2.3 连续波和线性调频信号：多普勒的呈现完全不同
如果只做单频连续波，很快会觉得题目太单薄。但换个信号波形，多普勒效应的样子完全不同。这一点非常值得在报告里展开，因为老师想看的就是“你能不能举一反三”。
| 发射信号 | 回波特征 | 多普勒的表现 | 能提取的信息 |
|---|---|---|---|
| 单频连续波CW | 仍为单频，相位连续变化 | 频谱峰值整体搬移 fd | 速度，无法测距 |
| 线性调频LFM | 频率偏移+时间延迟 | 去斜处理后中频出现偏移 | 距离与速度，但二者耦合 |
| 相干脉冲串 | 各脉冲间相位变化 | 慢时间维FFT出现峰值 | 速度，结合时延可测距 |
这个表放进报告，顺便分析一句：为什么单频连续波测不出距离？因为它没有时间标记，所有时刻看上去都一样；而LFM信号有频率随时间变化的规律，所以能同时携带距离和速度信息。这一层分析，就不是“查个公式”的水平了。
3. MATLAB仿真怎么搭：一套能跑又能答辩的框架
3.1 仿真参数怎么定才不会被老师问倒
我见过太多人一上来就把真实射频载频写进代码，比如2.4 GHz，然后按奈奎斯特取采样率，结果仿真1毫秒就需要上百万个点，MATLAB直接卡死，或者图出来全是锯齿。
正确的思路是：仿真演示的是规律，不是复现真实硬件。想展示多普勒频移，载频设成10 kHz、介质速度设成声速343 m/s，效果一样直观，但计算量小好几个数量级。关键参数就四组：
- 载频 f0：决定频移量级，选个便于观察的值。
- 采样率 fs：必须远大于 2(f0 + fd)，满足奈奎斯特。
- 观测时长 T：决定频率分辨率，Δf ≈ 1/T。
- 目标速度 v：决定频移大小，用来和理论值对照。
记住：频率分辨率不是代码精度问题，是观测时长的物理极限。想要分辨出100 Hz的频移，T至少得大于10 ms，否则FFT之后两个峰糊在一起，谁都救不了你。
3.2 一套可以直接改的MATLAB骨架
下面这段代码是我的简化框架，注释写得很详细，你可以直接复制改参数：
%% 多普勒效应基带仿真骨架 clear; clc; c = 343; % 介质速度，空气中声速；模拟雷达可改 3e8 f0 = 10e3; % 载频，这里用中频演示 v = 20; % 目标接近速度 m/s fs = 200e3; % 采样率，远大于 2*(f0 + fd) T = 0.2; % 观测时长 200 ms t = 0 : 1/fs : T - 1/fs; N = length(t); R0 = 50; % 初始距离，单位 m tau0 = 2 * R0 / c; % 初始时延 fd = 2 * v * f0 / c; % 理论多普勒频移 tau = tau0 - (2*v/c) * t; % 时变时延 % 发射信号和回波 x = cos(2*pi*f0 * t); r = cos(2*pi*f0 .* (t - tau)); % 频域观察 X = fft(x); R = fft(r); f = (0:N-1) * fs / N; [~, i0] = max(abs(X(1:N/2))); [~, i1] = max(abs(R(1:N/2))); fprintf('理论fd = %.2f Hz, 实测峰值偏移 = %.2f Hz\n', ... fd, f(i1) - f(i0)); % 时频图：观察频率随时间的变化 spectrogram(r, hamming(512), 256, 1024, fs, 'yaxis');跑出来会得到两个核心结果：FFT频谱峰值向右搬移，时频图中频率线整体抬高。这两张图是答辩时的“主菜”，一定要会解释它们是怎么来的。
3.3 得把三张核心图放出来
很多同学只会贴一张时域波形图，这远远不够。建议至少放三张：
- 时域图：把发射和回波叠在一起，截取一小段显示，说明回波是发射信号的延迟副本，作为基础验证。
- 频域图：FFT后两个谱峰的位置对照，横轴标成Hz，旁边标注理论fd值，做到“仿真结果可被理论验证”。
- 时频图：用spectrogram画，频率线随时间平移，直观展示瞬时频率的概念。
时频图尤其重要，因为它背后是短时傅里叶变换STFT——给信号加窗、分段做FFT、把频率随时间的变化画出来。这本身就是信号与系统课程频谱分析思想的延伸。你放这张图，等于告诉老师你理解的不只是“频移”，而是“频率如何随时间变化”。
4. 我踩过的四个坑，每个都能让图没法看
4.1 坑一：把射频载波当真采样，仿真规模爆炸
我第一次做的版本用2.4 GHz载波，采样率取5 GHz，想仿真0.1秒，结果算了算要5亿个采样点，MATLAB直接内存爆掉。后来才想明白：工程上处理射频信号，标准做法是先用正交下变频搬到基带，而不是对着载波硬算。仿真同理，用中频代替射频，或者直接建复基带模型，可以在计算量降低几个数量级的同时得到完全一样的频移规律。这个思路如果写到报告里，还能顺带展示你了解雷达或通信接收机的架构。
4.2 坑二：FFT看不出频移，不是代码错了，是观测时间不够
有一版代码里，理论频移fd算出来是200 Hz，但FFT图上两个谱峰完全重合。我一度以为是信号写错了，查了半天才发现观测时长只设了5 ms，对应的频率分辨率Δf = 1/T = 200 Hz，频移恰好等于分辨率，当然分不出来。正确做法是让 T 至少取 1/fd 的5到10倍。从那以后我养成了一个习惯：先算理论fd，再反推观测时长，最后才写代码。参数之间是联动的，不能随便拍脑袋。
4.3 坑三：频谱泄漏把多普勒峰压没了
直接用整段信号做FFT，如果观测时长内的载频周期数不是整数，频谱峰会“拖尾巴”，矩形窗的第一旁瓣只比主瓣低约13 dB。当发射信号和回波的能量比较接近时，泄漏出来的旁瓣会把多普勒小峰掩盖得干干净净。
解决办法是加窗，我一般用汉明窗或布莱克曼窗。但要注意，加窗不是免费的：压制旁瓣的同时会稍微展宽主瓣，频率分辨率会变差。所以窗函数的选择本质是在“旁瓣压制”和“主瓣宽度”之间做权衡。答辩时老师特别喜欢问这个问题，你能答出来，就是明显的加分项。
4.4 坑四：fd的正负号搞反
关于正负号，一个简单的记忆方法：目标靠近时，回波时延越来越小，信号被“压缩”，周期变短，频率升高，所以fd为正；目标远离则相反，fd为负。写代码的时候，先确定R(t)是R0 - vt还是R0 + vt，再代入相位推导，不要最后才猜符号。检查也很简单：把v改成负值，看FFT峰是左移还是右移，和理论对照即可。
5. 报告和答辩里真正加分的“信号与系统味道”
5.1 从“出图”到“出结论”：加参数扫描
只给一组参数下的截图，老师看完觉得“哦，做了”，但不会留下深刻印象。拿高分的关键是让参数跑起来：把速度从5 m/s扫描到50 m/s，记录每个速度下FFT实测的谱峰偏移，画出“理论fd vs 实测fd”的对照曲线。这张图能说明你的仿真程序不是一次性脚本，而是可以反复做定量验证的工具。再进一步，可以分析误差来源：FFT分辨率不够？窗函数展宽了主瓣？采样率偏低？每一项都能追溯到信号与系统的核心概念上。
5.2 模糊函数和时频分布：课程含量一眼可见
如果想做出点别人做不出来的东西，强烈建议加一个模糊函数图。对有雷达基础的老师来说，看到模糊函数就知道你学过或主动查过更深的内容。单频连续波的模糊函数在时延维是一条水平脊，说明它完全没有距离分辨能力；LFM信号的模糊函数则会倾斜，说明距离和速度存在耦合。配合前面那张“不同波形下的多普勒表现”表，理论与仿真互相印证，整份作业的层次一下就上来了。
另外可以做一个“观测时间对多普勒分辨率影响”的实验：固定其他参数，T分别取1 ms、10 ms、100 ms，画出三张频谱图，说明随着观测时间增加，谱峰越来越窄，多普勒频率估计越来越准。这个实验的本质，就是课程里“时间分辨率与频率分辨率成反比”这条主线。
5.3 报告结构和答辩追问清单
报告别按“摘要、背景、原理、结果、心得”这种万能模板写，建议按“问题—模型—仿真—验证—结论”的逻辑：
- 问题：多普勒效应为什么可以看成时变系统。
- 模型：从时变时延出发，推导fd，用傅里叶变换解释频谱搬移。
- 仿真：基带等效建模，FFT与STFT观察频域与时频特征。
- 验证：理论fd与仿真结果对比、参数扫描、误差分析。
- 结论：用两三句话总结“时域相位变化如何映射到频域搬移”。
答辩前，把下面三个问题自己过一遍：为什么可以用中频代替射频做仿真？fd为什么是2vf0/c而不是vf0/c？观测时长为什么会影响频率分辨率？这三个问题覆盖了建模、推导和仿真验证三个环节，基本就是老师最爱问的方向。
最后再分享一点我个人的体会。做完这个大作业之后，我记住的其实不是那个频移公式，而是“遇到一个现象，先问它的频谱长什么样”这种思维习惯。多普勒效应在时域里只是回波相位在偷偷变化，映射到频域就成了肉眼可见的频移。信号与系统这门课真正留给你的，就是这套在时域和频域之间来回切换的直觉。所以做作业之前别急着打开MATLAB，先拿笔把模型推一遍，想清楚每张图要说明什么，图画出来心里才有底。这套从物理现象到数学模型再到仿真验证的流程，才是你这次作业最该带走的东西。
本文还有配套的精品资源，点击获取
