 上传我的文档
 下载
 收藏
粉丝量：87
大部分资料来源于网络,仅供大家参考学习,版权归原作者。若有侵权，敬请留言告知，本人会及时删除侵权文档,谢谢!
 下载此文档
数字信号处理课后答案第6章【PPT】
下载积分：2000
内容提示：数字信号处理课后答案第6章【PPT】
文档格式：PPT|
浏览次数：160|
上传日期： 07:34:12|
文档星级：
全文阅读已结束，如果下载本文需要使用
 2000 积分
下载此文档
该用户还上传了这些文档
数字信号处理课后答案第6章【PPT】
关注微信公众号使用Origin进行函数绘图的方法 - Office办公助手
最专业的Office办公软件学习和原创Office教程网站！
当前位置： >
使用Origin进行函数绘图的方法
发布时间：11-21　　来源：Office办公助手
函数绘图的特点是不需要知道具体的数据，只需知道x,y之间的函数关系即可，且操作简单。本文以Origin为例，图文讲解了使用Origin进行函数绘图的方法。
1、打开Origin 8，进入&File(文件)&菜单，点击&New(新建)&。
2、在New(新建)对话框中，选择Function(函数)，点击OK按钮。
3、在弹出的函数编辑对话框中，编写所需要的函数5*exp(-x/3)。注意：不需要写等号!
点击OK按钮，即可得到该函数的曲线。
4、在生成的图片中，双击纵坐标的数字，进入坐标体系的编辑页面。
5、选择&Scale&，点击左侧的&Vertical&，对纵坐标进行编辑。
From(坐标起始值) To(坐标终止值)
Increment(坐标增量，表示每个大刻度坐标的增加值)
# Minor(两个大刻度之间，小刻度的数量)
6、按照步骤4-5的方法，对横坐标进行修改。
7、坐标修改完后，得到下图。然后点击&New Function&按钮，增加新的函数。
8、同步骤3，在新建的函数F2(x)中， 编写所需要的函数F1(x)*sin(2*x)。点击OK按钮完成。
9、继续新建函数F3(x)，编写所需要的函数-F1(x)。点击OK按钮完成。
10、完成后，对曲线样式进行修改。选中要修改的曲线，通过工具栏(红色方框标记)修改曲线样式。
11、进行&图示&编辑，鼠标右键点击&图示&，选择&Properties...&，调出图示编辑对话框。
12、按下图进行编辑，上面为编辑区，下面为显示预览区。
13、编辑完成后，得到最终的效果图。完成!
上一篇： 下一篇：
Origin绘图教程
Origin软件应用
专题教程推荐
(window.slotbydup=window.slotbydup || []).push({
id: '3577897',
container: s,
size: '300,450',
display: 'inlay-fix'
热门资源下载他的最新文章
他的热门文章
您举报文章：
举报原因：
原文地址：
原因补充：
(最多只允许输入30个字)matlab filter函数使用的滤波器是什么滤波器？_百度知道
matlab filter函数使用的滤波器是什么滤波器？
FIR滤波器还是IIR滤波器？如何用Verilog实现
我有更好的答案
Fs,2*fsu/Fs];rp=0;20;fm2=fc2/10;;subplot(2,1,2)说明。（3）为了滤波器阶数最低，选用椭圆滤波器。（之后，再依次实现巴特沃斯、切比雪夫1、切比雪夫2数字滤波器）按照图2 所示的程序框图编写的实验程序为exp1.m。2、实验程序清单%实验1程序exp1,T;tplot(y2t;xt2=cos(2*pi*fm2*t);Tp=N*T;t=0:T:(N-1)*T;k=0:N-1;f=k/Tp,rs,wp); %调用ellip计算椭圆带通DF系统函数系数向量B和Ay2t=filter(B;plot(t,xn);xlabel('t/s');ylabel(yn).1,Tp/8,min(st),max(st)]);title('(a) s(t)的波形')subplot(3,1,2)stem(f,abs(fxt)/max(abs(fxt)),'.');title('(b) s(t)的频谱')axis([0,Fs/5,0,1.2]);xlabel('f/Hz');ylabel('幅度');.*cos(2*pi*fc2*t);xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t),rs);;subplot(2,1,2),rs,Fs].1,ws,rp,wp]=ellipord( %调用绘图函数tplot绘制滤波器输出波形%高通滤波器设计与实现================================================fp=890;fs=450;wp=2*fp/Fs;t=n*T; %调用ellip计算椭圆带通DF系统函数系数向量B和Ay3t=filter(B,A,st);
%滤波器软件实现% 高低通滤波器设计与实现绘图部分figure(4);subplot(2,1,1);myplot(B,A); %调用绘图函数tplot绘制滤波器输出波形%带通滤波器设计与实现====================================================fpl=440;fpu=560;fsl=275;fsu=900;wp=[2*fpl/Fs;ws=2*fs/myplot(B,A);Fs;
%采样频率%调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st st=%低通滤波器设计与实现=========================================fp=280;损耗函数曲线');function tplot(xn,T,yn)%时域序列连续曲线绘图函数%tplot(y3t;;subplot(2,1,2):绘图信号的纵坐标名称（字符串）% T为采样间隔n=0:length(xn)-1,t(end),min(xn),1.2*max(xn)])
function st=mstgN=2000Fs=10000;T=1/Fs;
%滤波器软件实现% 带通滤波器设计与实现绘图部分figure(3);subplot(2,1,1);tplot(y1t,A]=ellip(N,rp,wp); %调用绘图函数tplot绘制滤波器输出波形
function myplot(B; %调用ellipord计算椭圆DF阶数N和通带截止频率wp[B,A]=ellip(N;)axis([0,1,-80,5]);title(&#39,T,A]=ellip(N,pi');ylabel('幅度(dB)&#39.1,st);ws=2*fs/axis([0;
%调用ellipord计算椭圆DF阶数N和通带截止频率wp[B.m% IIR数字滤波器设计及软件实现
%调用绘图函数myplot绘制损耗函数曲线yt='y_1(t)'xlabel('rs=60;),rs);
%调用绘图函数myplot绘制损耗函数曲线yt='y_3(t)'close allFs=10000;10;xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t);plot(W/pi,yt);xlabel('t/s');ylabel('s(t)'rs=60,A,A)%时域离散系统损耗函数绘图%B为系统函数分子多项式系数向量%A为系统函数分母多项式系数向量[H;grid on,ws,rp,yt););rs=60,ws,
%DF指标（低通滤波器的通、阻带边界频）[N,W]=freqz(B,A,1000);m=abs(H),wp]=ellipord(wp,st);
%调用ellipord计算椭圆DF阶数N和通带截止频率wp[B,rs);Fs;rp=0;
%调用绘图函数myplot绘制损耗函数曲线yt='y_2(t)&#39,st);T=1&#47,yt),N);subplot(3,1,1)plot(t,rp,fm1=fc1&#47,2*fpu&#47,20*log10(m/max(m)));ws=[2*fsl/
%滤波器软件实现% 低通滤波器设计与实现绘图部分figure(5);subplot(2,1,1);
%调用ellip计算椭圆带通DF系统函数系数向量B和Ay1t=filter(B;Fs;rp=0,A;fc3=Fs/40。（2）与信号产生函数mstg相同，采样频率Fs=10kHz;st=xt1+xt2+xt3;fxt=fft(st,wp]=ellipord(fm3=fc3/\omega/\ [N:信号数据序列，myplot(B,A);axis([0;10;fc2=Fs/fc1=Fs/10;fs=600;wp=2*fp/Fs,'high&#39,T：（1）为了使滤波器阶数尽可能低;
%DF指标（低通滤波器的通、阻带边界频）[N，每个滤波器的边界频率选择原则是尽量使滤波器过渡带宽尽可能宽
为您推荐：
其他类似问题
matlab的相关知识
换一换
回答问题，赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。文档分类：
下载后只包含 1 个 DOC 格式的文档，没有任何的图纸或源代码，
下载前请先预览，预览内容跟原文是一样的，在线预览图片经过高度压缩，下载原文更清晰。
您的浏览器不支持进度条
下载文档到电脑，查找使用更方便
还剩?页未读，继续阅读
播放器加载中，请稍候...
该用户其他文档
下载所得到的文件列表《数字信号处理(第三版)》第10章实验程序(验证版)(西安电子科技大学出版_主编：高西全_丁玉美)介绍.doc
文档介绍：
《数字信号处理(第三版)》第10章实验程序
注意:每次实验前,请新建以下8个m文件并依实验需要(已注明)移动至matlab当前工作目录下,便于主程序调用。各程序均已验证,请直接复制即可。
文件名:tstem.m(实验一、二需要)
function tstem(xn,yn)
%时域序列绘图函数
%xn:被绘图的信号数据序列,yn:绘图信号的纵坐标名称(字符串)
n=0:length(xn)-1;
stem(n,xn,'.');
xlabel('n');ylabel('yn');
axis([0,n(end),min(xn),1.2*max(xn)]);
文件名:tplot.m(实验一、四需要)
function tplot(xn,T,yn)
%时域序列连续曲线绘图函数
%xn:信号数据序列,yn:绘图信号的纵坐标名称(字符串)
%T为采样间隔
n=0;length(xn)-1;t=n*T;
plot(t,xn);
xlabel('t/s');ylabel(yn);
axis([0,t(end),min(xn),1.2*max(xn)]);
文件名:myplot.m(实验一、四需要)
%(1)计算时域离散系统损耗函数并绘制曲线图。
function myplot(B,A)
%B为系统函数分子多项式系数向量
%A为系统函数分母多项式系数向量
[H,W]=freqz(B,A,1000)
plot(W/pi,20*log10(m/max(m)));
xlabel('\omega/\pi');ylabel('幅度(dB)')
axis([0,1,-80,5]);title('损耗函数曲线');
文件名:mstem.m(实验一、三需要)
function mstem(Xk)
%mstem(Xk)绘制频域采样序列向量Xk的幅频特性图
M=length(Xk);
k=0:M-1;wk=2*k/M;%产生M点DFT对应的采样点频率(关于pi归一化值)
stem(wk,abs(Xk),'.');%绘制M点DFT的幅频特性图
xlabel('w/\pi');ylabel('幅度');
axis([0,2,0,1.2*max(abs(Xk))]);
文件名:mpplot.m(实验一需要)
%(2)计算时域离散系统损耗函数和相频特性函数,并绘制曲线图。
function mpplot(B,A,Rs)
%mpplot(B,A,Rs)
%时域离散系统损耗函数和相频特性绘图
%B为系统函数分子多项式系数向量
%A为系统函数分母多项式系数向量
%Rs为滤波器阻带最小衰减,省略则幅频曲线最小值取-80dB
if nargin&3 ymin=-80;else ymin=-Rs-20;%确定幅频曲线纵坐标最小值
[H,W]=[H,W]=freqz(B,A,1000)
subplot(2,2,1);
plot(W/pi,20*log10(m/max(m)));
xlabel('\omega/\pi');ylabel('幅度(dB)')
axis([0,1,ymin,5]);title('损耗函数曲线');
subplot(2,2,3);
plot(W/pi,p/pi);
xlabel('\omega/\pi');ylabel('相位/\pi');
title('(b)相频特性曲线');
文件名:mfftplot.m(实验一需要)
function mfftplot(xn,N)
%mfftplot(xn,N)计算序列向量xn的N点fft并绘制其幅频特性曲线
Xk=fft(xn,N);%计算信号xn的频谱的N点采样
%===以下为绘图部分====
k=0:N-1;wk=2*k/N;
m=abs(Xk);mm=max(m);
plot(wk,m/mm);
xlabel('\omega/\pi');ylabel('幅度(dB)');
axis([0,2,0,1.2]);
title('幅度特性曲线');
文件名:mstg.m(实验四需要)
function st=mstg
%产生信号序列向量st,并显示st的时域波形和频谱
%st=mstg 返回三路调幅信号相加形成的混合信号,长度N=1600
N=1600; %N为信号st的长度
Fs=10000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz,Tp为采样时间
t=0:T:(N-1)*T;k=0:N-1;f=k/Tp;
fc1=Fs/10; %第1路调幅信号的载波频率fc1=1000Hz
fm1=fc1/10; %第1路调幅信号的调制信号频率fm1=100Hz
fc2=Fs/20; %第2路调幅信号的载波频率fc2=500Hz
fm2=fc2/10; %第2路调幅信号的调制信号频率fm2=50Hz
fc3=Fs/40; %第3路调幅信号的载波频率fc3=250Hz
fm3=fc3/10; %第3路调幅信号的调制信号频率fm3=25Hz
xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); %产生第1路调幅信号
xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); %产生第2路调幅信号
xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %产生第3路调幅信号
st=xt1+xt2+xt3; %三路调幅信号相加
fxt=fft(st,N); %计算信号st的频谱
%====以下为绘图部分,绘制st的时域波形和幅频特性曲线====
subplot(3,1,1);
plot(t,st);xlabel('t/s');ylabel('s(t)');
axis([0,Tp/8,min(st),max(st)]);title('(a) s(t)的波形');
subplot(3,1,2);
stem(f,abs(fxt)/max(abs(fxt)),'.');title('(b) s(t)的频谱');
axis([0,Fs/5,0,1.2]);
xlabel('f/Hz');ylabel('幅度');
文件名:xtg.m(实验五需要)
function xt=xtg(N)
%实验五信号x(t)产生,并显示信号的幅频特性曲线
%xt=xtg(N) 产生一个长度为N,有加性高频噪声的单频调幅信号xt,采样频率Fs=1000Hz
%载波频率fc=Fs/10=100Hz,调制正弦波频率f0=fc/10=10Hz.
Fs=1000;T=1/Fs;Tp=N*T;
内容来自淘豆网www.taodocs.com转载请标明出处.
文件大小：87 KB
下载次数：