关闭

关闭

封号提示

内容

首页 matlab上机实验答案.doc

matlab上机实验答案.doc

matlab上机实验答案.doc

上传者: 这个夏天有dian冷 2017-09-18 评分 0 0 0 0 0 0 暂无简介 简介 举报

简介:本文档为《matlab上机实验答案doc》,可适用于高等教育领域,主题内容包含matlab上机实验答案第一次上机内容:习题:A=cos(pi)(sqrt())^()习题:x=:*pi:*piy=cos(x)*(*sin(x)(符等。

matlab上机实验答案第一次上机内容:习题:A=cos(pi)(sqrt())^()习题:x=:*pi:*piy=cos(x)*(*sin(x)(x*x))plot(x,y)习题:x=::y=::X,Y=meshgrid(x,y)Z=X*X*exp(X^Y^)mesh(Z)习题:N=t=::*pix=cos(t)y=sin(N*t)figure()subplot(,,)plot(x,y)title('a=')y=sin(N*tpi)subplot(,,)plot(x,y)title('a=pi')y=sin(N*tpi)subplot(,,)plot(x,y)title('a=pi')y=sin(N*tpi)subplot(,,)plot(x,y)title('a=pi')选择一幅彩色图片,将其变作灰度图再做上下颠倒、左右颠倒、逆时针顺时针旋转度A=imread('jpg','jpg')A=double(A)A=A(:,:,)A(:,:,)A(:,:,)imshow(uint(A))h,w=size(A)B=Afori=:hB(i,:)=A(hi,:)endimshow(uint(B))C=Aforj=:wC(:,j)=A(:,wj)endimshow(uint(C))D=eye(w,h)forj=:wD(j,:)=A(:,wj)endimshow(uint(D))E=Aforj=:wfori=:hE(i,j)=A(hi,wj)endendimshow(uint(E))第二次上机内容鼠标点击空白图面实时显示鼠标点和轨迹figure()axis(,,,)x,y,b=ginput()plot(x,y,'r')string=sprintf('(f,f)',x,y)text(x,y,string)holdonwhile()ifb==breakendx=xy=yx,y,b=ginput()plot(x,y,'')string=sprintf('(f,f)',x,y)text(x,y,string)plot(xx,yy,'r')end正弦的轨迹动态表示axis(,*pi,,)x=y=sin(x)holdonx=xy=yforx=pi:pi:*pipause()y=sin(x)plot(xx,yy)x=xy=yend第三次上机内容:制作一界面上有二个输入框、一个计算按钮、一个显示框、一个下拉选择框做二数的可选择的四则运算。str=get(handlesedit,'string')number=strnum(str)str=get(handlesedit,'string')number=strnum(str)v=get(handlespopupmenu,'value')switchvcasem=numbernumbercasem=numbernumbercasem=number*numbercasem=numbernumberendset(handlesedit,'string',m)第四次上机内容:P第、、、题习题:求代数方程x^x^x^x^x=的根。p=,,,,,r=roots(p)习题:方程的根为x=,,,,求对应的x多项式的系数。x=,,,pp=poly(x)习题:a=,,,,,,,求其特征根和特征向量。A=,,,,,,lambda=eig(A)V,D=eig(A)习题:生成一个*的随机矩阵并对其进行三角分解和正交分解。A=rand(,)L,U=lu(A)Q,R=qr(A)编写子函数可对任意两个多项式进行加减操作(自动补零)functionc=myfun(x,y,p)i=length(x)j=length(y)ifi>jy=zeros(,ij),yelsex=zeros(,ji),xendif(p==)c=xyelsec=xyendreturn实现一副灰度图像的倍放大并在硬盘上保存该放大图像w=imread('F:pictures各类图片起司猫jpg')figure()imshow(uint(w))title('原彩图')b=rgbgray(w)figure()imshow(uint(b))title('灰度图')b=double(b)h,w=size(b)width=:wdepth=:hwi=:(w)(*w):wdi=:(h)(*h):hc=interp(width,depth',b,wi,di','cubic')figure()imshow(uint(c))title('放大图')c=uint(c)imwrite(c,'D:ProgramFilesmatlabbinwork起司猫放大图jpg','jpg')第五次上机内容:作业:已知描述某连续系统的微分方程为:y''(t)y'(t)y(t)=f(t)()用留数法求解脉冲响应和阶跃响应并画出响应的波形。b=a=,,t=::r,p=residue(b,a)y=len=length(r)fori=:leny=yr(i)*exp(p(i)*t)endsubplot(,,)plot(t,y)title('脉冲响应')a=,,,r,p=residue(b,a)y=len=length(r)fori=:leny=yr(i)*exp(p(i)*t)endsubplot(,,)plot(t,y)title('阶跃响应')()绘出该系统在,秒范围内并以时间间隔秒取样的冲激响应和阶跃响应的时域波形,以及频率响应波形b=a=,,t=::figure()subplot(,,)impulse(b,a,t)title('脉冲响应')subplot(,,)step(b,a,t)title('阶跃响应')w=::*pi*figure()freqs(b,a,w)title('频率响应')()求出系统在,秒范围内并以时间间隔秒取样的冲激响应和阶跃响应的数值解以及频率响应数值解。并用数值解画出波形。b=a=,,t=::y=impulse(b,a,t)figure()subplot(,,)plot(t,y)title('脉冲响应')y=step(b,a,t)subplot(,,)plot(t,y)title('阶跃响应')w=::*pi*H=freqs(b,a,w)mag=abs(H)phase=angle(H)*pifigure()subplot(,,)plot(w,mag)gridonxlabel('Frequency(rads)'),ylabel('Magnitude'),title('幅频曲线')subplot(,,)plot(w,phase)gridonxlabel('Frequency(rads)')ylabel('Phase(degrees)'),title('相频曲线')

职业精品

精彩专题

上传我的资料

热门资料

资料评价:

/ 9
所需积分:0 立即下载

意见
反馈

返回
顶部

Q