当前位置:文档之家› Matlab与通信仿真实验报告册答案

Matlab与通信仿真实验报告册答案

1.在Command Window里面计算(1)(358)510++÷⨯; >>(3+5+8)/5*10 ans=32(2)sin(3)π >>sin(3*pi)/sqrt(9/5) ans=2.7384e-16(3)123456789A⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,789456123B⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,计算:,,\,/C A BD A B A C C B=⨯=+;>> a=[1 2 3;4 5 6;7 8 9];>> b=[7 8 9;4 5 6;1 2 3];>> c=a*bc =18 24 3054 69 8490 114 138>> d=a+bd =8 10 128 10 128 10 12>> e=a\cWarning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 1.541976e-018.e =34.0000 22.0000 62.0000-50.0000 -23.0000 -100.000028.0000 16.0000 56.0000>> f=c/bWarning: Matrix is singular to working precision.f =NaN NaN NaNNaN -Inf InfNaN NaN NaN(4)3 1.247.5 6.6 3.15.4 3.46.1A⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦,求1,,A A A-';(求矩阵的行列式)>> a=[3 1.2 4;7.5 6.6 3.1;5.4 3.4 6.1]; >> a'ans =3.0000 7.5000 5.40001.2000 6.6000 3.40004.0000 3.1000 6.1000>> inv(c)c =2.1555 0.4555 -1.6449-2.1040 -0.2393 1.5013-0.7354 -0.2698 0.7833>> det(a)ans =13.7880(5)12345678i iZi i++⎡⎤=⎢⎥++⎣⎦,输入复数矩阵;>> z=[1+2*i 3+4*i;5+6*i 7+8*i] or>>z=[1 3;5 7]+[2 4;6 8]i z =1.0000 +2.0000i3.0000 +4.0000i5.0000 +6.0000i7.0000 +8.0000i2.建立.m文件,用for循环语句生成5×5的矩阵A:125236569⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦,将A矩阵进行水平和垂直翻转得到矩阵B和C。

将A矩阵的前3行,3列变成0并赋值给D。

>> a=zeros(5);for m=1:5for n=1:5a(m,n)=m+n-1endenda =1 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0a =1 2 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 01 2 3 4 52 3 4 5 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 7 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 7 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 7 8 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 7a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 7 85 6 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 7 85 6 7 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 7 85 6 7 8 0a =1 2 3 4 52 3 4 5 65 6 7 8 9c=flipud(a)c =5 6 7 8 9 4 5 6 7 8 3 4 5 6 7 2 3 4 5 6 1 2 3 4 5 d=fliplr(a)d =5 4 3 2 16 5 4 3 27 6 5 4 38 7 6 5 49 8 7 6 5 >> d=a;>> d(1:3,1:3)=0d =0 0 0 4 5 0 0 0 5 6 0 0 0 6 7 4 5 6 7 8建立的m文件a=zeros(5)for m=1:5for n=1:5a(m,n)=m+n-1endendd=a;d(1:3,1:3)=0建立的m文件,存盘时,不能存成1.m。

不运行,运行就为1.我存为mydesign.m默认工作目录是work文件夹edit pathdef 改路径最后一句加上cd d:\matlabBFile---Set path3.建立.m文件,随机产生一个5×5的矩阵,元素值为从0到255,要求用0和255对该矩阵进行标记,元素值大于等于128的标记为255,元素值小于128的标记为0。

a=255*rand(5);a =193.2237 180.0418 209.9817 111.8798 124.8899189.4988 8.1174 177.1813 97.2974 113.6245100.0179 70.6154 80.8604 195.2068 164.8098167.1469 11.7737 242.3066 202.7760 180.888043.6526 24.7686 8.7838 47.6525 192.4451a=255*rand(5);for i=1:5for j=1:5if a(i,j)>128a(i,j)=255;elsea(i,j)=0 end end enda =*********** 255 0 255 0 0 *********** 255 255 255 0 255 ***********4.编写程序使用0.618搜索法(即:近似黄金分割法)求给定函数的极值。

搜索法求解min ()t t ϕ≥的基本过程如下:给出[a,b],使得t 在[a,b]中,其中[a,b]称为搜索区间,通过迭代缩短[a,b]的长度。

当[a,b]的长度小于某个预设的值,或者导数的绝对值小于某个预设的正数,则迭代终止。

以函数3min ()21t t t t ϕ≥=-+作为处理对象,其中搜索区间定为[0,3],精度ε定为0.1。

图1-1 程序流程图b=3;er=0.5;while(1)m=a+0.382*(b-a); n=a+0.618*(b-a); p=m^3-2*m+1;q=n^3-2*n+1;if p<=qif n-a<ermbreak; elseb=n;endelseif b-m<ernbreak; elsea=m;endendn =0.8753实验二 绘图和确知信号分析实验1.假设N=12,对于M=4、5、7、10时,在0≤n≤2N -1区间使用plot 和stem 分别绘制信号)2sin(][NnM n x M π=,并添加适当标注。

clc;clear;M=[4,5,7,10]; N=12; n=0:2*N-1; for i=1:length(M) xn=sin(2*pi*M(i)*n/N); subplot(length(M),1,i) plot(n,xn); hold on ; stem(n,xn,'r'); s=int2str(M(i)) title(['M=' s]); endclc,clear; N=12; M=[4,5,7,10]; n=0:2*N-1;y1=sin(2*pi*M(1)*n/N); y2=sin(2*pi*M(2)*n/N); y3=sin(2*pi*M(3)*n/N); y4=sin(2*pi*M(4)*n/N); subplot(221)plot(n,y1);title('M=4');xlabel('n');ylabel('y1'); subplot(222)plot(n,y2);title('M=5');xlabel('n');ylabel('y2'); subplot(223)plot(n,y3);title('M=7');xlabel('n');ylabel('y3'); subplot(224)plot(n,y4);title('M=10');xlabel('n');ylabel('y4'); figure subplot(221)stem(n,y1);title('M=4');xlabel('n');ylabel('y1'); subplot(222)stem(n,y2);title('M=5');xlabel('n');ylabel('y2'); subplot(223)stem(n,y3);title('M=7');xlabel('n');ylabel('y3'); subplot(224)stem(n,y4);title('M=10');xlabel('n');ylabel('y4');2.设周期信号一个周期[0,T]的波形为()1,0/20/2t T s t T t T≤≤⎧=⎨<≤⎩,,其中T=1。

相关主题