CCS查看DSP程序运行时间(精)
CCS查看DSP程序运行时间
1．进入CCS环境，装载已有工程，并load生成的.out文件，并找到要察看代码执行周期的代码处。如图1所示。
图1
1．选择ccs的菜单Profiler中的enable clock,如图所示。
图2
2．选择Profiler菜单下的clock setup子菜单，并在Instruction Cycle中输入你的DSP时钟周期，它的单位为纳秒，例如，2407的系统时钟为40MHz,你就该
填入25,如果是2812系统时钟为150MHz,就该填入6.67ns，其他配置不动，然后确定。如图3所示。
图3
3．选择Profiler菜单下的Start New Session子菜单，出现如图4所示的对话框，可以改名字，也可以不改，本例中不修改，直接确定。
图4
4．通过第四部设定后就出现了如图5所示的一个窗体。
这个窗体中，有四个选项卡，其中Files为以源文件列出统计数据，Functions选项卡用于剖析程序中的函数，Ranges用于剖析一段连续的代码，Setup用于设置开始点和结束点，用于剖析不连续的代码。
窗体的左边按钮的含义为：（这里介绍主要的）
剖析所有的函数。
建立剖析区域。
设置开始点。
设置结束点。
在窗体中剖析数据有一个表格，用红框圈起来的，每个表格的字段名的含义为：Code size:剖析代码的大小，以程序存储器最小可寻址单元为单位，此值在剖析过程中不会发生变化。
Incl. Count:在统计过程中，程序运行进入剖析代码段的次数
Incl. Total: 在统计工程中剖析代码段消耗的所有时钟周期（如果是统计时钟周期的话，CCS还可以统计子程序调用等其他计数，统计其他特性则显示相应的值）。Incl. Maximum: 执行剖析代码段一遍（包括在剖析代码段中对子程序的调用）消耗的最大时钟周期（由于每次进入剖析代码段的初始条件不同等原因，每次运行剖析代码段消耗的时钟周期可能不同）；
Incl. Minimum: 执行剖析代码段一遍（包括在剖析代码段中对子程序的调用）消耗的最小时钟周期
Incl. Average: 剖析代码段执行一遍（包括在剖析代码段中对子程序的调用）消耗的平均时钟周期。--
以上这三个就是用户关心的代码执行的时钟周期。
Excl. Count:在统计过程中，程序运行进入剖析代码段的次数，与Incl.Count的值相同。
Excl. Maximum: 剖析代码段执行一遍（不包括在剖析代码段中对子程序的调用）消耗的最大时钟周期。
Excl. Minimum: 剖析代码段执行一遍（不包括在剖析代码段中对子程序的调用）消耗的最小时钟周期。
Excl. Average: 剖析代码段执行一遍（不包括在剖析代码段中对子程序的调用）消耗的平均时钟周期。
5．以剖析函数为例，找到该函数，然后将光标放在该函数的函数名上，选择建立剖析区域按钮，图中用红框框起来的那个按钮。如图6所示。
图六。
6．出现对话框，如图7所示，因为我们做的是function，所以不用修改，如果
做的是一段代码，只要把下拉菜单里的function改成Range即可。
图7
7．点击OK后出现，如图8所示，各个字段已经被赋予了初值。
图8
8．接下来运行程序RUN,就可以剖析出你所选中的代码的执行周期了。如图所示。且这些值是随着程序运行的时间而变化的，动态显示。
图9
值得注意的是，图中显示的是时钟周期，不是时间，要看时间的化，用时钟周期乘以前面设定的时钟周期ns值，如6.67ns，就是最终函数执行的时间了
数字信号处理实验一
一、实验目的 1. 通过本次实验回忆并熟悉MATLAB这个软件。 2. 通过本次实验学会如何利用MATLAB进行序列的简单运算。 3. 通过本次实验深刻理解理论课上的数字信号处理的一个常见方法——对时刻n的样本附近的一些样本求平均，产生所需的输出信号。 3. 通过振幅调制信号的产生来理解载波信号与调制信号之间的关系。 二、实验内容 1. 编写程序在MATLAB中实现从被加性噪声污染的信号中移除噪声的算法，本次试验采用三点滑动平均算法，可直接输入程序P1.5。 2. 通过运行程序得出的结果回答习题Q1.31-Q1.33的问题，加深对算法思想的理解。 3. 编写程序在MATLAB中实现振幅调制信号产生的算法，可直接输入程序P1.6。 4. 通过运行程序得出的结果回答习题Q1.34-Q1.35的问题，加深对算法思想的理解。 三、主要算法与程序 1. 三点滑动平均算法的核心程序： %程序P1.5 %通过平均的信号平滑 clf; R=51; d=0.8*(rand(R,1)-0.5);%产生随噪声 m=0:R-1; s=2*m.*(0.9.^m);%产生为污染的信号 x=s+d';%产生被噪音污染的信号 subplot(2,1,1); plot(m,d','r-',m,s,'g--',m,x,'b-.');
xlabel('时间序号n');ylabel('振幅'); legend('d[n]','s[n]','x[n]'); x1=[0 0 x];x2=[0 x 0];x3=[x 0 0]; y=(x1+x2+x3)/3; subplot(2,1,2); plot(m,y(2:R+1),'r-',m,s,'g--'); legend('y[n]','s[n]'); xlabel('时间序号n');ylabel('振幅'); 2. 振幅调制信号的产生核心程序：（由于要几个结果，因此利用subplot函数画图） %程序P1.6 %振幅调制信号的产生 n=0:100; m=0.1;fH=0.1;fL=0.01; m1=0.3;fH1=0.3;fL1=0.03; xH=sin(2*pi*fH*n); xL=sin(2*pi*fL*n); y=(1+m*xL).*xH; xH1=sin(2*pi*fH1*n); xL1=sin(2*pi*fL1*n); y1=(1+m1*xL).*xH; y2=(1+m*xL).*xH1; y3=(1+m*xL1).*xH; subplot(2,2,1); stem(n,y); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.1;fH=0.1;fL=0.01;'); subplot(2,2,2); stem(n,y1); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.1;fL=0.01;'); subplot(2,2,3); stem(n,y2); grid; xlabel('时间序号n');ylabel('振幅');title('m=0.3;fH=0.3;fL=0.01;'); subplot(2,2,4); stem(n,y3); grid;
显示系统时间运行时间的程序详解
S7200+TD400C显示系统时间、运行时间的程序详解 楼主??发帖时间：2007-7-21 15:46:00 ?????????? 看见论坛上有些朋友对西门子TD文本显示器显示时钟的问题比较关心，在这个帖子里笔者给出一个已经应用于工程上的程序例子，并作出详细分析，希望对关心这个问题的朋友有些帮助。 ????这个程序是S7200+TD400C显示系统时间、当班运行时间、累计运行时间的例子。 ??1楼回复时间：2007-7-21 16:02:00 系统硬件配置如下： ????PLC：西门子S7-200?CN；CPU?226?CN?REL?02.01；AC100~230V电源/DC24V 输入/继电器输出；订货号6ES7?216-2BD23-0XB8；固件02.01?Build?2；ASIC：01.00。 ????文本显示器：TD400C；订货号6AV6 6640-0AA00-0AX0；自带9芯TD/CPU电缆；版本：1.0.0.3。 ????S7-200与TD400C通过TD400C随机配置的TD/CPU通信电缆连接，实现电源供给和通信（因为当TD400C与S7-200?CPU之间的距离小于?2.5米时，采用 TD/CPU电缆的方式进行供电；当TD400C与S7-200?CPU之间的距离超过2.5米时，使用外部电源供电并使用PROFIBUS组件连接网络）。 ??2楼回复时间：2007-7-21 16:12:00 系统软件配置如下： ????S7-200参数设置：在“系统块”的“断电数据保持”中设置VW1600以后2000个单元为断电数据保持。 ????TD400C参数设置： ????????TD400C地址：1 ????????CPU地址：2 ????????参数块地址：0 ????????波特率：9.6K ????????HSA：31 ????????GUF：10 ????????键盘声音反馈：开 ????????屏幕保护时间：10分钟 S7-200的系统时钟调整为准确的北京时间。 ??3楼回复时间：2007-7-21 16:15:00 下面先写一下该程序所用到的存储器的意义，以便于理解程序。 本班运行时间：小时VW1600、分钟VW1604、秒VW1608； 累计运行时间：小时VW1620、分钟VW1624、秒VW1628； 系统时间（BCD码字节）：年VB2000、月VB2001、日VB2002、小时VB2003、分钟VB2004、秒VB2005、星期VB2007； 系统时间（整数）：?????年VW1644、月VW1648、日VW1652、小时VW1656、分钟VW1660、秒VW1664； 本班设备开始运行时间：小时VW1680、分钟VW1684、秒VW1688； 本班设备最后运行时间：小时VW1700、分钟VW1704、秒VW1708； 上班时间设置：小时VW1720、分钟VW1724；
Matlab中计算程序运行时间的三种方法
Matlab中计算程序运行时间的三种方法 经常我们需要计算我们程序到底运行多长时间，这样可以比较程序的执行效率。当然这个对于只有几秒钟的小程序没有什么意义，但是对于大程序就有很重要的意义了。 下面我们就说说Matlab中计算程序运行时间的三种常用方法吧！ 注意：三种方法由于使用原理不一样，得到结果可能有一定的差距！ 1、tic和toc组合(使用最多的) 计算tic和toc之间那段程序之间的运行时间，它的经典格式为 1. tic 2. 。。。。。。。。。。 3. toc 复制代码 换句话说程序，程序遇到tic时Matlab自动开始计时，运行到toc时自动计算此时与最近一次tic之间的时间。这个有点拗口，下面我们举个例子说明 1. % by dynamic of Matlab技术论坛 2. % see also https://www.360docs.net/doc/949935034.html, 3. % contact me matlabsky@https://www.360docs.net/doc/949935034.html, 4. % 2009-08-18 12:08:47 5. clc 6. tic;%tic1 7. t1=clock; 8. for i=1:3 9. tic ;%tic2 10. t2=clock; 11. pause(3*rand) 12. %计算到上一次遇到tic的时间，换句话说就是每次循环的时间 13. disp(['toc计算第',num2str(i),'次循环运行时间：',num2str(toc)]); 14. %计算每次循环的时间 15. disp(['etime计算第',num2str(i),'次循环运行时间：',num2str(etime(clock,t2))]); 16. %计算程序总共的运行时间 17. disp(['etime计算程序从开始到现在运行的时间:',num2str(etime(clock,t1))]); 18. disp('======================================') 19. end 20. %计算此时到tic2的时间，由于最后一次遇到tic是在for循环的i=3时，所以计算 的是最后一次循环的时间 21. disp(['toc计算最后一次循环运行时间',num2str(toc)]) 22. disp(['etime程序总运行时间：',num2str(etime(clock,t1))]); 复制代码 运行结果如下，大家可以自己分析下 1. toc计算第1次循环运行时间： 2.5628 2. etime计算第1次循环运行时间：2.562
数字信号处理实验二
实验报告（本科） 学号 2015141443002 姓名柏冲 专业通信工程 日期 2017/12/4 实验题目时域采样和频域采样 一、实验目的
时域采样理论与频域采样理论是数字信号处理中重要的理论。要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使得采样后的信号不丢失信息；要求掌握频率采样会引起时域周期化的概念，以及频域采样定理及其对频域采样点数选择的指导作用。 二、实验过程 附：源程序 (1)时域采样 Tp=64/1000; %观察时间Tp=64毫秒 %产生M长采样序列x(n) Fs=1000; T=1/Fs; M=Tp*Fs; n=0:M-1; A=444.128; a=pi*50*2^0.5; omega=pi*50*2^0.5; xnt=A*exp(-a*n*T).*sin(omega*n*T); Xk=T*fft(xnt,M); %M点FFT[(xnt)] subplot(3,2,1); stem(xnt,'.'); %调用编绘图函数stem绘制序列图 box on;title('(a) Fs=1000Hz'); k=0:M-1;fk=k/Tp; subplot(3,2,2);stem(fk,abs(Xk),'.');title('(a) T*FT[xa(nT)],Fs=1000Hz'); xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))]); % Fs=300Hz和 Fs=200Hz的程序与上面Fs=1000Hz完全相同。 Tp=64/1000; %观察时间Tp=64毫秒 %产生M长采样序列x(n) Fs=300; T=1/Fs; M=Tp*Fs; n=0:M-1; A=444.128; a=pi*50*2^0.5; omega=pi*50*2^0.5; xnt=A*exp(-a*n*T).*sin(omega*n*T); M1=fix(M); Xk=T*fft(xnt,M1); %M点FFT[(xnt)] subplot(3,2,3); stem(xnt,'.'); %调用自编绘图函数stem绘制序列图 box on;title('(b) Fs=300Hz'); k=0:M-1;fk=k/Tp; subplot(3,2,4);stem(fk,abs(Xk),'.');title('(b) T*FT[xa(nT)],Fs=300Hz'); xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))]); Tp=64/1000; %观察时间Tp=64毫秒 %产生M长采样序列x(n) Fs=200; T=1/Fs; M=Tp*Fs; n=0:M-1; A=444.128; a=pi*50*2^0.5; omega=pi*50*2^0.5; xnt=A*exp(-a*n*T).*sin(omega*n*T); M2=fix(M);
单片机C延时时间怎样计算
C程序中可使用不同类型的变量来进行延时设计。经实验测试，使用unsigned char类型具有比unsigned int更优化的代码，在使用时 应该使用unsigned char作为延时变量。以某晶振为12MHz的单片 机为例，晶振为12M H z即一个机器周期为1u s。一. 500ms延时子程序 程序: void delay500ms(void) { unsigned char i,j,k; for(i=15;i>0;i--) for(j=202;j>0;j--) for(k=81;k>0;k--); } 计算分析: 程序共有三层循环 一层循环n:R5*2 = 81*2 = 162us DJNZ 2us 二层循环m:R6*(n+3) = 202*165 = 33330us DJNZ 2us + R5赋值 1us = 3us 三层循环: R7*(m+3) = 15*33333 = 499995us DJNZ 2us + R6赋值 1us = 3us
循环外: 5us 子程序调用 2us + 子程序返回2us + R7赋值 1us = 5us 延时总时间 = 三层循环 + 循环外 = 499995+5 = 500000us =500ms 计算公式:延时时间=[(2*R5+3)*R6+3]*R7+5 二. 200ms延时子程序 程序: void delay200ms(void) { unsigned char i,j,k; for(i=5;i>0;i--) for(j=132;j>0;j--) for(k=150;k>0;k--); } 三. 10ms延时子程序 程序: void delay10ms(void) { unsigned char i,j,k; for(i=5;i>0;i--) for(j=4;j>0;j--) for(k=248;k>0;k--);
c++程序计时
C++运行时间的代码 如何获取代码运行时间在调试中，经常需要计算某一段代码的执行时间，下面给出两种常用的方式： 第一种：使用GetTickCount函数 #include CCS查看DSP程序运行时间 2008-09-04 19:35 1．进入CCS环境，装载已有工程，并load生成的.out文件，并找到要察看代码执行周期的代码处。如图1所示。 图1 1．选择ccs的菜单Profiler中的enable clock,如图所示。 图2 2．选择Profiler菜单下的clock setup子菜单，并在Instruction Cycle中输入你的DSP时钟周期，它的单位为纳秒，例如，2407的系统时钟为40MHz,你就该填入25,如果是2812系统时钟为150MHz,就该填入6.67ns，其他配置不动，然后确定。如图3所示。 图3 3．选择Profiler菜单下的Start New Session子菜单，出现如图4所示的对话框，可以改名字，也可以不改，本例中不修改，直接确定。 图4 4．通过第四部设定后就出现了如图5所示的一个窗体。 这个窗体中，有四个选项卡，其中Files为以源文件列出统计数据，Functions选项卡用于剖析程序中的函数，Ranges 用于剖析一段连续的代码，Setup用于设置开始点和结束点，用于剖析不连续的代码。 窗体的左边按钮的含义为：（这里介绍主要的） 剖析所有的函数。 建立剖析区域。 设置开始点。 设置结束点。 在窗体中剖析数据有一个表格，用红框圈起来的，每个表格的字段名的含义为： Code size:剖析代码的大小，以程序存储器最小可寻址单元为单位，此值在剖析过程中不会发生变化。 Incl. Count:在统计过程中，程序运行进入剖析代码段的次数 Incl. Total: 在统计工程中剖析代码段消耗的所有时钟周期（如果是统计时钟周期的话，CCS还可以统计子程序调用等其他计数，统计其他特性则显示相应的值）。 Incl. Maximum: 执行剖析代码段一遍（包括在剖析代码段中对子程序的调用）消耗的最大时钟周期（由于每次进入剖析代码段的初始条件不同等原因，每次运行剖析代码段消耗的时钟周期可能不同）； Incl. Minimum: 执行剖析代码段一遍（包括在剖析代码段中对子程序的调用）消耗的最小时钟周期 Incl. Average: 剖析代码段执行一遍（包括在剖析代码段中对子程序的调用）消耗的平均时钟周期。-- 以上这三个就是用户关心的代码执行的时钟周期。 Excl. Count:在统计过程中，程序运行进入剖析代码段的次数，与Incl.Count的值相同。 Excl. Maximum: 剖析代码段执行一遍（不包括在剖析代码段中对子程序的调用）消耗的最大时钟周期。 Excl. Minimum: 剖析代码段执行一遍（不包括在剖析代码段中对子程序的调用）消耗的最小时钟周期。 Excl. Average: 剖析代码段执行一遍（不包括在剖析代码段中对子程序的调用）消耗的平均时钟周期。 5．以剖析函数为例，找到该函数，然后将光标放在该函数的函数名上，选择建立剖析区域按钮，图中用红框框起来的那个按钮。如图6所示。 图六。 6．出现对话框，如图7所示，因为我们做的是function，所以不用修改，如果做的是一段代码，只要把下拉菜单里的function改成Range即可。 武汉工程大学 数字信号处理实验报告 姓名：周权 学号：1204140228 班级：通信工程02 一、实验设备 计算机，MATLAB语言环境。 二、实验基础理论 1.序列的相关概念 2.常见序列 3.序列的基本运算 4.离散傅里叶变换的相关概念 5.Z变换的相关概念 三、实验内容与步骤 1.离散时间信号（序列）的产生 利用MATLAB语言编程产生和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形表示。 四实验目的 认识常用的各种信号，理解其数字表达式和波形表示，掌握在计算机中生成及绘制数字信号波形的方法，掌握序列的简单运算及计算机实现与作用，理解离散时间傅里叶变换，Z变换及它们的性质和信号的频域分 实验一离散时间信号（序列）的产生 代码一 单位样值 x=2; y=1; stem(x,y); title('单位样值 ') 单位阶跃序列 n0=0; n1=-10; n2=10; n=[n1:n2]; x=[(n-n0)>=0]; stem(n,x); xlabel('n'); ylabel('x{n}'); title('单位阶跃序列'); 实指数序列 n=[0:10]; x=(0.5).^n; stem(n,x); xlabel('n'); ylabel('x{n}'); title('实指数序列'); 正弦序列 n=[-100:100]; x=2*sin(0.05*pi*n); stem(n,x); xlabel('n'); ylabel('x{n}'); title('正弦序列'); 随机序列 n=[1:10]; x=rand(1,10); subplot(221); stem(n,x); xlabel('n'); ylabel('x{n}'); title('随机序列'); 内容： Q:如何获取时间？精度如何？ A: 1 使用time_t time( time_t * timer ) 精确到秒 计算时间差使用double difftime( time_t timer1, time_t timer0 ) 2 使用clock_t clock() 得到的是CPU时间精确到1/CLOCKS_PER_SEC秒 3 使用DWORD GetTickCount() 得到的是系统运行的时间精确到毫秒 4 如果使用MFC的CTime类，可以用CTime::GetCurrentTime() 精确到秒 5 要获取高精度时间，可以使用 BOOL QueryPerformanceFrequency(LARGE_INTEGER *lpFrequency)获取系统的计数器的频率 BOOL QueryPerformanceCounter(LARGE_INTEGER *lpPerformanceCount)获取计数器的值 然后用两次计数器的差除以Frequency就得到时间。 6 还有David的文章中提到的方法： Multimedia Timer Functions The following functions are used with multimedia timers. timeBeginPeriod/timeEndPeriod/timeGetDevCaps/timeGetSystemTime timeGetTime/timeKillEvent/TimeProc/timeSetEvent 精度很高 Q:GetTickCount()函数，说是毫秒记数，是真的吗，还是精确到55毫秒？ A: GetTickCount()和GetCurrentTime()都只精确到55ms(1个tick就是55ms)。如果要精确到毫秒，应该使用timeGetTime函数或QueryPerformanceCounter函数。具体例子可以参考QA001022 "VC++中使用高精度定时器"、QA001813 "如何在Windows实现准确的定时"和QA004842 "timeGetTime函数延时不准"。 Q:vc++怎样获取系统时间，返回值是什么类型的变量呢？ GetSystemTime返回的是格林威志标准时间 GetLocalTime，和上面用法一样，返回的是你所在地区的时间，中国返回的是北京时间VOID GetSystemTime( LPSYSTEMTIME lpSystemTime // address of system time structure ); 函数就可以获得了，其中LPSYSTEMTIME 是个结构体 含：年，月，日，周几，小时，分，秒，毫秒。 以下是Time的MSDN文档: Compatibility in the Introduction. Libraries LIBC.LIBSingle thread static library, retail versionLIBCMT.LIBMultithread static library, retail versionMSVCRT.LIBImport library for MSVCRT.DLL, retail version Return Value time returns the time in elapsed seconds. There is no error return. Parameter timer Storage location for time Remarks 实验一 离散时间信号分析 班级 信息131班 学号 201312030103 姓名 陈娇 日期 一、实验目的 掌握两个序列的相加、相乘、移位、反褶、卷积等基本运算。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列)}({n x 来表示，其中)(n x 代表序列的第n 个数字，n 代表时间的序列，n 的取值范围为+∞<<∞-n 的整数，n 取其它值)(n x 没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号)(t x a 进行等间隔采样，采样间隔为T ，得到)}({nT x a 一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）） （n δ、单位阶跃序列)(n u 、矩形序列)(n R N 、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 ∑∞ -∞==-= m n h n x m n h m x n y )(*)()()()( 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将)(n x 和)(n h 的变量n 换成m ，变成)(m x 和)(m h ，再将)(m h 以纵轴为对称轴反褶成)(m h -。 （2）移位：将)(m h -移位n ，得)(m n h -。当n 为正数时，右移n 位；当n 为负数时，左移n 位。 （3）相乘：将)(m n h -和)(m x 的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得)(n y 。 三、主要实验仪器及材料 微型计算机、Matlab6.5 教学版、TC 编程环境。 四、实验内容 （1）用Matlab 或C 语言编制两个序列的相加、相乘、移位、反褶、卷积等的程序； （2）画出两个序列运算以后的图形； （3）对结果进行分析； （4）完成实验报告。 五、实验结果 六、实验总结 Qt系统运行时间差计算 在网上查了很多资料，发觉网上很多用Qt写的系统运行时间差的例子写的都不是很全，今天自己研究了一下，可以成功得显示日时分秒，觉得不错，就与大家分享了 #include C与汇编程序运行时间比较 1. 程序简介 分别采用C语言和汇编语言多次对大量数据进行冒泡排序，比较并分析在数据量不同时二者所用时间的不同。 2. 流程图 3. 代码实现 #include"stdio.h" #include"stdlib.h" #include"time.h" #include"stdafx.h" //void srand (unsigned int seed); void bubble(int * d, long num); void *malloc(unsigned int size); void free(void *p); int main() { double t=0.00; int k,i,j,temp; int* p; clock_t start=0,end=0; printf("要求数组的大小：\n"); scanf("%d",&k); p=(int*)malloc(k*sizeof(int)); srand((int)time(0)); for(i=0;i 实验二离散时间系统的时域分析实验室名称: 实验时间: 六、实验记录(数据、图表、波形、程序等) Q2、1 程序代码: %产生输入信号 n = 0:100; s1 = cos(2*pi*0、05*n); %一个低频正弦 s2 = cos(2*pi*0、47*n); %一个高频正弦 x = s1+s2; %滑动平均滤波器的实现 M = input('Desired length of the filter = '); num = ones(1,M); y = filter(num,1,x)/M; clf; %显示输入与输出信号 subplot(2,2,1); plot(n, s1); axis([0, 100, -2, 2]); xlabel('Time index n'); ylabel('Amplitude'); title('Signal #1'); subplot(2,2,2); plot(n, s2); axis([0, 100, -2, 2]); xlabel('Time index n'); ylabel('Amplitude'); title('Signal #2'); subplot(2,2,3); plot(n, x); axis([0, 100, -2, 2]); xlabel('Time index n'); ylabel('Amplitude'); title('Input Signal'); subplot(2,2,4); plot(n, y); axis([0, 100, -2, 2]); xlabel('Time index n'); ylabel('Amplitude'); title('Output Signal'); axis; 显示结果: 实验三：离散LSI 系统的频域分析 一、实验内容 2、求以下各序列的z 变换： 12030() ()sin() ()sin()n an x n na x n n x n e n ωω-=== 程序清单如下： syms w0 n z a; x1=n*a^n;X1=ztrans(x1) x2=sin(w0*n);X2=ztrans(x2) x3= exp(-a*n)*sin(w0*n);X3=ztrans(x3) 程序运行结果如下： X1 =z/(a*(z/a - 1)^2) X2 =(z*sin(w0))/(z^2 - 2*cos(w0)*z + 1) X3 =(z*exp(a)*sin(w0))/(exp(2*a)*z^2 - 2*exp(a)*cos(w0)*z + 1) 3、求下列函数的逆z 变换 0 312342 1 1() () () ()() 1j z z z z X z X z X z X z z a z a z e z ω---= = = = ---- 程序清单如下： syms w0 n z a; X1=z/(z-a);x1=iztrans(X1) X2= z/(a-z)^2;x2=iztrans(X2) X3=z/ z-exp(j*w0);x3=iztrans(X3) X4=(1-z^-3)/(1-z^-1);x4=iztrans(X4) 程序运行结果如下： x1 =a^n x2 =n*a^n/a 课程名称 数字信号 实验成绩 指导教师 实 验 报 告 院系 信息工程学院 班级 学号 姓名 日期 x3 =charfcn[0](n)-iztrans(exp(i*w0),w0,n) x4 =charfcn[2](n)+charfcn[1](n)+charfcn[0](n) 4、求一下系统函数所描述的离散系统的零极点分布图，并判断系统的稳定性 （1） (0.3)()(1)(1) z z H z z j z j -= +-++ z1=[0,0.3]';p1=[-1+j,-1-j]';k=1; [b1,a1]=zp2tf(z1,p1,k); subplot(1,2,1);zplane(z1,p1); title('极点在单位圆外); subplot(1,2,2);impz(b1,a1,20); 由图可见:当极点位于单位圆内，系统的单位序列响应随着频率的增大而收敛；当极点位于单位圆上，系统的单位序列响应为等幅振荡；当极点位于单位圆外，系统的单位序列响应随着频率的增大而发散。由此可知系统为不稳定系统。 -1 -0.5 00.51 -2 -1.5-1-0.500.511.5 2Real Part I m a g i n a r y P a r t 极点在单位圆外 n (samples) A m p l i t u d e Impulse Response 计算程序运行时间time t clock t 计算程序运行时间(time_t,clock_t) 转载我们有时需要得到程序的运行时间，但我们也要知道，根本不可能精 确测量某一个程序运行的确切时间，文献中说的很明白，现摘录如下。 我们平时常用的测量运行时间的方法并不是那么精确的，换句话说，想精 确获取程序运行时间并不是那么容易的。也许你会想，程序不就是一条条指令么，每一条指令序列都有固定执行时间，为什么不好算?真实情况下，我们的计算机并不是只运行一个程序的，进程的切换，各种中断，共享的多用户，网络 流量，高速缓存的访问，转移预测等，都会对计时产生影响。 文献中还提到：对于进程调度来讲，花费的时间分为两部分，第一是计时 器中断处理的时间，也就是当且仅当这个时间间隔的时候，操作系统会选择， 是继续当前进程的执行还是切换到另外一个进程中去。第二是进程切换时间， 当系统要从进程A切换到进程B时，它必须先进入内核模式将进程A的状态保存，然后恢复进程B的状态。因此，这个切换过程是有内核活动来消耗时间的。具体到进程的执行时间，这个时间也包括内核模式和用户模式两部分，模式之 间的切换也是需要消耗时间，不过都算在进程执行时间中了。 那么有哪些方法能统计程序的运行时间呢?通过查找一些资料并结合自己的实践体会，摘录和总结了下面几种方法。 一、Linux的time命令 Linux系统下统计程序运行实践最简单直接的方法就是使用time命令，文 献[1,2]中详细介绍了time命令的用法。此命令的用途在于测量特定指令执行 时所需消耗的时间及系统资源等资讯，在统计的时间结果中包含以下数据： 实际时间(real time)：从命令行执行到运行终止的消逝时间； 用户CPU时间(user CPU time)：命令执行完成花费的系统CPU时间，即命令在用户态中执行时间的总和； 标题：如何用vc++获取系统时间和程序运行时间 出处：春天的事业 时间：Mon, 22 Jun 2009 17:34:26 +0000 作者：xiechunye 地址：https://www.360docs.net/doc/949935034.html,/read.php/612.htm 内容： Q:如何获取时间？精度如何？ A: 1 使用time_t time( time_t * timer ) 精确到秒 计算时间差使用double difftime( time_t timer1, time_t timer0 ) 2 使用clock_t clock() 得到的是CPU时间精确到1/CLOCKS_PER_SEC秒 3 使用DWORD GetTickCount() 得到的是系统运行的时间精确到毫秒 4 如果使用MFC的CTime类，可以用CTime::GetCurrentTime() 精确到秒 5 要获取高精度时间，可以使用 BOOL QueryPerformanceFrequency(LARGE_INTEGER *lpFrequency)获取系统的计数器的频率 BOOL QueryPerformanceCounter(LARGE_INTEGER *lpPerformanceCount)获取计数器的值 然后用两次计数器的差除以Frequency就得到时间。 6 还有David的文章中提到的方法： Multimedia Timer Functions The following functions are used with multimedia timers. timeBeginPeriod/timeEndPeriod/timeGetDevCaps/timeGetSystemTime timeGetTime/timeKillEvent/TimeProc/timeSetEvent 精度很高 Q:GetTickCount()函数，说是毫秒记数，是真的吗，还是精确到55毫秒？ A: GetTickCount()和GetCurrentTime()都只精确到55ms(1个tick就是55ms)。如果要精确到毫秒，应该使用timeGetTime函数或QueryPerformanceCounter函数。具体例子可以参考QA001022 "VC++中使用高精度定时器"、QA001813 "如何在Windows实现准确的定时"和QA004842 "timeGetTime函数延时不准"。 Q:vc++怎样获取系统时间，返回值是什么类型的变量呢？ GetSystemTime返回的是格林威志标准时间 GetLocalTime，和上面用法一样，返回的是你所在地区的时间，中国返回的是北京时间VOID GetSystemTime( LPSYSTEMTIME lpSystemTime // address of system time structure ); 函数就可以获得了，其中LPSYSTEMTIME 是个结构体 含：年，月，日，周几，小时，分，秒，毫秒。 以下是Time的MSDN文档: Compatibility in the Introduction. Libraries LIBC.LIBSingle thread static library, retail versionLIBCMT.LIBMultithread static library, retail versionMSVCRT.LIBImport library for MSVCRT.DLL, retail version 对输入x[n]滤波，求得y1[n] 。y[n]和y1[n]有差别吗？为什么要使用对x[n]补零后得到的x1[n]作为输入来产生y1[n] ？ 三、实验器材及软件 1?微型计算机1台 2. MATLAB 7.0 软件 四、实验原理 1.三点平滑滤波器是一个线性时不变的有限冲激响应系统，将输出延时一个抽样周期，可得到三点平滑滤波器的因果表达式，生成的滤波器表示为 y[n] 1 -(x[ n] x[n 1] x[n 2]) 3 归纳上式可得 1 M 1 y[n] x[n k] IVI k 0 此式表示了一个因果M点平滑FIR滤波器 2.对线性离散时间系统，若y1［n］和y2［n］分别是输入序列xqn］和X2［n］的响 应， 则输入 x[n] xdn] X2【n] 的输出响应为 y[n] yd n] y2【n] 则系统称为线性系统。 3.对于离散时不变系统，若y1[n]是x1[n]的响应，则输入 x[n]=x 1[n-n o] 的输出响应为 y[n]=y qn-n o] 则称系统为时不变系统。 五、实验步骤 2.1 首先利用MATLA产生一个高频正弦信号和一个低频正弦信号，利用两个信 个信号，再对输出信号进行分析。 2.2 在2.1的基础上编写num=ones[1-1]，运行程序得出结论。 2.4 分别用扫频信号通过2.1、2.2的系统，进行比较分析。 2.7 分别计算出y1[n]和y2[n]，得到yt[n]；再利用filter 函数求得y[n]，计算差值输出，比较y[n]和yt[n]。 2.20根据impz函数的调用方式，得到 num = [0.9 -0.45 0.35 0.002] den = [1 0.71 -0.46 -0.62] ，再调用impz 函数，画出图像。 2.23首先产生序列x[n]，把它作为四阶系统的输入，生成y[n]。然后将同样的输入x[n]应用到第一级得到y1[n]。接着用相同的方法得到y2[n]。最后求得两者的差，并画出图像。 2.28分别用conv函数和filter 函数求得输出，进行图像比较。 六、实验记录（数据、图表、波形、程序等） 2.1 对M=2运行上述程序，生成输入x[n]=s1[n]+s2[n] 的输出信号。输入x[n] 的哪个分量被该离散时间系统抑制? % Program P2_1 clf;CCS查看DSP程序运行时间
数字信号处理实验报告一
vc++获取系统时间和程序运行时间
数字信号处理实验一
Qt系统运行时间差计算
C与汇编程序运行时间比较
数字信号处理实验二
数字信号处理实验三
计算程序运行时间 time t clock t
如何用vc 获取系统时间和程序运行时间
数字信号处理实验二