《杭州电子科技大学matlab课内实验(共6页).docx》由会员分享，可在线阅读，更多相关《杭州电子科技大学matlab课内实验(共6页).docx（6页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上第一次试验1. 画出公式y=cos(x)*(5=8*sin(x)+x*e x 分别画在0-2x内取1000个的曲线图，和x=【-3,3】取0.01间隔的曲线图x1=linspace(0,2*pi,100);y1=cos(x1).*(5+8*sin(x1)+x1.*exp(-x1);subplot(2,1,1);plot(x1,y1);x2=-3:0.01:3;y2=cos(x2).*(5+8*sin(x2)+x2.*exp(-x2);subplot(2,1,2);plot(x2,y2)2.画出公式x=cos（t+a），y=sin（2t），其中a=0，/6，/3，/2，
2、在4个子图中画出曲线t=-10:0.001:10;c1=0;x1=cos(t+c1);y1=sin(2*t);subplot(2,2,1);plot(x1,y1);c2=pi/6;x2=cos(t+c2);y2=sin(2*t);subplot(2,2,2);plot(x2,y2);c3=pi/3;x3=cos(t+c3);y3=sin(2*t);subplot(2,2,3);plot(x3,y3);c4=pi/2;x4=cos(t+c4);y4=sin(2*t);subplot(2,2,4);plot(x4,y4)3.读取一幅彩色图像，将其转变成灰度图像，将灰度图顺时针旋转90，将该图吧噢存
3、到硬盘上i=imread(C:Documents and SettingsAdministrator桌面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,:);endfigure;imshow(uint8(t);imwrite(uint8(t),C:Documents and SettingsAdministrator桌面dddd2.bmp);4.读取一幅图像，并显示图像，在图
4、像上用鼠标取两个点，并在图上用*号显示这两个点，以这两个点为左上角和右下角在图上画一个矩形框，截取该矩形框内的子图像，显示并保存在硬盘上i=imread(psfHomeDesktopddd2.bmp);figure;imshow(i);hold on;x,y=ginput(2);for m=1:2plot(x(m),y(m),*);endn=imcrop(i,min(x(1),x(2),min(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(WindowButtonDownFcn,md);axes(XLimMode,manual,yLimMode,manual);hold on;function md(src,eventdata)ud=get(src,userdata);cp=get(gca,currentpoint);x=cp(1,1);y=cp(1,2);plot(x,y,r*);if isempty(ud) line(ud(1),x,ud(2),y);endset(src,userdata,x,y);2.音乐柱状图的动态标识for i=1:1000 pause
6、(0.3); y=rand(1,10); bar(y);end3.编写一个函数，以图像文件路径为输入参数（1）显示灰度图像（2）图像灰度值通过y=x2变换并拉伸到【0,255】；（3）图像分为上下左右四个字块，并显示这四个子块图像；（4）把左上子块和右下子块图像进行交换，把右上和左下进行交换，并显示交换后的图像w=imread(C:Documents and SettingsAdministrator桌面图片.bmp);figure(1); b=rgb2gray(w); figure(1);imshow(uint8(b);title(灰度图); c=b.*b; figure(2);imshow
7、(uint8(c);title(改变灰度图);hang lie=size(b);figure(3);title(分割灰度图);subplot(2,2,1);imshow(uint8(b(1:hang/2,1:lie/2);subplot(2,2,2);imshow(uint8(b(1:hang/2,lie/2:lie);subplot(2,2,3);imshow(uint8(b(hang/2:hang,1:lie/2);subplot(2,2,4);imshow(uint8(b(hang/2:hang,lie/2:lie);figure(4); title(交换分割灰度图); subplot(
8、2,2,4); imshow(uint8(b(1:hang/2,1:lie/2); subplot(2,2,3);imshow(uint8(b(1:hang/2,lie/2:lie);subplot(2,2,2); imshow(uint8(b(hang/2:hang,1:lie/2);subplot(2,2,1);imshow(uint8(b(hang/2:hang,lie/2:lie);第三次试验制作一个GUI，实现简易计算器第四次试验1.求方程2x5+6x4+11x3+5x2+9x+12=0的跟y=2,6,11,5,9,12;x=roots(y);disp(x);2.设方程的跟为x=【-
9、5,3,8,9】，求对应的x多项式系数x=-5 ,3 ,8 ,9;y=poly(x);disp(y);3.编写子函数可对任意两个多项式进行加减操作（自动补零）unction y=jiafa(x1,x2)n1=length(x1);n2=length(x2);if n1n2 x2=zeros(1,n1-n2),x2;elseif n1n2 x2=zeros(1,n1-n2),x2;elseif n1n2 x1=zeros(1,n2-n1),x1;endy=x1-x2;4.实现一副灰度图像的2倍放大，并在硬盘上保存该放大图像w=imread(C:UsersAdministratorDesktop未
10、命名.bmp);b=rgb2gray(w);figure(1);imshow(uint8(b);title(灰度图);b=double(b);hang lie=size(b);width=1:lie;depth=1:hang;wi=1:(lie-1)/(2*lie-1):lie;di=1:(hang-1)/(2*hang-1):hang;c=interp2(width,depth,b,wi,di,cubic);figure(2);imshow(uint8(c);title(放大图);imwrite(c,D:图片未命名放大图.bmp,bmp);第五次试验已知描述某连续系统的微分方程为2y（t）+
11、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=length(r);for i=1:len; 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=length(r);for i=1:len; y=y+r(i)*exp(p(i)*t);endsubplot(2,1,2);plot(t,y);title(脉冲响应);2.绘出该系
12、统在0-30秒范围内，并以时间间隔0.01秒取样的冲击响应和阶跃响应的时域波形，以及频率响应波形a=2,1,8;b=1;t=0:0.01:30;figure(1);subplot(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;
13、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=angle(H)*180/pi;figure(2);subplot(2,1,1);plot(w,mag);grid on;xlabel(Frequency(rad/s),ylabel(Magnitud频曲线);subplot(2,1,2);plot(w,phase);grid on;xlabel(Frequency(rad/s),ylabel(Phase(deg)频曲线);专心-专注-专业