北京工业大学信号处理工程应用训练

合集下载

“信号处理工程实践”课程与工程素质培养

ｔｎｉｅｒｎａｉｏｏｌｇｔｄｎｓｏＥｇｎｅｉｇＱｕｌｙｆｒＣｌｅＳｕｅｔｔｅ
Ｌｕｗｅ，ＩＲ・ｉＷＵＬ－ｎ，ＵＨｕ－ｎｉａｇＤＯｉｉｇｆｊ
（ｃｏｌｏｌｃｏｉＩｏｍｔｎａｄＣｎｒｌｎｉｅｒｇ，ｅｉｎｅｉＴｃｎｌｙＢｉｎ０１４，ｈｎ）ＳｈｏｆＥｅｒｎｃｎｒａｉｎｏｔｇｎｅｎＢｉｎＵｉｒｔｏｅｈｏｏ，ｅｉｇ１０２ＣｉａｔｆｏｏＥｉｊｇｖｓｙｆｇｊ
四个实验，来巩固理论教学。为了进一步提高学用生应用理论知识、决问题的能力和工程素质，年解去
开始，我院开设了“ 信号处理工程实践” 课程。
音检测；基于Ｌ ② ＭＳ最小均方误差法的语音降噪；
１ “ 号处理工程训练＂信的课题来源
第３４卷
第３期
电气电子教学学报
ＪＵＲＮＡＬＯＦＥＥＯＥ
Ｖｏ．４ＮＯ３１３．
２１０２年６月
Ｊｎ２２ｕ．０１
“ 号处理工程实践＂信课程与工程素质培养
李如玮，毋立芳，窦慧晶
（北京工业大学电子信息与控制工程学院，京１０２）北０程是电子信息课
类大学生必须掌握的基础知识。为了帮助学生掌握
工作的几名教师，自己的研究方向中提炼出几个从
典型的小课题，最终汇成 “ 信号处理工程训练 ” 程课

北京工业大学2018年全日制硕士研究生招生学科目录

原则上只接收建筑学、城乡规划学专业的考生报考。
0833城乡规划学
_01城乡规划与设计
_02住区规划与社区发展
_03城市设计与园林景观设计
_04古城与历史街区保护更新
_05城乡防灾减灾规划
11
①101思想政治理论
②201英语一
③633城市规划原理
④503城市规划学术快速设计
原则上只接收建筑学、城乡规划学(含园林景观)专业的考生报考。
_01可再生能源利用及先进环境能源理论与技术
_02强化传热传质理论与工程应用
_03制冷低温系统及其环保节能理论与技术
_04车辆及动力系统节能、净化与控制
25
①101思想政治理论
②201英语一
③301数学一
④851传热学Ⅱ或852工程热力学
0817化学工程与技术
_01绿色化学与精细有机化工
_02工业催化与纳米科学
③662普通物理Ⅰ
④828激光原理或861量子力学
0803光学工程
_01激光先进制造技术
_02非金属材料激光制造与改性
_03先进光纤激光技术
_04超短脉冲激光技术
_05高功率固体激光技术
_06高功率半导体激光技术
_07微纳光学与微纳制造
_08生物光子学
_09光纤光子学
33
①101思想政治理论
②201英语一
20
①101思想政治理论
②204英语二
③303数学三
④432统计学
009材料科学与工程学院
154
0805材料科学与工程
_01生态环境材料与资源循环技术
_02稀土、难熔金属等功能材料
_03高性能结构材料技术
_04先进材料加工技术

北京工业大学信号与信息处理研究室团队简介

维普资讯
北京工业大学信号与信息处理研究室团队简介
北京工业大学信号与信息处理研究室是北京市自然科学基金委重点建设的基础性研究室北京工业大学２１１工程重点建设的研究室，有在校博士硕士生３余名，现０主要研究方向为：图像／频编视码、输、传处理及相关应用，究领域涉及图像／频信息的采集、研视处理、分析、传输、存储、压缩、建重以及医学图像分析处理，５是中医舌象采集与分析技术等。特３１研究室主任博士生导师沈兰荪教授为ＩＥＥＥ学会高级会员、国际欧亚科学院院士、北京市突出贡献专家。历任中国科技大学电子工程系主任、国科技大学高技术学院副院长、京工业大学副校中北长、北京市自然科学基金会副会长等。现任中国大百科全书（二版）电子学学科副主编、国微弱信号中检测学会副理事长、国家自然科学基金委信息科学部专家评审组成员，以及中国科技大学等校兼职
教授等职。
近年来，该研究室在沈兰荪教授带领下先后承担和完成了国家自然科学基金、家８３目、国６项国家９３目、７项国家攻关及部委等５个科研项目，表论文４）多篇，０发（（）申请专利１余项，０获得北京市科技进步奖、国家发明展览会奖等多项奖励。先后出版了《数据采集技术》《、图像编码与异步传输》《、小波编码与网络视频传输》《医舌象的采集与分析》、中等著作，共培养博士３余人，－５余人。０硕ｋ０－－研究室先后与中国科技大学、悉尼大学、港理工大学、香北京中医医院、中国中医科学院、北京中医药大学等多家研究单位开展合作研究，并与早稻田大学、台湾中山大学、香港浸会大学等建立了合作关

2023大学_最新电子信息工程专业大学排名

2023最新电子信息工程专业大学排名电子信息工程专业大学排名专业概述电子信息工程是一门应用计算机等现代化技术进行电子信息控制和信息处理的学科，主要研究信息的获取与处理，电子设备与信息系统的设计、开发、应用和集成。

现在，电子信息工程已经涵盖了社会的诸多方面。

电子信息工程专业是集现代电子技术、信息技术、通信技术于一体的.专业。

本专业是一个电子和信息工程方面的专业。

本专业学生主要学习信号的获取与处理、电厂设备信息系统等方面的专业知识，受到电子与信息工程实践的基本训练，具备设计、开发、应用和集成电子设备和信息系统的能力。

该专业是前沿学科，现代社会的各个领域及人们日常生活等都与电子信息技术有着紧密的联系。

全国各地从事电子技术产品的生产、销售和应用的企事业单位很多.，随着改革步伐的加快，这样的企事业单位会越来越多。

为促进市场经济的发展，培养一大批具有大专层次学历，能综合运用所学知识和技能，适应现代电子技术发展的要求，从事企事业单位与本专业相关的产品及设备的生产、安装调试、运行维护、销售及售后服务、新产品技术开发等应用型技术人才和管理人才是社会发展和经济建设的客观需要，市场对该类人才的需求越来越大。

为此电子信息工程专业的人才有着广泛的就业前景，毕业生可从事电子设备、信息系统和通信系统的研究、设计、制造、应用和开发工作。

注重培养电子信息技术基础知识与能力;具有电子产品的装配、调试及设计的基本能力，具有一般电子设备的安装、调试、维护与应用能力;具有对办公自动化设备的安装、调试、维修和维护管理能力;具有对通信设备、家用电子产品电路图的阅读分析及安装、调试、维护能力;具有对机电设备进行智能控制的设计和组织能力;具有阅读相关专业英语资料能力;计算机技术应用能力达到计算机等级四级要求水平。

课程设置主干学科：电子科学与技术、信息与通信工程、计算机科学与技术。

主要课程：电路理论系列课程、计算机技术系列课程、信息理论与编码、信号与系统、数字信号处理、电磁场理论、自动控制原理、感测技术等。

通信工程

通信工程通信工程(也作信息工程，电信工程，旧称远距离通信工程、弱电工程)是电子工程的一个重要分支，同时也是其中一个基础学科。

该学科关注的是通信过程中的信息传输和信号处理的原理和应用。

一、专业简介通信专业实验室通信工程（Communication Engineering）专业是信息科学技术发展迅速并极具活力的一个领域，尤其是数字移动通信、光纤通信、Internet网络通信使人们在传递信息和获得信息方面达到了前所未有的便捷程度。

通信工程具有极广阔的发展前景，也是人才严重短缺的专业之一。

本专业学习通信技术、通信系统和通信网等方面的知识，能在通信领域中从事研究、设计、制造、运营及在国民经济各部门和国防工业中从事开发、应用通信技术与设备。

毕业后可从事无线通信、电视、大规模集成电路、智能仪器及应用电子技术领域的研究，设计和通信工程的研究、设计、技术引进和技术开发工作。

近年来的毕业生集中在通信系统、高科技开发公司、科研院所、设计单位、金融系统、民航、铁路及政府和大专院校等。

本专业本着加强基础、拓宽专业、跟踪前沿、注重能力培养的指导思想，培养德、智、体全面发展，具有扎实的理论基础和开拓创新精神，能够在电子信息技术、通信与通信技术、通信与系统和通信网络等领域中，从事研究、设计、运营、开发的高级专门人才。

二、专业综合介绍引子细心的你是否留意到，十年前港台电影中黑帮大佬手里可以用来砸人的“大哥大”，早已变得如此纤细轻巧、色彩缤纷，并且飞入寻常百姓之手；从前只有数月飞鸽传书才能联系的国外亲友可以用简单方便快捷的伊妹儿（E-mail）互致问候、即时聊天，甚至装上摄像头开个网络会议！这一切都应该归功于通信工程（Communication Engineering）技术的迅猛发展。

如果让科学家们选出近十年来发展速度最快的技术，恐怕也是非通信技术莫属。

那么让我们来多了解一下这个年轻而又略显“神奇”的专业吧。

专业简介通信技术是以现代的声、光、电技术为硬件基础，辅以相应软件来达到信息交流目的。

信息类高年级本科生《信号处理新技术概论》教学实践与思考

摘要：《信号处理新技术概论》是信息类高年级本科生的一门重要选修课程。开设该课程的目的是让学生了解本专业领域的热点内容及其发展动态。笔者认为讲好这门课程，应做好以下几个方面的工作：注重演化理论和方法的历史，讲清楚其产生的来龙去脉；注重数学公式的物理概念讲解；理论与方法同实例结合；设计多种类型的实验供学生选做，以此巩固他们
如果硬性地让学生去记一些公式，这样容易增大学生学习
工作还是工程类工作都是大有裨益的。该课程属于进展类课程，如何开好该类课程对我们来说是一个小小的挑战。以电子信息类本科生为例，笔者在２００７年首次在我校开设了《信号处理新技术概论》课程，这些课程内容涵盖了当前信号处理的热点，之所以加上 “ 简介 ” 这个词是因为每一部分所包括的内容都非常丰富，而在讲课中不可能做到面面俱到。在本课程大纲修订时，参与讨论的老师们一致认为讲这门课程的起点高、难度大、内容广，应该重点介绍基本概念、基本原理和基本方法。例如神经网络部分我们介绍典型的ＢＰ神经网络。笔者第一次上课时对学生做了一下调查，主要了解一下他们学习该课程的目的和认识，以及通过学习该课程他们最想得到什么。综合绝大多数同学的答案，希望学到和日常生活紧密相连、有实践性、能在以后的生活中用到的知识。由此看出学生学习该课程的目的性很强，很多同学希望将所学的内容亲自实践一下。笔者认为在讲授课程的过程中，在“ 授之以鱼 ” 的同时，更要注重 “ 授之以渔” ，逐步培养学生的创新思维和创新素质。使学生学会掌握学习的方法并运用于实践中。本着这一基本思路，笔者认为讲好这门课程，要处理好以下几点。

2022年北京工业大学信息管理与信息系统专业《计算机网络基础》科目期末试卷A(有答案)

2022年北京工业大学信息管理与信息系统专业《计算机网络基础》科目期末试卷A（有答案）一、选择题1、当数据由主机A送传至主机B时，不参与数据封装工作的是（）。

A.物理层B.数据链路层C.网络层D.传输层2、若甲向乙发起一个TCP连接，最大段长MSS-1KB，RTT-5ms，乙开辟的接收缓存为64KB，则甲从连接建立成功至发送窗口达到32KB，需经过的时间至少是（）。

A.25msB.30msC.160msD.165ms3、使用两种编码方案对比特流01100111进行编码的结果如图所示，编码1和编码2分别是（）A.NRZ 和曼彻斯特编码B.NRZ 和差分曼彻斯特编码C.NRZ-I和曼彻斯特编码D.NRZ-I和差分曼彻斯特编码4、无法隔离冲突域的网络互连设备是（）A.路由器B.交换机C.集线器D.网桥5、为了纠正2比特的错误，编码的海明距应该为（）。

A.2B.3C.4D.56、X台计算机连接到一台YMbit/s的集线器上，则每台计算机分得的平均带宽为（）。

A.XMbit/sB.YMbit/sC.Y/XMbit/sD.XYMbit/s7、使用集线器连接局域网是有限制的，任何两个数据终端设备之间允许的传输通路中可使用的集线器个数最多是（）。

A.1个B.2个C.4个D.5个8、某局域网采用SNMP进行网络管理，所有被管设备在15min内轮询一次，网络没有明显拥塞，单个轮询时间为0.4s，则该管理站最多可支持（）个设备。

A.18000B.3600C.2250D.900009、下面关于POP3，（）是错误的。

A.由客户端选择接收后是否将邮件保存在服务器上B.登录到服务器后，发送的密码是加密的C.协议是基于ASCII码的，不能发送二进制数据D.一个账号在服务器上只能有一个邮件接收目录10、ICMP报文的传输方式是（）。

A.无连接的UDP数据报形式传送B.面向连接的TCP报文形式传送C.放在IP数据报的首部字段中传送D.放在IP数据报的数据字段中传送二、填空题11、HTML 语言是一种____________语言。

北京工业大学信号与系统考研导师介绍资料

29研究与开发3智能手机应用软件开发智能手机应用软件开发4物联物联传感网络协议研究开发传感网络协议研究开发5lteadvanced移动通信关键技术研究移动通信关键技术研究对考生的要求1具有通信信号网络方面的基础知识具有通信信号网络方面的基础知识2勤奋踏实学习能力强勤奋踏实学习能力强3能严格要求自我能严格要求自我博硕士生导师介绍博硕士生导师介绍导师姓名张新峰副教授硕导性别男联系方式北京工业大学信号与信息处理研究室10012467391587802学历学位研究生博士研究方向介绍信号与信息处理方向l图像处理分析与
在数字图像处理研究方面有专门的“光电成像跟踪实验室”，从事光电目标检测、识别、跟踪，远地点目标跟踪等技术研究。
对考生的要求
• 首先要有严谨做学问的精神，求知欲强，工作勤奋，有合作精神，读研期间不能工作实习。
• 要求有扎实的数字图像处理、计算机网络、数值分析等方面知识，熟练的C语言编程能力，及一定的DSP硬件知识，动手能力强。
主要研究方向包括：1、多媒体信息安全；2、移动流媒体系统；3、图像/视频编码与网络传输；4、无线视频传感器网络；5、遥感图像的压缩与处理；6、图像/视频超分辨率复原；7、三维数字处理技术与应用
1、英语、数学基础好；
2、能熟练进行C语言编程（或者C＋＋、.NET等）；
3、动手能力强；
4、有团队精神，做人踏实、勤奋努力。
志愿的排序，因此填报志愿请慎重！如果你的三个志愿都未被导师选上，我们会通过导师抽签的方法确定你的导师。
关于导师选择学生的说明
1. 流程控制与管理：吴强，按照以下流程的顺序进行并控制时间
2. 流程安排
（1）第一阶段：导师根据学生填报的第一志愿选择自己的学生, 仅第一志愿可选；时间长度20分钟

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

北京工业大学通信系统工程应用训练报告专业：通信工程学生姓名：刘莹莹指导教师：席大林完成时间：2016年4月29日目录训练十一 DFT性质研究 (1)训练十二 DFT及抽样定理研究 (13)训练十三数字滤波器制作 (20)训练十四 IIR数字滤波器设计与实现 (25)训练十五线性卷积计算 (46)训练十六 FIR数字滤波器设计与实现 (55)训练十一 DFT性质研究验证dft函数正确性设置原始输入信号为x[8]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,0},{8,0}}，将输入信号x[8]进行DFT正变换，dft(X,x,8,1)，输出保存在X[8]，如下：可以看到，输入信号x(n)已经变换到频域X(k)，且仍为8位。

再对X[8]进行DFT反变换，dft(x,X,8,-1)，重新得到x[8]，观察得到的输出与原始输入数据是否相同。

结果如下：可以看到，输出的x[8]取值仍为x[8]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,0},{8,0}}，证明经过DFT正反变换后，信号能够恢复原始信号。

根据帕塞瓦尔定理，应有时域、频域总能量相等：。

经过计算，时域、频域能量和分别为，证明时域、频域能量和相同，符合帕塞瓦尔定理。

综上，证明DFT变换正确。

A、补0效应研究原数组：x[8]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,},{8,0}}示例程序中补0后数组为：x2[16]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,0},{8,0},{0,0},{0,0},{0,0},{0,0} ,{0,0},{0,0},{0,0},{0,0}}补0方式我使用的补0方式为：for(i=8;i<13;i++)x2[i]=COMPLEX(0,0);补0后数组为：x2[13]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,0},{8,0},{0,0},{0,0},{0,0},{0,0} ,{0,0}}结果分析与图在时域中，信号长度增加，由于所增加的项均为零，波形仍与未补0时相同未补零时的信号时域图补5个零后的信号时域图补8个零后的信号时域图经过DFT变换后，X(k)长度也会随着x(n)长度的增加而增加，且增加的值非零未在末端补零时，信号频谱图在末端补5个零时，信号频谱图在末端补8个零时，信号频谱图可以看到，经过补0，经过DFT变换的频谱与未补零时形状基本相同，只是在长度上进行扩展，且补零数量越多，扩展越长。

可以理解为经过补0效应，增加了频域采样频率，但是由于信号未增加新的信息，因此不能提高物理分辨率。

在能量上，补5/8个零时，信号能量时域、频域能量和如下：时域能量和、频域能量和始终相等，符合帕塞瓦尔定理，且能量与未插值时的相同。

B、插值效应研究原数组：x[8]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,},{8,0}}示例程序中插值后数组为：x3[16]={{1,0},{8,0},{2,0},{7,0},{3,0},{6,0},{4,0},{5,0},{5,0},{4,0},{6,0},{3,0} ,{7,0},{2,0},{8,0},{1,0}}插值方式我使用的插值方式为：for(i=0;i<16;i=i+2){x3[i]=COMPLEX(1+i/2,0);x3[i+1]=COMPLEX(i*0.5+2.5,0);}插值后数组为：x[16]={{1,0},{3,0},{2,0},{4,0},{3,0},{5,0},{4,0},{6,0},{5,0},{7,0},{6,0},{8,0}, {7,0},{9,0},{8,0},{10,0}}结果分析与图(1)在示例程序中，在x[8]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,0},{8,0}}中反向插入原序列，使原序列变为x3[16]={{1,0},{8,0},{2,0},{7,0},{3,0},{6,0},{4,0},{5,0},{5,0},{4,0},{6,0},{3,0} ,{7,0},{2,0},{8,0},{1,0}}，再进行DFT变换，观察频谱，对比时域、频域能量和。

反向插值后，时域、频域图可以看到，反向插值后，信号频谱有了很大的直流分量，且近乎左右对称。

从三维频谱图上可以看出，高频、低频部分实际上是共轭反对称：反向插值后，三维频域图。

符合帕塞瓦尔定理，且能量是未插值时的2倍。

(2)在x[8]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,0},{8,0}}中插入序列{{3,0},{4,0},{5,0},{6,0},{7,0},{8,0},{9,0},{10,0}}，使原序列变为x3[16]= {{1,0},{3,0},{2,0},{4,0},{3,0},{5,0},{4,0},{6,0},{5,0},{7,0},{6,0},{8,0},{7,0}, {9,0},{8,0},{10,0}}，再进行DFT变换，观察频谱，对比时域、频域能量和。

插值后，时域、频域图可以看到，插值后，信号频谱有了很大的直流分量，且近乎左右对称。

，符合帕塞瓦尔定理。

(3)在x[8]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,0},{8,0}}中正向插入原序列，使原序列分别变为x2[16]={{1,0},{1,0},{2,0},{2,0},{3,0},{3,0},{4,0},{4,0},{5,0},{5,0},{6,0},{6,0},{7,0}, {7,0},{8,0},{8,0}}，再进行DFT变换，观察频谱，对比时域、频域能量和。

正向插值后，时域、频域图可以看到，正向插值后，信号频谱有了很大的直流分量，且近乎左右对称。

符合帕塞瓦尔定理，且能量是未插值时的2倍。

C、插0效应研究原数组：x[8]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,},{8,0}}示例程序中插0后数组为：x4[16]={{1,0},{0,0},{2,0},{0,0},{3,0},{0,0},{4,0},{0,0},{5,0},{0,0},{6,0},{0,0} ,{7,0},{0,0},{8,0},{0,0}}插0方式我使用的插0方式为：for(i=0;i<16;i=i+3){x4[i]=COMPLEX(1+i/2,0);x4[i+1]=COMPLEX(2+i/2,0);x4[i+2]=COM PLEX(0,0);}插0后数组为：x4[12]={{1,0},{2,0},{0,0},{3,0},{4,0},{0,0},{5,0},{6,0},{0,0},{7,0},{8,0},{0,0} }结果分析与图(1)在示例程序中，在x[8]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,0},{8,0}}中，每隔一点，插入1个0值，使原序列分别变为x1[16]={{1,0},{0,0},{2,0},{0,0},{3,0},{0,0},{4,0},{0,0},{5,0},{0,0},{6,0},{0,0},{7,0}, {0,0},{8,0},{0,0}}，再进行DFT变换，观察频谱，对比时域、频域能量和。

插0前，时域、频域图插0后，时域、频域图可以看到，插0后的频谱是对原始信号频谱的周期延拓。

画出三维图像，可以更直观地看出周期延拓关系：未插入零/插入一个零后的三维频谱图通过对插零后图像进行DFT运算，可以证明插零后的DFT是原信号DFT的周期延拓。

符合帕塞瓦尔定理，且能量与未插值时的相同。

（2）在x[8]={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,0},{8,0}}中，每隔两点，插入1个0值，使原序列变为x4[16]={{1,0},{2,0},{0,0},{3,0},{4,0},{0,0},{5,0},{6,0},{0,0},{7,0},{8,0},{0,0}}，再进行DFT变换，观察频谱，对比时域、频域能量和。

插0后，时域、频域图符合帕塞瓦尔定理源程序：// 11yy.cpp : Defines the entry point for the console application.//#include "stdafx.h"#include"D:\xhclgcyy\x_math.cpp"#include"D:\xhclgcyy\x_graph.cpp"void plotgri2(COLORREF gridcolor,COLORREF linecolor,COMPLEX p[],int N){int i;HPENpen1=CreatePen(PS_SOLID,1,gridcolor),oldpen=(HPEN)SelectObject(win3.hdc,pen1); HPEN pen2=CreatePen(PS_SOLID,1,linecolor);for(i=0;i<N;i++)line2(i,0,i,abs(p[i]));SelectObject(win3.hdc,pen2);moveto2(0,p[0].r);for(i=0;i<N;i++)lineto2(i,abs(p[i]));SelectObject(win2.hdc,oldpen);DeleteObject(pen1);DeleteObject(pen2);}void plotgri3(COLORREF gridcolor,COLORREF linecolor,COMPLEX p[],int N){int i;HPENpen1=CreatePen(PS_SOLID,1,gridcolor),oldpen=(HPEN)SelectObject(win3.hdc,pen1); HPEN pen2=CreatePen(PS_SOLID,1,linecolor);for(i=0;i<N;i++)line3(i,0,0,i,p[i].r,p[i].i);SelectObject(win3.hdc,pen2);moveto3(0,p[0].r,p[0].i);for(i=0;i<N;i++)lineto3(i,p[i].r,p[i].i);SelectObject(win2.hdc,oldpen);DeleteObject(pen1);DeleteObject(pen2);}void main(){int i;double sumT,sumF;COMPLEX x[8],//{{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,},{8,0}}X[8],x2[13],//={{1,0},{2,0},{3,0},{4,0},{5,0},{6,0},{7,0},{8,0},{0,0},{0,0},{0,0 },{0,0},{0,0}}X2[16],x3[16],//={{{1,0},{1,0},{2,0},{2,0},{3,0},{3,0},{4,0},{4,0},{5,0},{5,0},{6, 0},{6,0},{7,0},{7,0},{8,0},{8,0}}}X3[16],x4[12],//={{1,0},{2,0},{0,0},{3,0},{4,0},{0,0},{5,0},{6,0},{0,0},{7,0},{8,0 },{0,0}}X4[16];//给待变换的复数数组赋值：for(i=0;i<8;i++){x[i]=COMPLEX(i+1,0);X[i]=COMPLEX(0,0);}for(i=0;i<8;i++)x2[i]=COMPLEX(i+1,0);for(i=8;i<13;i++)x2[i]=COMPLEX(0,0);for(i=0;i<16;i=i+2){x3[i]=COMPLEX(1+i/2,0);x3[i+1]=COMPLEX(1+i/2,0);}for(i=0;i<16;i=i+3){x4[i]=COMPLEX(1+i/2,0);x4[i+1]=COMPLEX(2+i/2,0);x4[i+2]=COM PLEX(0,0);}//第1步：验证dft函数正确性dft(X,x,8,1);for(i=0;i<8;i++)printf("X[%d]=%f+%f\n",i,X[i].r,X[i].i);getch();dft(x,X,8,-1);for(i=0;i<8;i++)printf("x[%d]=%f+%f\n",i,x[i].r,x[i].i);getch();for(sumT=0,sumF=0,i=0;i<8;i++){sumT=sumT+x[i].r*x[i].r;sumF=sumF+X[i].r*X[i].r+X[i].i*X[i].i;}printf("时域能量和=%f,频域能量和=%f\n",sumT,sumF/8.0);window2("函数图形显示",-20,40,20,-20,"t","f(t)");xy2(BLUE);plotgri2(BLUE,RED,X,8);getch();frame2(win2.xstr,win2.ystr);xy2(BLUE);plotgri2(BLUE,RED,x,8);getch();window3("周期信号频谱图",-1,-12,-12,20,12,12,"N","r","i");xyz3(BLUE);plotgri3(BLUE,RED,X,8);getch();frame3();xyz3(BLUE);plotgri3(BLUE,RED,x,8);getch();//第2步：补0效应dft(X2,x2,13,1);for(i=0;i<13;i++)printf("X2[%d]=%f+%f\n",i,X2[i].r,X2[i].i);getch();dft(x2,X2,13,-1);for(i=0;i<13;i++)printf("x2[%d]=%f+%f\n",i,x2[i].r,x2[i].i);getch();for(sumT=0,sumF=0,i=0;i<13;i++){sumT=sumT+x2[i].r*x2[i].r;sumF=sumF+X2[i].r*X2[i].r+X2[i].i*X2[i].i;}printf("时域能量和=%f,频域能量和=%f\n",sumT,sumF/16.0);window2("函数图形显示",-20,40,20,-20,"t","f(t)");xy2(BLUE);plotgri2(BLUE,RED,X2,13);getch();frame2(win2.xstr,win2.ystr);xy2(BLUE);plotgri2(BLUE,RED,x2,13);getch();window3("周期信号频谱图",-1,-12,-12,20,12,12,"N","r","i"); xyz3(BLUE);plotgri3(BLUE,RED,X,13);getch();frame3();xyz3(BLUE);plotgri3(BLUE,RED,x,13);getch();//第3步：插值效应dft(X3,x3,16,1);for(i=0;i<16;i++)printf("X3[%d]=%f+%f\n",i,X3[i].r,X3[i].i); getch();dft(x3,X3,16,-1);for(i=0;i<16;i++)printf("x3[%d]=%f+%f\n",i,x3[i].r,x3[i].i); getch();for(sumT=0,sumF=0,i=0;i<16;i++){sumT=sumT+x3[i].r*x3[i].r;sumF=sumF+X3[i].r*X3[i].r+X3[i].i*X3[i].i;}printf("时域能量和=%f,频域能量和=%f\n",sumT,sumF/16.0);window2("函数图形显示",-20,20,20,-20,"t","f(t)");xy2(BLUE);plotgri2(BLUE,RED,X3,16);getch();frame2(win2.xstr,win2.ystr);xy2(BLUE);plotgri2(BLUE,RED,x3,16);getch();window3("周期信号频谱图",-1,-12,-12,20,12,12,"N","r","i"); xyz3(BLUE);plotgri3(BLUE,RED,X3,16);getch();frame3();xyz3(BLUE);plotgri3(BLUE,RED,x3,16);getch();//第4步：插0效应dft(X4,x4,12,1);for(i=0;i<12;i++)printf("X4[%d]=%f+%f\n",i,X4[i].r,X4[i].i); getch();dft(x4,X4,12,-1);for(i=0;i<12;i++)printf("x4[%d]=%f+%f\n",i,x4[i].r,x4[i].i); getch();for(sumT=0,sumF=0,i=0;i<12;i++){sumT=sumT+x4[i].r*x4[i].r;sumF=sumF+X4[i].r*X4[i].r+X4[i].i*X4[i].i;}printf("时域总和=%f,频域总和=%f\n",sumT,sumF/12.0);window2("函数图形显示",-20,40,20,-20,"t","f(t)");xy2(BLUE);plotgri2(BLUE,RED,X4,12);getch();frame2(win2.xstr,win2.ystr);xy2(BLUE);plotgri2(BLUE,RED,x4,12);getch();window3("周期信号频谱图",-1,-12,-12,20,12,12,"N","r","i"); xyz3(BLUE);plotgri3(BLUE,RED,X4,12);getch();frame3();xyz3(BLUE);plotgri3(BLUE,RED,x4,12);getch();}训练十二 DFT及抽样定理研究D、给定单频信号抽样1. 对给定信号x(t)=sin(2πfct),fc==50,N=264进行抽样，抽样频率分布为100Hz，110Hz，200Hz，230Hz，250Hz。