南邮MATLAB软件设计
通信与信息工程学院
2015 / 2016 学年第 2 学期
软件设计实验报告
模块名称MATLAB软件设计
专业通信工程
学生班级B130102
学生学号B13010229
学生姓名郭智超
指导教师孟田甜
报告内容
一、要求练习的实验部分
1.题目：在时间区间 [0,10]中，绘出
t e y t
2cos 15.0--=曲线。 程序代码：
t=0:0.02:10;
y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'b'); xlabel('time');
ylabel('amplitude'); title('figure 1'); grid
实验结果：
24
6810
00.51
1.5
time
a m p l i t u d e
figure 1
2.题目：写出生成如图E2-1所示波形的MATLAB 脚本M 文件。图中虚线为正弦波，要求它的负半波被置零，且在
2
3
处被削顶。
程序代码：
t=linspace(0,10,500); %产生线性间隔的向量（在0-10之间取500点）
y=sin(t); a=sin(pi/3); z=(y>=0).*y;
z=(y>=a).*a+(y hold on; %保持当前图形 plot(t,z,'b-') xlabel('t'); ylabel('z=f(t)'); title('逐段解析函数'); legend('y=sin(t)','z=f(t)',4); %确定向量长度 axis([0 10 -1.2 1.2]); %限定坐标轴范围 hold off; grid 实验结果： 2 4 6 8 10 01t z =f (t ) 逐段解析函数 3.题目： 令??????????=987654321A ，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题： （1） sum(diag(S)>tol) - rank(A) 是多少 ？ （2） S(1,1) - norm(A) = 0 是多少 ? （3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么 ? （4） S(1,1)/S(3,3) - cond(A) 的结果是什么 ? （5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？ （6） V(:,1)'*null(A) 得到什么结果 ? （7） abs(A*null(A)) < tol 得到什么结果 ? （8） U(:,1:2) = = orth(A) 的运行结果是什么 ? 程序： A=[1 2 3;4 5 6;7 8 9] [U,S,V]=svd(A);tol=S(1,1)*3*eps sum(diag(S)>tol) - rank(A) S(1,1) - norm(A) sqrt(sum(diag(S*S))) - norm(A,'fro') S(1,1)/S(3,3) - cond(A) S(1,1)*S(2,2)*S(3,3) - det(A) < tol V(:,1)'*null(A) abs(A*null(A)) < tol U(:,1:2) == orth(A) 运行结果： （1） sum(diag(S)>tol) - rank(A) 的结果是0 （2） S(1,1) - norm(A) = 0的结果是0 （3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是3.5527e-015 （4） S(1,1)/S(3,3) - cond(A) 的结果是 -8 （5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的结果是1 （6） V(:,1)'*null(A) 的结果是 0 （7） abs(A*null(A)) < tol 的结果是1；1；1 （8） U(:,1:2) = = orth(A) 的结果是 [1 1 ；1 1 ；1 1] 4.题目：求积分dt t x ?0]sin[cos ，π20≤≤x 。 提示：abs,sin,cos,cumsum,disp,plot,grid 程序代码： t=0:pi/100:2*pi; f=abs(cos(sin(t))); F=cumsum(f)*pi/100; %求和 plot(t,F); xlabel('t'); ylabel('f'); title('积分曲线'); grid on 运行结果： 1 2 3 4 5 6 7 00.511.522.533.544.5 5t f 积分曲线 5.题目： 求方程 2,122==+xy y x 的解。提示：solve 程序代码： syms x syms y %变量定义 A=solve('x*x+y*y=1','x*y=2',x,y); Ax=A.x Ay=A.y 运行结果： Ax = -1/2*(1/2*5^(1/2)+1/2*i*3^(1/2))^3+1/4*5^(1/2)+1/4*i*3^(1/2) -1/2*(1/2*5^(1/2)-1/2*i*3^(1/2))^3+1/4*5^(1/2)-1/4*i*3^(1/2) -1/2*(-1/2*5^(1/2)+1/2*i*3^(1/2))^3-1/4*5^(1/2)+1/4*i*3^(1/2) -1/2*(-1/2*5^(1/2)-1/2*i*3^(1/2))^3-1/4*5^(1/2)-1/4*i*3^(1/2) Ay = 1/2*5^(1/2)+1/2*i*3^(1/2) 1/2*5^(1/2)-1/2*i*3^(1/2) -1/2*5^(1/2)+1/2*i*3^(1/2) -1/2*5^(1/2)-1/2*i*3^(1/2) 6.题目： 在某激励条件下，二阶系统归一化响应可表示为 ) sin(1 1)(θββ ξ+- =-t e t y t ，其中ξ为 阻尼系数， 2 1ξβ-=， ) 1( 2 ξ ξθ-=arctg 。请用不同的颜色或线型，在同一张图上，绘制 0.2:2.0:2.0=ξ取值下系统在]18,0[∈t 区间内的响应曲线，并要求用2.0=ξ和0.2=ξ对它们相应的两条曲线进行醒目的文字标志。 代码段1： clc clear syms k; t=0:0.1:18; hold on; %kz={'-r' '.g' 'b' 'c' 'm' 'y' 'k' 'w' ':r' ':g'}; kz='r g :b c m y k w r g r:'; i=0; for k=0.2:0.2:2 i=i+1; clear y; a=sqrt(1-k^2); b=atan(sqrt(1-k^2)/k); if k==1 y=(exp(t)-t-1)./exp(t); plot(t,y,kz(i)) else y=1-1/a*exp(-k*t).*sin(a*t+b); plot(t,y,kz(i)) end % plot(t,y,kz(i)) if k==0.2 text(8,1.2,'阻尼系数是0.2') end if k==2 text(14,0.9,'阻尼系数是2') end end title('二阶系统归一化响应'); xlabel('t'); ylabel('y(t)'); hold off grid 运行结果： 5 1015 20 -0.2 00.20.40.60.8 1 1.21.41.6 二阶系统归一化响应 t y (t ) 代码段2（求极限）： syms k limit('1-1/sqrt(1-k^2)*exp(-k*t)*sin(sqrt(1-k^2)*t+atan(sqrt(1-k^2)/k))',k,1) 运行结果： ans = (exp(t)-t-1)/exp(t) 7.题目： 构建一个简单的全波整流模型，并用示波器分两路同时观察原信号和整流后的信号波形。要求：信源模块被重新命名为“输入正弦波”；信宿模块被重新命名为“示波器”；连接到信宿的信号线上分别标注“原信号”和“整流信号”；在模型窗中添加注释文本。 方法1：用sine wave 产生正弦波，用abs 进行整流，用Scope 进行信号的观察。 示波器观测结果： 方法2： 示波器观测结果： 结果分析：两种方法都较好的实现了全波整流。 8.题目： 利用 SIMULINK 及其标准模块设计一个低通滤波器，从受噪声干扰的多频率混合信号)(t x 中获取10Hz 的信号。在此)()1002cos()102sin()(t n t t t x +??+??=ππ，而)2.0,0(~)(N t n 实现方法: 用Sine Wave 产生正弦波，用Random Number 产生随机信号，用sum 将信号相加，用Analog Filter Design 设计低通滤波器，用Mux 模块将信号合并在一起,用Scope 观察信号； 实验模型： 示波器观测结果： 结果分析：LPC 使用Butterworth 方法，实现了不错的低通滤波效果。 9.题目：已知系统的状态方程为：???=--=1222 211)1(y y y y y y ，其中25.0)0(,25.0)0(21==y y ，请构建该 系统的仿真模型，并用XY Graph 模块观察21,y y 相轨迹。 实现方法：利用乘法器，加法器，积分器搭建模型，利用XY Graph 模块观察21,y y 相轨。 实验模型： 运行结果： 10.题目：CD74HC00仿真（数字电路） 数字电路实验模型： 封装后功能检测模型： 运行结果： 结果分析：前两路分别为产生“1010”和“1100”的脉冲发生器的波形，后四路波形则分别为四个与非门的输出；系统对于“11”、“01”、“10”、“00”四种输入的输出符合与非逻辑关系，设计无误，可以实现预期功能。 11.题目： 2FSK信号的生成与接收滤波 要求： 1)产生二进制[0，1]随机信号波，信号波特率为1000B 2)对此信号进行2FSK调制，f1=1500Hz, f2=3050Hz 3)加入高斯白噪声 4)在接收端进行带通滤波 5)抽样速率为20KHz. 观察滤波前后信号的波形和频谱。 实现方法： 利用Bernoulli Binary Generator模块产生随机信号作为输入信号； 利用Sine Wave模块产生f1=1500Hz, f2=3050Hz 的载波，结合Switch模块进行2FSK调制； 利用AWGN Channel模块构建高斯白噪声信道，其参数设置为：信噪比：10dB，输入信号功率：1W； 利用Digital Filter Design设计带通滤波器，分别为1500Hz和3050Hz。 实验模型： 频谱观测结果： 示波器观测结果如下： 完整2FSK 信号生成与接收滤波实验模型： 为便于观察，设置信号波特率为50B ，示波器观测结果： 结果分析：略有误差与延迟，但大体上成功实现信号的滤波接收。 12.题目：创建一个简单的离散多速率系统：单位阶跃信号经过具有不同速率的采样后分 别用作两个离散传递函数的输入。这两个离散传递函数有相同的有理分式2.01 .0-+z z ，但采 样时间和时间偏置二元对分别设为[1 , 0.1]和[0.7 , 0]。要求：观察这两个离散传递函数的输出有什么不同；用不同的颜色标帜不同采样速率系统。（通信系统） 实验模型： 示波器观测结果： 二、M文件实现部分 题目：模拟调制---DSB调制与解调 要求：对正弦信号（300Hz）进行DSB调制和解调。信道采用高斯白噪声信道。 用MATLAB编写程序计算并画出调制信号波形和频谱。 程序代码： t=0:1/(400*1000-1):0.01; s=sin(300*2*pi*t); %原始信号 f=300Hz正弦信号 figure(1); subplot(3,1,1); plot(t,s); xlabel('t'); ylabel('amplitude'); %坐标轴标注 grid; %网格 title('原始信号'); %画出原始信号 z=ammod(s,1*10^5,4*10^5); %使用ammod函数进行调制 z = ammod(x,Fc,Fs) subplot(3,1,2); plot(t,z); xlabel('t'); ylabel('SDSB1(t)'); title('DSB-SC调制信号信号图'); %画出调制信号 y=awgn(z,5); %加入高斯白噪声信噪比SNR=5 subplot(3,1,3); plot(t,y); xlabel('t'); ylabel('SDSB2(t)'); title('加入噪声后的DSB-SC调制信号图，SNR=5'); %画出加入噪声后的DSB-SC图%%%%%%%%%%%%%%%%%%%%%%%以上是调制部分%%%%%%%%%%%%%%%%%%%%%%%%%% figure(2); f=(0:length(fft(z))-1)*(400*1000-1)/length(fft(z)); %FFT 离散时间傅里叶变换subplot(3,1,1); plot(f,abs(fft(z))); xlabel('f'); ylabel('M(f)'); grid; title('调制后的频谱图'); %以下为带通滤波器参数 As=100; Ap=1; fs=4*10^5; %采样频率 fp1=200; fp2=1.1*10^5; %通频带 fs1=100; fs2=1.2*10^5; %截止带 wp=[fp1 fp2]/(fs/2); ws=[fs1 fs2]/(fs/2); [n,wn]=cheb1ord(wp,ws,Ap,As); %带通滤波 b=fir1(n,wn,chebwin(n+1,As)); x=fftfilt(b,y); z1=amdemod(x,1*10^5,4*10^5);%滤波法解调 subplot(3,1,2); plot(t,z1); xlabel('t'); ylabel('S(t)'); title('解调波波形') f=(0:length(fft(z1))-1)*(400*1000-1)/length(fft(z1)); %FFT subplot(3,1,3); plot(f,abs(fft(z1))); axis([0 1000 0 2000]); %坐标轴长度设置 xlabel('f'); ylabel('M(f)'); grid; title('解调后的频谱图'); 运行结果： 0.002 0.004 0.006 0.008 0.01 -10 1t a m p l i t u d e 原始信号 0.002 0.0040.0060.0080.01 -10 1t S D S B 1(t ) DSB-SC 调制信号信号图 0.002 0.004 0.006 0.008 0.01 -50 5t S D S B 2(t ) 加入噪声后的DSB-SC 调制信号图，SNR=5 0.5 1 1.5 2 2.5 3 3.5 4 x 10 5 0500 1000 f M (f ) 调制后的频谱图 0.002 0.004 0.006 0.008 0.01 -20 2 t S (t ) 解调波波形 200 400 600 800 1000 01000 2000 f M (f ) 解调后的频谱图 南京邮电大学本科生毕业设计（论文）模板 南京邮电大学 毕业设计（论文）题目 专业 学生姓名 班级学号 指导教师 指导单位 日期：年月日至年月日 毕业设计（论文）原创性声明 本人郑重声明：所提交的毕业设计（论文），是本人在导师指导下，独立进行研究工作所取得的成果。除文中已注明引用的内容外，本毕业设计（论文）不包含任何其他个人或集体已经发表或撰写过的作品成果。对本研究做出过重要贡献的个人和集体，均已在文中以明确方式标明并表示了谢意。 论文作者签名： 日期：年月日 摘□□要 （“摘要”之间空两格，采用三号字、黑体、居中，与内容空一行）□□×××××××××（内容采用小四号宋体） 采用小四号、宋体、接排 小四号、黑体、顶格 ABSTRACT （采用三号字、Times New Roman字体、加黑、居中、与内容空一行）□□×××××××××（内容采用小四号Times New Roman字体） ：×××××；×××××；×××××；×××××； ××××× 采用小四号、Times New Roman字体、接排小四号、Times New Roman、加黑、顶格 目录 （三号、黑体、居中、目录两字空两格、与正文空一行） 第一章（空两格）☆☆☆（四号、宋体）……………………………×1.1☆☆☆☆（小四号宋体）………………………………………………………×1.2☆☆☆☆………………………………………………………………………×1.3☆☆☆☆………………………………………………………………………×……………… 第四章（空两格）☆☆☆（四号、宋体）……………………………×4.1☆☆☆☆………………………………………………………………………×4.2☆☆☆☆………………………………………………………………………×……………… 结束语（四号、宋体）………………………………………………×致谢（四号、宋体）……………………………………………………×参考文献（四号、宋体）………………………………………………×附录（四号、宋体）……………………………………………………× 不标页码 南京邮电大学 毕业设计(论文)任务书题目：OFDMA系统信道估计的研究 专业通信工程 学生姓名施祥威 班级学号10001417 指导教师(签字) 指导教师职称副教授 指导单位通信与信息学院 负责人签字 院(系)领导签字 日期2013.12 题目OFDMA系统信道估计的研究 课题类型理论研究■工程设计□ 产品开发□实验研究□ 课题类别 软件■硬件□ 软硬结合□其他□ 设计 内容与技术要求、成果形式 研究和学习无线移动通信原理，探讨数字通信系统的信号解调接收，针对信道对信号造成的失真，要进行信道估计，然后进行均衡（补偿）。具体到OFDMA系统，要力求准确估算多频率、时间点信道特性。 要求探讨无线信道的特点和模型，参考WiMAX系统标准，建立OFDMA系统模型。分析上、下行链路信道特点，分析不同导频符号安排下的信道估计算法。分析多种不同的信道估计算法，并作比较。在解调中用信道估计值进行补偿。 用Matlab仿真，绘出 ~,~ b MSE SNR BER E n曲线，分析其性能。提交程序清单，仿真数据。 设计进度11月中旬～12月初：查阅资料，学习无线通信和OFDMA的基础知识，消化课题内容，掌握仿真工具的使用方法，拟定系统总体实现方案。写开题报告。3月1日~3月底：依照总体方案，逐个编写程序，调试。联接程序，总体调试。 4月1日~5月底：调试程序，修改完善，测试系统数据、曲线，分析性能。撰写毕业论文。 6月1日~6月20日：完成毕业论文，毕业答辩。 参考资料[1]IEEE 802.16 IMT-Advanced Evaluation Group Coordination Meeting, Overview of IEEE P802.16m Technoledge and Candidate Rit for IMT-Advanced[S], IEEE L802.16-10/0002, IEEE,2010 [2](芬) 霍尔马 (Holma.H.) . (芬) 托斯卡拉 (Toskala.A.), UMTS中的LTE-基于OFDMA和SC-FDMA的无线接入, [M]机械工业出版社，2010 [3]Ramjee Prasad l Fernando J. Velez，WiMAX Networks: Techno-Economic Vision and Challenges[M], Springer press, 2010 [4]W-C.Huang,C-H Pan, C-P Li,H-J Li, Subspace-Based Semi-Blind Channel Estimation in Uplink OFDMA Systems[J], IEEE Trans. on Broadcasting, vol. 56, no. 1, March 2010 [5]沈嘉，下一代无线宽带通信的核心OFDM/OFDMA[J]，计算机世界，2007.1 [6]Yong Soo Cho , Jaekwon Kim, Won Young Yang , Chung G.Kang, MIMO-OFDM无线通信技术及MATLAB实现[M]，电子工业出版社，2013.6 [7]徐明远、邵玉斌，MATLAB仿真在通信与电子工程中的应用(第2版)[M]，西安电子科技大学出版社，2010 通信与信息工程学院 2015 / 2016 学年第 2 学期 软件设计实验报告 模块名称MATLAB软件设计 专业通信工程 学生班级B130102 学生学号B13010229 学生姓名郭智超 指导教师孟田甜 报告内容 一、要求练习的实验部分 1.题目：在时间区间 [0,10]中，绘出 t e y t 2cos 15.0--=曲线。 程序代码： t=0:0.02:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'b'); xlabel('time'); ylabel('amplitude'); title('figure 1'); grid 实验结果： 24 6810 00.51 1.5 time a m p l i t u d e figure 1 2.题目：写出生成如图E2-1所示波形的MATLAB 脚本M 文件。图中虚线为正弦波，要求它的负半波被置零，且在 2 3 处被削顶。 程序代码： t=linspace(0,10,500); %产生线性间隔的向量（在0-10之间取500点） y=sin(t); a=sin(pi/3); z=(y>=0).*y; z=(y>=a).*a+(y 课 程 设 计学院：数学学院 学号：20106496 姓名：黄星奕 辅导老师：陈晓红殷明 实验一 1．1 水手、猴子和椰子问题 一、问题描述 1．1 水手、猴子和椰子问题：五个水手带了一只猴子来到南太平洋的一个荒岛上，发现那里有一大堆椰子。由于旅途的颠簸，大家都很疲惫，很快就入睡了。第一个水手醒来后，把椰子平分成五堆，将多余的一只给了猴子，他私藏了一堆后便又去睡了。第二、第三、第四、第五个水手也陆续起来，和第一个水手一样，把椰子分成五堆，恰多一只猴子，私藏一堆，再去入睡，天亮以后，大家把余下的椰子重新等分成五堆，每人分一堆，正好余一只再给猴子，试问原先共有几只椰子？ 二、思考与实验 试分析椰子数目的变化规律，利用逆向递推的方法求解这一问题。 三、问题分析 用递推算法。首先分析椰子数目的变化规律，设最初的椰子数为p 0，即第一个水手所处理之前的椰子数，用p 1、p 2、p 3、p4、p 5分别表示五个水手对椰子动了手脚以后剩余的椰子数目，则根据问题有 再用x表示最后每个水手平分得到的椰子数，于是有 所以 p5 = 5x +1 利用逆向递推的方法，有 但由于椰子数为一正整数，用任意的x作为初值递推出的p0数据不一定是合适的。在实验中可以用for 循环语句结合break语句来寻找合适的x和p0，对任意的x递推计算出p0，当计算结果为正整数时，结果正确，否则选取另外的x再次重新递推计算，直到计算出的结果p0为正整数为止。 四、源程序 n=input('input n:'); for x=1:n p=5*x+1; for k=1:5 p=5*p/4+1; end if p==fix(p) break; end end disp([x,p]); 五、实验结果 六、结果分析 从理论上分析，由于 所以 “”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> ('(x)-3*x^2',0) = -2*(-1/6*3^(1/2)) -2*(-11/6*3^(1/2)) -2*(1/6*3^(1/2)) 3、求解下列各题： 1）30 sin lim x x x x ->- >> x; >> (((x))^3) = 1/6 2） (10)cos ,x y e x y =求 >> x; >> ((x)*(x),10) = (-32)*(x)*(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> x; >> ((((x^2),0,1/2)),17) = 0.54498710418362222 4）4 2 254x dx x +? >> x; >> (x^4/(25^2)) = 125*(5) - 25*x + x^3/3 5）求由参数方程arctan x y t ??=? =??dy dx 与二阶导 数22 d y dx 。 >> t; >> ((1^2))(t); >> ()() = 1 6）设函数(x)由方程e所确定，求y′(x)。>> x y; *(y)(1); >> ()() = (x + (y)) 7） sin2 x e xdx +∞- ? >> x; >> ()*(2*x); >> (y,0) = 2/5 8） 08x =展开（最高次幂为） >> x (1); taylor(f,0,9) = - (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 1 9） 1sin (3)(2)x y e y =求 >> x y; >> ((1)); >> ((y,3),2) = 邮电大学通信学院 软件课程设计 实验报告 模块名称：___MATLAB软件设计 专业班级：__通信工程 姓名：____ ____ 学号：___ 实验日期：2013年 6 月 17—28日实验报告日期： 2013 年 7 月 1 日 一、要求练习的实验部分 1．在时间区间 [0,10]中，绘出t e y t 2cos 15.0--=曲线。 程序： t=0:0.1:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); shg 结果： 2. 写出生成如图E2-1所示波形的MATLAB 脚本M 文件。图中虚线为正弦波，要求它的负半 波被置零，且在2 3处被削顶。 程序： t=linspace(0,3*pi,500); y=sin(t); a=sin(pi/3); z=(y>=0).*y; z=(y>=a).*a+(y xlabel('t'),ylabel('z=f(t)'),title('逐段解析函数') legend('y=sin(t)','z=f(t)',4),hold off 结果： 3. 令???? ??????=987654321A ，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题： （1） sum(diag(S)>tol) - rank(A) 是多少 ？ （2） S(1,1) - norm(A) = 0 是多少 ? （3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么 ? （4） S(1,1)/S(3,3) - cond(A) 的结果是什么 ? （5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？ （6） V(:,1)'*null(A) 得到什么结果 ? （7） abs(A*null(A)) < tol 得到什么结果 ? （8） U(:,1:2) = = orth(A) 的运行结果是什么 ? 程序： clear; clc; disp('设 A=[1,2,3;4,5,6;7,8,9]，得') 数值代数实验 数值线性代数实验一 一、实验名称：矩阵的LU分解. 二、实验目的：用不选主元的LU分解和列主元LU分解求解线性方程组Ax=b, 并比较这 两种方法. 三、实验内容与要求 （1）用所熟悉的计算机语言将不选主元和列主元LU分解编成通用的子程序，然后用编写的程序求解下面的84阶方程组 将计算结果与方程组的精确解进行比较，并就此谈谈你对Gauss消去法的看法. （2）写出追赶法求解三对角方程组的过程，并编写程序求该实验中的方程组 Gauss消去法： 用消去法解方程组的基本思想是用逐次消去未知数的方法把原来方程组Ax=b化为与其等价的三角方程组，而求解三角方程组就容易了。换句话说，上述过程就是用行的初等变换将原方程组系数矩阵化为简单形式，从而将求解原方程组的问题转化为求解简单方程组的问题。 利用Gauss消去法对线性方程组Ax=b进行求解。 用MATLAB建立m文件DelGauss.m,程序如下： function x=DelGauss(a,b) [n,m]=size(a); nb=length(b); det=1; x=zeros(n,1); for k=1:n-1 for i=k+1:n if a(k,k)==0 return end m=a(i,k)/a(k,k); for j=k+1:n a(i,j)=a(i,j)-m*a(k,j); end b(i)=b(i)-m*b(k); end det=det*a(k,k); end det=det*a(n,n); for k=n:-1:1 for j=k+1:n b(k)=b(k)-a(k,j)*x(j); end x(k)=b(k)/a(k,k); End 在matlab中输入如下： 结果如下： 南京邮电大学本科生毕业（论文） 南京邮电大学 毕业论文 题目 专业 学业姓名 班级学号 指导教师 指导单位 日期：年月日 摘要 本论文主要内容是基于本高校的论文格式要求，利用论文排版机器人PaperYes 的人工智能引擎排版出来的论文模板。各位同学可以通过此模板撰写论文，只需将自己撰写的论文内容填充到相应的标题下面就可以了，帮助同学们解决论文排版的困扰。如果有的同学不会用此模板编辑论文，可以选择PaperYes论文排版，根据您所在高校的格式要求,自动生成目录和页眉页脚,解决奇偶页难题，智能编排标题、图、表和公式,适配字体字号和段落等，10秒搞定，安全高效。 关键词:论文排版；机器人；PaperYes；论文模板 ABSTRACT The main content of this paper is based on the requirement of paper format in our university, using the artificial intelligence engine of PaperYes, a paper typesetting robot, to compose the paper template. Students can use this template to write papers, just fill the content of the papers under the corresponding headings, to help students solve the problem of paper typesetting. If some students do not use this template to edit papers, you can choose PaperYes paper typesetting. According to the format requirements of your university, automatically generate directories and headers and footers, solve the odd and even page problems, intelligently arrange titles, graphs, tables and formulas, and adapt font size and paragraph, etc., 10 seconds to complete, safe and efficient. Key words：PaperYes;Layout Technology;AIrtificial Intelligence; Research; Based on; “MATLAB”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6) 1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408 3）2sin cos 0x x x -= 所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans = 毕业设计（论文）模板 摘要 (“摘要”之间空两格，采用三号字、黑体、居中，与内容空一行) （内容采用小四号宋体）××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××× 关键词：(小四号、黑体、顶格)××；×××；×××；××× （内容采用小四号、宋体、接排、各关键词之间有1个分号及空格） ABSTRACT (采用三号字、Times New Roman字体、加黑居中、与内容空一行) （内容采用小四号Times New Roman字体）××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××× Keywords：(小四号、Times New Roman、黑体、顶格)××；×××；×××；××（内容采用小四号、Times New Roman字体、接排、各关键词之间有1个分号及空格） 目录 （三号、黑体、居中、目录两字空两格、空一行） 第一章题目（四号、宋体） (1) 1.1 题目（小四号、宋体）……………………………………………………………× 1.1.1 题目…………………………………………………………………………× 1.1.2 题目…………………………………………………………………………× 1.1.3 题目…………………………………………………………………………× 1.2 题目…………………………………………………………………………………× 1.3 题目…………………………………………………………………………………× 1.4 题目…………………………………………………………………………………× 1.5 题目…………………………………………………………………………………×第二章题目………………………………………………………………………………×第三章题目………………………………………………………………………………×结束语……………………………………………………………………………………×致谢……………………………………………………………………………………×参考文献……………………………………………………………………………………×附录X……………………………………………………………………………………× MATLAB 程序设计实验报告 一、实验目的 1. 通过实验熟悉MATLAB 仿真软件的使用方法； 2. 掌握用MATLAB 对连续信号时域分析、频域分析和s 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征； 3. 掌握用MATLAB 对离散信号时域分析、频域分析和z 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征； 4. 通过绘制信号运算结果的波形，了解这些信号运算对信号所起的作用。 二、实验设备 1. 计算机 ： 2. MATLAB R2007a 仿真软件 三、实验原理 对系统的时域分析 信号的时域运算包括信号的相加、相乘，信号的时域变换包括信号的平移、反折、倒相及信号的尺度变换。 （1）信号的相加和相乘：已知信号)(1t f 和)(2t f ，信号相加和相乘记为 )()(1t f t f =)(2t f +；)()(1 t f t f =)(2t f *。 （2）信号的微分和积分：对于连续时间信号，其微分运算是用diff 函数来完成的，其语句格式为：diff(function,’variable’,n)，其中function 表示需要进行求导运算的信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；n 为求导的阶数，默认值为求一阶导数。连续信号的积分运算用int 函数来完成，语句格式为：diff(function,’variable’,a,b)，其中function 表示需要进行被积信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；a,b 为积分上、下限，a 和b 省略时为求不定积分。 （3）信号的平移、翻转和尺度变换 信号的平移包含信号的左移与右移，信号的翻转包含信号的倒相与折叠，平移和翻转信号不会改变信号)(t f 的面积和能量。信号的尺度变换是对信号)(t f 在时间轴上的变化，可使信号压缩或扩展。)(at f 将原波形压缩a 倍，)/(a t f 将原波形扩大a 倍。 ￥ 对系统频率特性的分析 通信与信息工程学院 2014 /2015 学年第2学期 软件设计实验报告 模块名称MATLAB软件设计专业通信工程 学生班级通工2班 学生学号 学生姓名 指导教师孟田甜 报告内容 一、要求练习的实验部分 1． 题目：在时间区间 [0,10]中，绘出t e y t 2cos 15.0--=曲线。 实现方法：plot 画图 程序：clc clear t=0:0.01:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); xlabel('t'); ylabel('y'); title('绘制曲线'); 结果及结 果分析： 1234 5 678910 00.51 1.5 t y 绘制曲线 对源程序的改进：1.加了坐标轴标示，标题 2.将步长细化 2． 题目：写出生成如图E2-1所示波形的MATLAB 脚本M 文件。图中虚线为正弦 波，要求它的负半波被置零，且在 2 3 处被削顶。 图 E2-1 题15所要表示的波形（逐段解析函数） 实现方法：plot 画图 程序 t=0:0.01*pi:10; y=sin(t); a=sin(pi/3); z=(y>=0).*y; z=(y>=a).*a+(y 【关键字】模板 南京邮电大学本科生毕业设计（论文）模板 南京邮电大学 毕业论文 (自考助学本科二学历) 题目 专业 学生姓名 准考证号 指导教师 指导单位 日期：年月日至年月日 摘□□要 （“摘要”之间空两格，采用三号字、黑体、居中，与内容空一行）□□×××××××××（内容采用小四号宋体） 关键词：×××××；×××××；×××××；×××××； ××××× 采用小四号、宋体、接排 小四号、黑体、顶格 目录 （三号、黑体、居中、目录两字空两格、与正文空一行） 第一章（空两格）☆☆☆（四号、宋体）……………………………× 1.1☆☆☆☆（小四号宋体）………………………………………………………× 1.2☆☆☆☆………………………………………………………………………× 1.3☆☆☆☆………………………………………………………………………×……………… 第四章（空两格）☆☆☆（四号、宋体）……………………………× 4.1☆☆☆☆………………………………………………………………………× 4.2☆☆☆☆………………………………………………………………………×……………… 结束语（四号、宋体）………………………………………………× 致谢（四号、宋体）……………………………………………………× 参照文献（四号、宋体）………………………………………………× 附录（四号、宋体）……………………………………………………× 不标页码 第一章□□☆☆☆☆☆（居中、小三号、黑体） 1.1☆☆☆（四号、黑体、顶格） □□☆☆☆☆☆☆☆☆☆正文（小四号、宋体）宋体五号，居中，位于表上 表1.3□□☆☆☆ □□☆☆☆☆☆☆☆☆☆☆☆☆☆ （下一章另起一页） 第二章□□☆☆☆☆☆（居中、小三号、黑体） 2.1☆☆☆（四号、黑体、顶格） □□☆☆☆☆☆☆☆☆☆正文（小四号、宋体、空2格、用1.25倍行间距） 图2.5□□×××组织结构图 宋体五号居中，位于图下 图与下文空一行 正文开始标注页码;位置：页面底端（页脚）; 对齐方式：居中 结束语 （三号、黑体、居中、与正文空一行） （空两格，内容采用小四号宋体，行距 1.25倍）××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××× 摘 要 本文概述了信号仿真系统的需求、总体结构、基本功能。重点介绍了利用Matlab 软件设计实现信号仿真系统的基本原理及功能，以及利用Matlab 软件提供的图形用户界面（Graphical User Interfaces ,GUI ）设计具有人机交互、界面 友好的用户界面。本文采用 Matlab 的图形用户界面设计功能, 开发出了各个实 验界面。在该实验软件中, 集成了信号处理中的多个实验, 应用效果良好。本系 统是一种演示型软件,用可视化的仿真工具,以图形和动态仿真的方式演示部分 基本信号的传输波形和变换，使学习人员直观、感性地了解和掌握信号与系统的 基本知识。 随着当代计算机技术的不断发展，计算机逐渐融入了社会生活的方方面面。 计算机的使用已经成为当代大学生不可或缺的基本技能。 信号与系统课程具有传 统经典的基础内容，但也存在由于数字技术发展、计算技术渗入等的需求。在教 学过程中缺乏实际应用背景的理论学习是枯燥而艰难的。 为了解决理论与实际联系起来的难题国内外教育人士目光不约而同的投向 一款优秀的计算机软件——MATLAB 。 通过它可用计算机仿真， 阐述信号与系统 理论与应用相联系的内容，以此激发学习兴趣，变被动接受为主动探知，从而提 升学习效果，培养主动思维、学以致用的思维习惯。 以 MATLAB 为平台开发的信号与系统教学辅助软件可以充分利用其快速运 算，文字、动态图形、声音及交互式人机界面等特点来进行信号的分析及仿真。 运用 MATLAB 的数值分析及计算结果可视化、信号处理工具箱的强大功能将 信号与系统课程中较难掌握和理解的重点理论和方法通过概念浏览动态演示及典型例题分析等方式，形象生动的展现出来，从而使学生对所学知识理解更加透彻。同时运用教学软件中的数值计算工具箱，将学生从大量繁琐的手工数学运算中解放出来，将更多时间留给对基本概念和基本方法的思考。关键词：关键词：信号与系统，matlab，应用分析 ABSTRACT As the modern computer technology development, computer gradually merged with the various aspects of social life. the use of the computer essential part of the basic skills of students. the signal and systems of traditional classic of course, but there is also due to a digital technology development, the technology in the demand. in the teaching process of the practical application of theoretical study is boring and difficult. Theory with practice in order to solve the problem of educational circles and looked into a simultaneous of computer software —— matlab. it can be used by computer simulations, signals with the system theory and application related to the content, it aroused interest in learning and became passive acceptance of the initiative and thereby elevate learning, training of active in the make full use of its rapid operation of graphic, text, dynamic and interactive voice man-machine interface to the characteristics of the analysis and emulation. To the platform of the development of matlab signal system of teaching and assistive software can make full use of its rapid operation of graphic, text, dynamic and interactive voice man-machine interface to the characteristics of the 第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =，求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4). 和县移动分公司销售渠道的策略研究 摘要 自电信产业重组以来，全业务竞争带给整个行业内极大的改变，各家运营商在更新各自产品的同时，渠道代理商的抢夺更加激烈，本文将结合和县移动分公司自身的渠道发展情况，运用SOWT法对和县分公司的销售渠道进行分析，从而在培养核心竞争力、加大核心渠道保有以及优化渠道发展方面提出一定的策略。 关键词：渠道；SOWT；核心竞争力；优质渠道 ABSTRACT (采用三号字、Times New Roman字体、加黑居中、与内容空一行) （内容采用小四号Times New Roman字体）××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××××× Keywords：(小四号、Times New Roman、黑体、顶格)××；×××；×××；×× （内容采用小四号、Times New Roman字体、接排、各关键词之间有1个分号及空格） 目录 （三号、黑体、居中、目录两字空两格、空一行） 第一章题目（四号、宋体） (1) 1.1题目（小四号、宋体）……………………………………………………………× 1.1.1题目…………………………………………………………………………× 1.1.2题目…………………………………………………………………………× 1.1.3题目…………………………………………………………………………× 1.2题目…………………………………………………………………………………× 1.3题目…………………………………………………………………………………× 1.4题目…………………………………………………………………………………× 1.5题目…………………………………………………………………………………×第二章题目………………………………………………………………………………×第三章题目………………………………………………………………………………×结束语……………………………………………………………………………………×致谢……………………………………………………………………………………×参考文献……………………………………………………………………………………×附录X……………………………………………………………………………………- MATLAB入门教程 1．MATLAB的基本知识 1-1、基本运算与函数 在MATLAB下进行基本数学运算，只需将运算式直接打入提示号（>>）之後，并按入Enter键即可。例如： >> (5*2+1.3-0.8)*10/25 ans =4.2000 MATLAB会将运算结果直接存入一变数ans，代表MATLAB运算後的答案（Answer）并显示其数值於萤幕上。 小提示：">>"是MATLAB的提示符号（Prompt），但在PC中文视窗系统下，由於编码方式不同，此提示符号常会消失不见，但这并不会影响到MATLAB的运算结果。 我们也可将上述运算式的结果设定给另一个变数x： x = (5*2+1.3-0.8)*10^2/25 x = 42 此时MATLAB会直接显示x的值。由上例可知，MATLAB认识所有一般常用到的加（+）、减（-）、乘（*）、除（/）的数学运算符号，以及幂次运算（^）。小提示：MATLAB将所有变数均存成double的形式，所以不需经过变数宣告（Variable declaration）。MATLAB同时也会自动进行记忆体的使用和回收，而不必像C语言,必须由使用者一一指定.这些功能使的MATLAB易学易用，使用者可专心致力於撰写程式，而不必被软体枝节问题所干扰。 若不想让MATLAB每次都显示运算结果，只需在运算式最後加上分号（；）即可，如下例： y = sin(10)*exp(-0.3*4^2); 若要显示变数y的值，直接键入y即可： >>y y =-0.0045 在上例中，sin是正弦函数，exp是指数函数，这些都是MA TLAB常用到的数学函数。 下表即为MATLAB常用的基本数学函数及三角函数： 小整理：MATLAB常用的基本数学函数 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数z的相角(Phase angle) sqrt(x)：开平方 real(z)：复数z的实部 imag(z)：复数z的虚部 conj(z)：复数z的共轭复数 round(x)：四舍五入至最近整数 fix(x)：无论正负，舍去小数至最近整数 floor(x)：地板函数，即舍去正小数至最近整数 南京邮电大学通信学院 软件课程设计 实验报告 模块名称：___MATLAB软件设计 专业班级：__通信工程 姓名：____ ____ 学号：___ 实验日期：2013年 6 月 17—28日实验报告日期： 2013 年 7 月 1 日 一、要求练习的实验部分 1．在时间区间 [0,10]中，绘出t e y t 2cos 15.0--=曲线。 程序： t=0:0.1:10; y=1-exp((-0.5)*t).*cos(2*t); plot(t,y,'r-'); shg 结果： 2. 写出生成如图E2-1所示波形的MA TLAB 脚本M 文件。图中虚线为正弦波，要求它的 负半波被置零，且在2 3 处被削顶。 程序： t=linspace(0,3*pi,500); y=sin(t); a=sin(pi/3); z=(y>=0).*y; z=(y>=a).*a+(y xlabel('t'),ylabel('z=f(t)'),title('逐段解析函数') legend('y=sin(t)','z=f(t)',4),hold off 结果： 3. 令???? ??????=987654321A ，运行 [U,S,V]=svd(A);tol=S(1,1)*3*eps; 然后回答以下问题： （1） sum(diag(S)>tol) - rank(A) 是多少 ？ （2） S(1,1) - norm(A) = 0 是多少 ? （3） sqrt(sum(diag(S*S))) - norm(A,'fro') 的结果是什么 ? （4） S(1,1)/S(3,3) - cond(A) 的结果是什么 ? （5） S(1,1)*S(2,2)*S(3,3) - det(A) < tol 的运行结果是什么？ （6） V(:,1)'*null(A) 得到什么结果 ? （7） abs(A*null(A)) < tol 得到什么结果 ? （8） U(:,1:2) = = orth(A) 的运行结果是什么 ? 程序： clear; clc; disp('设 A=[1,2,3;4,5,6;7,8,9]，得')南京邮电大学毕设格式模板
南京邮电大学
南邮MATLAB软件设计
Matlab的实际应用设计(经典)
MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案
南京邮电大学matlab软件设计(超详细).
南京邮电大学 数值代数实验
南京邮电大学-本科-毕业论文-格式模板范文
MATLAB实验练习题(计算机) 南邮 MATLAB 数学实验大作业答案
南邮毕业设计(论文)模板
MATLAB程序设计实验报告
南邮matlab软件设计实验报告
【模板】南京邮电大学本生毕业设计论文模板
(完整版)matlab毕业设计论文
南邮MATLAB数学实验答案(全)
南邮毕业设计
matlab入门教程文献
南京邮电大学matlab软件设计(超详细).doc