MATLAB复习2011

第一次上机内容：1.P33-48的命令上机实践2.习题2、11、13、153.选择一幅彩色图片,将其变作灰度图，再做上下颠倒、左右颠倒、逆时针90，顺时针旋转180度a2=cos(pi/3)-(9-2.^(1/2)).^(1/3)%题2x=0:pi/124:2*pi;%题11y=cos(x)*(0.5+3*sin(x)/(1+x.^2));figure;plot(y);title('y=cos(x)*(0.5+3*sin(x)/(1+x.^2))' );xlabel('x') ;ylabel('y') ;x=-2:0.1:2; %题13y=x;[X,Y]=meshgrid(x,y);Z=(X.^2).*exp(-(X.^2+Y.^2));figure;surf(X,Y,Z)t=0:pi/50:4*pi;%题15N=2;x=cos(t);a=0;y=sin(N*t+a);subplot(2,2,1);plot3(x,y,t);title('a=0');grid on;a=pi/3;y=sin(N*t+a);subplot(2,2,2);plot3(x,y,t);title('a=pi/3');grid on;a=pi/2;y=sin(N*t+a);subplot(2,2,3);plot3(x,y,t);title('a=pi/2');grid on;a=pi;y=sin(N*t+a);subplot(2,2,4);plot3(x,y,t);title('a=pi');grid on; w=imread('C:\Documents and Settings\jin\My Documents\My Pictures\hello.jpg'); figure;subplot(3,3,1);imshow(uint8(w));title('原彩图');b=rgb2gray(w);subplot(3,3,2);imshow(uint8(b));title('灰度图');c=fliplr(b)subplot(3,3,3);imshow(uint8(c));title('灰度左右翻转图');d=flipud(b);subplot(3,3,4);imshow(uint8(d));title('灰度上下翻转图');[hang lie]=size(b);for i=1:hang %逆转90°for j=1:liem(lie+1-j,i)=b(i,j);% y(c+1-j,i)=y(i,j);endendsubplot(3,3,5);imshow(uint8(m));title('逆转90');for i=1:hang%顺转180°for j=1:lien(hang+1-i,lie+1-j)=b(i,j);%y(b+1-i,c+1-j)=y(i,j) endendsubplot(3,3,6);imshow(uint8(n));title('顺转180°');[hang lie]=size(b);for i=1:hang %顺转90°for j=1:liep(j,hang+1-i)=b(i,j);% y(c+1-j,i)=y(i,j);endendsubplot(3,3,7);imshow(uint8(p));title('顺转90');第二次上机内容：1.鼠标点击空白图面，实时显示鼠标点和轨迹。

2.正弦的轨迹动态表示axis([0 10 -10 10]);disp('单击鼠标左键点取气车下一个运动到的点')disp('单击鼠标右键确定气车运动的终点')[x1,y1,button]=ginput(1);while button==1plot(x1,y1,'g*');[x2,y2,button]=ginput(1);plot(x2,y2,'g*');plot([x1 x2],[y1 y2],'r');x1=x2;y1=y2;hold on;disp('单击鼠标左键点取气车下一个运动到的点')disp('单击鼠标右键确定气车运动的终点')endclose all;作业二axis([0 4*pi -2 2]);for x=0:0.1:4*piy=sin(x);hold on;plot(x,y,'r*');pause(0.01);End第三次上机内容：制作一界面，上有二个输入框、一个计算按钮、一个显示框、一个下拉选择框，做二数的可选择的四则运算。

%/***********GUI两数+ - * / ************/作业3str1=get(handles.edit1,'string');number1=str2num(str1);str2=get(handles.edit2,'string');number2=str2num(str2); v=get(handles.popupmenu1,'value');switch vcase 1m=number1+number2;case 2m=number1-number2;case 3m=number1*number2;case 4m=number1/number2;end;set(handles.edit3,'string',m)第四次上机内容：1.(旧版）p120，（新版）P124第4、5 、8、14题2.编写子函数可对任意两个多项式进行加减操作(自动补零)3.实现一副灰度图像的3倍放大，并在硬盘上保存该放大图像CHAP4function c=jiafa(x,y)i=length(x);j=length(y);if i>jy=[zeros(1,i-j) y];elsex=[zeros(1,j-i) x];endc=x+y;return c;w=imread('C:\Documents and Settings\jin\My Documents\图片\未命名.bmp');figure(1);imshow(uint8(w));title('原彩图');b=rgb2gray(w);figure(2);imshow(uint8(b));title('灰度图'); [hang lie]=size(b);width=1:lie;depth=1:hang;wi=1:(lie-1)/(3*lie-1):lie;di=1:(hang-1)/(3*hang-1):hang;c=interp2(width,depth',b,wi,di','cubic');figure(3);imshow(uint8(c));title('放大图');第五次上机内容：作业：已知描述某连续系统的微分方程为：2y' ' (t)+y' (t)+8y(t)=f(t)，试用MATLAB：（1）绘出该系统在0～30秒范围内，并以时间间隔0.01秒取样的冲激响应和阶跃响应的时域波形,以及频率响应波形；（2）求出系统在0～30秒范围内，并以时间间隔0.01秒取样的冲激响应和阶跃响应的数值解，以及频率响应数值解。

%注意运行后要等好几秒才会出结果num=[1];den=[2,1,8];t=0:0.01:30;figure(1);%冲激响应subplot(2,1,1);impulse(num,den,t);title('冲激响应');Y1=impulse(num,den,t)%数值解%阶跃响应subplot(2,1,2);step(num,den,t);title('阶跃响应');Y2=step(num,den,t)%数值解%频率响应a=[1];b=[2,1,8];w=0:0.1:2*pi*5;figure(2);freqs(a,b,w)[H,M]=freqs(a,b,w);H %数值解figure(3);subplot(2,1,1);plot(log(M),log(abs(H))),grid on;title('幅频曲线');subplot(2,1,2);c=atan2(imag(H),real(H))./pi.*180;semilogx(M,c);grid on;legend('ψ(H(w))');xlabel('w(rad/s)'),ylabel('ψ(H(w))(degrees)'),title('相频曲线'); 未做：Ex1编写一个M文件，完成如下功能：1. 读入一幅图像（图像位置在Matlab 安装目录下面的\toolbox\images\imdemos\cameraman.tif) ，并显示图像；2. 在步骤1的图像上用鼠标取两个点，并在图上用’*’号显示这两个点；3. 以步骤2的点为左上角和右下角在图上绘制一个矩形框；4. 截取步骤3矩形框内的子图像，并显示结果。

w=imread('C:\Documents andSettings\jin\My Documents\My Pictures\hello.jpg');figure(1);imshow(uint8(w));title('原彩图');hold on;[x,y,but]=ginput(2);plot(x,y,'*r');for i=y(1):y(2)plot(x(1),i,'r');plot(x(2),i,'r');endfor j=x(1):x(2)plot(j,y(1),'r');plot(j,y(2),'r');endfigure(2);imshow(uint8(w(y(1):y(2),x(1):x(2) ,:)));title('原彩图');%行，列，y是行,x是列Ex2编写一个函数，以图像文件路径为输入参数，1. 显示灰度图像；2. 图像灰度值通过y=x2变换并拉伸到[0 255]；3. 图像分成上下左右四个子块，并显示这四个子块图像；4. 把左上子块和右下子块图像进行交换，把右上子块和左下子块图像进行交换，并显示交换了子块后的图像。