杭州电子科技大学matlab课内实验#精选

1、第一次试验1.画出公式y=cos(x)*(5=8*sin(x)+x*e x分别画在0-2x内取1000个的曲线图, 和x=【-3,3】取0.01间隔的曲线图x1=li nspace(0,2*pi,100);y仁 cos(x1).*(5+8*si n(x1)+x1.*exp(-x1);subplot(2,1,1);plot(x1,y1);x2=-3:0.01:3;y2=cos(x2).*(5+8*si n( x2)+x2.*exp(-x2);subplot(2,1,2);plot(x2,y2)2. 画出公式 x=cos (t+a ), y=sin (2t),其中 a=0 , n /6 , n /

2、3 , n /2，在 4 个子图中 画出曲线t=-10:0.001:10;c1=0;x1=cos(t+c1);y1=si n( 2*t); subplot(2,2,1);plot(x1,y1);c2=pi/6; x2=cos(t+c2);y2=si n(2*t); subplot(2,2,2);plot(x2,y2);c3=pi/3; x3=cos(t+c3);y3=si n( 2*t); subplot(2,2,3);plot(x3,y3);c4=pi/2; x4=cos(t+c4);y4=si n(2*t); subplot(2,2,4);plot(x4,y4)3.读取一幅彩色图像，将其转

3、变成灰度图像，将灰度图顺时针旋转90 ,将该图吧噢存到硬盘上桌面 dddd.bmp');a=double(i);pr=a(:,:,1);pg=a(:,:,2);pb=a(:,:,3); a1=(pr+pg+pb)/3;x=size(a1);t=zeros(x(2),x(1);for i=1:x(1)t(:,x(1)-i+1)=a1(i,:)'end |figure;imshow( uin t8(t);桌面 dddd2.bmp');4.读取一幅图像，并显示图像，在图像上用鼠标取两个点，并在图上用*号显示这两个点，以这两个点为左上角和右下角在图上画一个矩形框，截取该矩形框内

4、的子图像，显示并保存在硬盘上i=imread('psfHomeDesktopddd2.bmp');figure;imshow(i);hold on;x,y=gi nput(2);for m=1:2plot(x(m),y(m),'*');endn=imcrop(i,mi n(x(1),x(2),mi n( y(1),y (2),abs(x(1)-x (2) ),abs(y(1)-y (2); figure;imshow( n);imwrite( n,'psfHomeDesktopddd3.bmp');第二次试验1. 鼠标点击空白图面，实时显示鼠标点

5、和轨迹function dsafigure('Wi ndowButto nDow nFcn',md); axes('XLimMode','ma nual','yLimMode','ma nual'); hold on;function md(src,eve ntdata) ud=get(src,'userdata');cp=get(gca,'curre ntpo in t'); x=cp(1,1);y=cp(1,2);plot(x,y,'r*');if isempty

6、(ud)lin e(ud(1),x,ud (2),y);endset(src,'userdata',x,y);2. 音乐柱状图的动态标识for i=1:1000pause(0.3);y=ra nd(1,10);bar(y);end(4 )把左上子块和右下子块图像进行交换，把右上和左下进行交换，并显示交换后的图像桌面图片.bmp');figure(1);b=rgb2gray(w);figure(1);imshow(uint8(b);title('灰度图');c=b.*b;figure(2);imshow(uint8(c);title('改变灰度图&

7、#39;);hang lie=size(b);figure(3);title('分割灰度图');subplot(2,2,1);imshow(ui nt8(b(1:ha ng/2,1:lie/2);subplot(2,2,2);imshow(ui nt8(b(1:ha ng/2,lie/2:lie);subplot(2,2,3);imshow(ui nt8(b(ha ng/2:ha ng,1:lie/2);subplot(2,2,4);imshow(ui nt8(b(ha ng/2:ha ng,lie/2:lie);figure(4);title('交换分割灰度图'

8、;);subplot(2,2,4);imshow(ui nt8(b(1:ha ng/2,1:lie/2);subplot(2,2,3);imshow(ui nt8(b(1:ha ng/2,lie/2:lie);subplot(2,2,2);imshow(ui nt8(b(ha ng/2:ha ng,1:lie/2);subplot(2,2,1);imshow(ui nt8(b(ha ng/2:ha ng,lie/2:lie);第三次试验制作一个GUI，实现简易计算器第四次试验1. 求方程 2x5+6x4+11x3+5x2+9x+12=0的跟y=2,6,11,5,9,12;x=roots(y);

9、disp(x);2. 设方程的跟为x=【-5,3,8,9】，求对应的x多项式系数x=-5 ,3 ,8 ,9;y=poly(x);disp(y);3. 编写子函数可对任意两个多项式进行加减操作(自动补零)un cti on y=jiafa(x1,x2)n1=le ngth(x1);n2=le ngth(x2);if n1> n2x2=zeros(1, n1- n2),x2;elseif n1< n2x1=zeros(1, n2-n1),x1;endy=x1+x2;function y=jia nfa(x1,x2)n1=le ngth(x1);n2=le ngth(x2);if n1&

10、gt; n2x2=zeros(1, n1- n2),x2;elseif n1< n2x1=zeros(1, n2-n1),x1;endy=x1-x2;4. 实现一副灰度图像的 2倍放大，并在硬盘上保存该放大图像未命名.bmp'); b=rgb2gray(w);figure(1);imshow( uin t8(b);title('灰度图');b=double(b);hang lie=size(b);width=1:lie;depth=1:ha ng;wi=1:(lie-1)/(2*lie-1):lie;di=1:(ha ng-1)/(2*ha ng_1):ha ng

11、;c=in terp2(width,depth',b,wi,di','cubic');figure(2);imshow(ui nt8(c);title('放大图');图片未命名放大图.bmp','bmp');第五次试验已知描述某连续系统的微分方程为2y ''( t) +y '(t) +8y(t) =f (t)1. 用留数法求解脉冲响应和阶跃响应，并画出响应的被波形a=2,1,8;b=1;t=0:0.01:30;r,p=residue(b,a);y=0;len=le ngth(r);for i=1:l

12、e n;y=y+r(i)*exp(p(i)*t);endsubplot(2,1,1);plot(t,y);title('脉冲响应');a=2,1,8,0;r,p=residue(b,a);y=0;len=le ngth(r);for i=1:le n;y=y+r(i)*exp(p(i)*t);endsubplot(2,1,2);plot(t,y);title('脉冲响应');2. 绘出该系统在0-30秒范围内，并以时间间隔 0.01秒取样的冲击响应和阶跃响应的时域 波形，以及频率响应波形a=2,1,8;b=1;t=0:0.01:30;figure(1);subp

13、lot(2,1,1);impulse(b,a,t); title('脉冲响应'); subplot(2,1,2); step(b,a,t);title('阶跃响应'); w=0:0.1:2*pi*5; figure(2); freqs(b,a,w);title('频率响应');3. 求出系统在0-30秒范围内，并以时间间隔 以及频率响应数值解，并用数值解画出波形0.01秒取样的冲激响应和阶跃响应的数值解,a=2,1,8;b=1;t=0:0.01:30;y1=impulse(b,a,t);figure(1);subplot(2,1,1);plot(t,y1);title('脉冲响应');y2=step(b,a,t);subplot(2,1,2);plot(t,y2);title('阶跃响应');w=0:0.1:2*pi*5;H=freqs(b,a,w);mag=abs(H);phase=a ngle(H)*180/pi;figure(2);subplot(2,1,1);plot(w,mag);gri