北科大MATLAB第五次作业
【VIP专享】MATLAB实验报告实验五
实验五1.M文件如下：函数fushu.M文件：function [e,l,s,c] = fushu(z)%fushu 复数的指数，对数，正弦，余弦的计算%e 复数的指数函数值%l 复数的对数函数值%s 复数的正弦函数值%c 复数的余弦函数值e=exp(z);l=log(z);s=sin(z);c=cos(z);命令文件M：z=input('请输入一个复数z=');[a,b,c,d]=fushu(z)请输入一个复数z=1+ia =1.4687 +2.2874ib =0.3466 + 0.7854ic =1.2985 + 0.6350id =0.8337 - 0.9889i2.函数fc.M文件:function X= fc(A,B)%fc fc是求解线性方程的函数%A A是未知矩阵的系数矩阵X=A\B;命令文件M：m1=input('请输入m1的值m1=');m2=input('请输入m2的值m2='); theta=input('请输入theta=');x=pi/180*theta;g=9.8;A=[m1*cos(x) -m1 -sin(x) 0m1*sin(x) 0 cos(x) 00 m2 -sin(x) 00 0 -cos(x) 1];B=[0;m1*g;0;m2*g];X=fc(A,B)请输入m1的值m1=10请输入m2的值m2=10请输入theta=60X =9.69952.424928.0000112.00003.函数文件：function [p]=prime(p)m=p(length(p));for i=2:sqrt(m)n=find(rem(p,i)==0&p~=i);p(n)=[];endp命令文件：p=10:99;p=prime(p);p=10*rem(p,10)+fix(p/10);p=prime(p)p =Columns 1 through 2011 13 17 19 23 29 31 37 41 43 47 53 59 61 67 71 73 79 83 89Column 2197p =11 31 71 13 73 17 37 97 794.函数文件：function f=fx(x)A=(x-2).^2+0.1;B=(x-3).^4+0.01;f=1./A+1./B;命令文件：x=input('请输入矩阵x：'); f=fx(x)请输入矩阵x：[1,2;3,4]f =0.9716 10.9901 100.9091 1.23405.(1)函数文件：function f=f(x)f=x+10.*log(x.^2+5);命令文件：y=f(40)./(f(30)+f(20))y =0.6390(2)函数文件：function g=g(x)sum=0;for i=1:xj=x.*(x+1);sum=sum+j;g=sum;end命令文件：y=g(40)./(g(30)+g(20)) y =1.8072。
matlab习题及答案
matlab习题及答案2. ⽤MATLAB 语句输⼊矩阵A 和B3．假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数⾏提取出来，赋给B 矩阵，⽤magic(8)A =命令⽣成A 矩阵，⽤上述命令检验⼀下结果是不是正确。
4．⽤数值⽅法可以求出∑=++++++==6363622284212i i S ，试不采⽤循环的形式求出和式的数值解。
由于数值⽅法是采⽤double 形式进⾏计算的，难以保证有效位数字，所以结果不⼀定精确。
试采⽤运算的⽅法求该和式的精确值。
5．选择合适的步距绘制出下⾯的图形。
（1）)/1sin(t ，其中)1,1(-∈t ；（2）)tan(sin )sin(tan t t -，其中),(ππ-∈t6. 试绘制出⼆元函数2222)1(1)1(1),(yx yx y x f z ++++-==的三维图和三视图7. 试求出如下极限。
（1）xxxx 1)93(lim +∞→；（2）11lim0-+→→xy xy y x ；（3）22)()cos(1lim222200y x y x ey x y x +→→++-8. 已知参数⽅程-==tt t y t x sin cos cos ln ，试求出x y d d 和3/22d d π=t x y9. 假设?-=xyt t e y x f 0d ),(2，试求222222y fy x f x f y x ??+-?? 10. 试求出下⾯的极限。
（1）-++-+-+-∞→1)2(1161141121lim 2222n n ；（2）)131211(lim 2222ππππn n n n n n n ++++++++∞→ 11. 试求出以下的曲线积分。
（1）?+ls y x d )(22，l 为曲线)sin (cos t t t a x +=，)cos (sin t t t a y -=，)20(π≤≤t 。
（2）?-+++ly y y xe x e yx )dy 2(xy d )(33，其中l 为22222c y b x a =+正向上半椭圆。
第5章MATLAB绘图_习题答案
第5章MATLAB绘图习题5一、选择题1．如果x、y均为4×3矩阵，则执行plot(x,y)命令后在图形窗口中绘制（）条曲线。
DA．12B．7C．4D．32．下列程序的运行结果是（）。
Ax=0:pi/100:2*pi;forn=1:2:10plot(n*sin(x),n*cos(x))holdonendaxissquareA．5个同心圆B．5根平行线C．一根正弦曲线和一根余弦曲线D．5根正弦曲线和5根余弦曲线3．命令text(1,1,'{\alpha}+{\beta}')执行后，得到的标注效果是（）。
C A．{\alpha}+{\beta}B．αβ}C．α+βD．αβ4．subplot(2,2,3)是指（）的子图。
AA．两行两列的左下图B．两行两列的右下图C．两行两列的左上图D．两行两列的右上图x的曲线绘制成直线，应采用的绘图函数是（）。
C5．要使函数y=2eA．polarB．semilogxC．semilogyD．loglog6．下列程序的运行结果是（）。
B[x,y]=meshgrid(1:5);surf(x,y,5*ones(size(x)));A．z=x+y平面B．与xy平面平行的平面C．与xy平面垂直的平面D．z=5x平面7．下列函数中不能用于隐函数绘图的是（）。
DA．ezmeshB．ezsurfC．ezplotD．plot38．下列程序运行后，看到的图形（）。
Ct=0:pi/20:2*pi;[x,y]=meshgrid(-8:0.5:8);z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);surf(x,y,z)view(0,90);axisequalA．像墨西哥帽子B．是空心的圆C．边界是正方形D．是实心的圆9．下列程序运行后得到的图形是（）。
A[x,y]=meshgrid(-2:2);z=x+y;i=find(abs(x)<1&abs(y)<1);z(i)=NaN;surf(x,y,z);shadinginterpA．在一个正方形的正中心挖掉了一个小的正方形B．在一个正方形的正中心挖掉了一个小的长方形C．在一个正方形的上端挖掉了一个小的正方形D．在一个正方形的下端挖掉了一个小的正方形10．在使用MATLAB“绘图”选项卡中的命令按钮绘图之前，需要（）。
matlab课后习题及答案详解
matlab课后习题及答案详解第1章练习题1.安装matlab时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？在安装matlab时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即matlab选项）必须安装。
第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。
2.matlab操作方式桌面存有几个窗口？如何并使某个窗口瓦解桌面沦为单一制窗口？又如何将瓦解过来的窗口再次置放至桌面上？与其他计算机语言相比较，matlab语言注重的特点就是什么？matlab系统由那些部分共同组成？在matlab操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的close按钮，一个是可以使窗口成为独立窗口的undock按钮，点击undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择dock……菜单项就可以将独立的窗口重新防止的桌面上。
matlab具备功能强大、使用方便、输出简便、库函数多样、开放性弱等特点。
matlab系统主要由开发环境、matlab数学函数库、matlab语言、图形功能和应用程序接口五个部分组成。
3.如何设置当前目录和搜寻路径，在当前目录上的文件和在搜寻路径上的文件存有什么区别？命令历史窗口除了可以观测前面键入的命令外，除了什么用途？当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的setpath菜单项来完成。
在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被matlab运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。
命令历史窗口除了用作查阅以前键入的命令外，还可以轻易执行命令历史窗口中选取的内容、将选取的内容拷贝到剪贴板中、将选取内容轻易拷贝到m文件中。
数字图像处理第五次作业
数字图像处理第五次作业摘要本次报告主要记录第五次作业中的各项任务完成情况。
本次作业以Matlab 2013为平台，结合matlab函数编程实现对所给图像文件的相关处理：1. 频域低通滤波器：设计低通滤波器包括 butterworth and Gaussian (选择合适的半径，计算功率谱比),平滑测试图像test1和2；2.频域高通滤波器：设计高通滤波器包括butterworth and Gaussian，在频域增强边缘。
选择半径和计算功率谱比，测试图像test3,4；3. 其他高通滤波器：拉普拉斯和Unmask，对测试图像test3,4滤波；4.比较并讨论空域低通高通滤波（Project4）与频域低通和高通的关系。
以上任务完成后均得到了预期的结果。
1. 频域低通滤波器：设计低通滤波器包括 butterworth and Gaussian (选择合适的半径，计算功率谱比),平滑测试图像test1和2。
（1）实验原理及方法低通滤波是要保留图像中的低频分量而除去高频分量。
图像中的边缘和噪声都对应图像傅里叶频谱中的高频部分，所以低通滤波可以除去或消弱噪声的影响并模糊边缘轮廓。
理想低通滤波器具有传递函数：其中为制定的非负数，D(u,v)为点(u,v)到滤波器中心的距离。
功率谱定义：为滤波前图像的功率谱,为滤波后图像的功率谱。
滤波器的功率谱理解为：1>Butterworthn阶Butterworth低通滤波器（BLPF）的传递函数（截止频率距原点的距离为）定义如下：其中，。
BLPF变换函数在通带与被滤除的频率之间没有明显的截断。
对于有平滑传递函数的滤波器，定义一个截止频率的位置并使H(u，v)幅度降到其最大值的一部分。
在上式中，当时，H（u，v）=0.5（从最大值降到它的50%）。
一阶的巴特沃斯滤波器没有振铃，在二阶中振铃通常很微小，这是因为与理想低通滤波器相比，它的通带与阻带之间没有明显的跳跃，在高低频率间的过渡比较光滑。
MATLAB教程课后实验报告题目及解答[第一至第五章]
![MATLAB教程课后实验报告题目及解答[第一至第五章]](https://img.taocdn.com/s3/m/512b4436f111f18583d05a96.png)
MATLAB教程实验报告实验项目名称实验一 Matlab基本操作学生姓名汪德旺专业班级 09数教（1）班学号 0301090131实验成绩日期一. 实验目的和要求1、了解MATLAB 的开发环境。
2、熟悉Matlab的基本操作。
3、掌握建立矩阵的方法。
4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。
5、填写实验报告，实验报告文件取名为report1.doc。
6、于邮件附件形式将实验报告文件report1.doc 发到邮箱*******************，邮件主题为班级学号姓名，如：09数教1班15号张三。
二、实验内容1、先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。
（1）0 122sin851ze =+（2）2212 1ln(0.4552i z x x+⎡⎤==⎢⎥-⎣⎦其中（3）0.30.330.3sin(0.3)ln,22a ae e az a--+=++a=-3.0,-2.9,-2.8,…, 2.8, 2.9,3.0（4）2242,011,12,0:0.5:2.521,23t tz t tt t t⎧≤<⎪=-≤<=⎨⎪-+≤<⎩其中t2.已知：1234413134787,2033657327 A B--⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：（1）A+6*B和A-B+I(其中I为单位矩阵) （2）A*B和A.*B（3）A^3和A.^3（4）A/B和B\A（5）[A,B]和[A([1,3],:);B^2]3、设有矩阵A 和B123453016678910A=,B=17-691112131415023-41617181920970212223242541311⎡⎤⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦（1）求它们的乘积C 。
（2）将矩阵C 的右下角3*2子矩阵赋给D 。
（3）查看MATLAB 工作空间的使用情况。
matlab实验报告(实验5)
一、 实验目的
开课学院及实验室：机电学院
2012 年 12 月 28 日
机电学
学院
院
实验课程名称
年级、专 业、班
1、 掌握 M 脚本文件和函数文件的创建方法，理解两者在形式上和运行上的区别； 2、 理解基本空间和函数空间的含义。 3、 掌握以控制流进行 MATLAB 程序设计的方法以及向量化程序设计方法，了解两者之间的异同； 二、
为频域抽样间隔， 2 1 为带宽。 K
四、实验报告要求
1、简述实验目的、实验原理及其要点； 2、把所有实验内容涉及的指令都以 M 文件方式创建并运行，保存实验结果图形。 3、总结实验中的主要结论。
按照实验内容 1 定义的 f(t)及时域采样方式，并取频域抽样点数 K=100，利用 MATLAB 编程，使用 for 循 环结构计算 F ( ) 在 [8 ,8 ] 在各个频域采样点的值，并画出其频谱图形。
T N 1 F (1 k ) f (t1 nt )e j (1 k )(t1 nt ) N n 0
其中，
（2.1）
根据以上两个等式，利用 MATLAB 编程，使用向量化计算的方法重新完成实验内容 2 和实验内容 3，并把计算结 果与实验内容 2 和实验内容 3 的结果进行比较是否相同（通过计算两者之间的最大绝对误差，观察其值是否接 近 0） 。 6、利用实验内容 5 的代码，把傅里叶变换和傅里叶逆变换分别写为两个 M 函数文件，函数定义形式为： function [F, w] = Fourier_lab2(f, N, T, K) function [f, t] = iFourier_lab2(F, K, W, N) 其中，f 为时域信号，t 为对时域信号 f 的采样点（即采样时刻） 为对时域信号 f 的采用点数，T 为 f 的 ，N 持续时间，K 为频域采样点数，F 为 f 的傅里叶变换，w 为对 F 的频域采样点，W 为 F 的频域带宽。 利用以上两个函数，写出 M 脚本文件重新完成实验内容 5。
5章习题答案matlab
张卫华 MATLAB课堂
例13：
分析以下程序， 并运行观察。
clf;x=3*pi*(-1:0.05:1);y=x; [X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps; Z=sin(R)./R; h=surf(X,Y,Z);colormap(jet); axis off n=12;mmm=moviein(n); for i=1:n rotate(h,[0 0 1],25); mmm(:,i)=getframe; end movie(mmm,5,10)
张卫华 MATLAB课堂
例5程序：
fplot('cos(tan(pi*x))',[ 0,1],1e-4)
张卫华 MATLAB课堂
例6：
绘制r=sin(t)cos(t)的极坐标图，并标记 数据点。
张卫华 MATLAB课堂例6 Nhomakorabea序：t=0:pi/50:2*pi; r=sin(t).*cos(t); polar(t,r,'-*');
title('y1=0.2e^{-0.5x}cos(4\pix) 和y2=2e^{0.5x}cos(4\pix)比较') text(2,2*exp(-1),'\fontsize{20}\bf \leftarrow x_{1}=2,y_{1}=2e^{-1}')
张卫华 MATLAB课堂
例5：
绘制f(x)=cos(tan(πx))的曲线
张卫华 MATLAB课堂
例11：
绘制柱形，剪切掉x、y小于零 的部分
张卫华 MATLAB课堂
例11程序：
t=linspace(0,2*pi,100); r=1-exp(-t/2).*cos(4*t); [X,Y,Z]=cylinder(r,60); ii=find(X<0&Y<0); Z(ii)=NaN; surf(X,Y,Z);colormap(spring), shading interp light('position',[-3,-1,3],'style','local')
MATLAB基础习题第五章习题答案
disp('平手');
end
switch(you-compute)
case{1,-2}
disp('您输了');
case{-1,2}
disp('您赢了');
end
%%
%是否继续判别模块
n=input('是否继续玩该游戏？否（0），是（1）');
while n~=0&n~=1
imwrite(I_gray,'1_gray.bmp', 'bmp');
subplot(1,3,2);
imshow(I_gray);
title('灰度图')
%%变换后的灰度图
I1=I_gray;
I1(I_gray>200)=200;
I1(I_gray<100)=100;
subplot(1,3,3);
该模块实现石头剪刀布的功能石头剪刀布由用户输入电脑实际出然后判断输赢然后根据用户判断是否继续根据用户判断进行下一步的操作clc
第五章习题
1.读取一副图像，如果是彩色图则转化为灰度图；对得到的灰度图做如下处理：
（1）将灰度值大于200的像素值全部设置为200；
（2）将灰度值小于100的像素值全部设置为100；
close all;
%%
%方程输入模块
x=inputdlg({'第一个方程系数','第二个','第三个','第四个','第五个','第六个'});
A=cell2mat(x);
实验五 matlab基础知识(简单)
本次实验注意：《实验五MALTAB基础知识（简单）》《实验五基于Matlab的信号频谱分析（复杂）》选作一个即可实验五MALTAB基础知识（一）实验目的 (2)（二）实验设备 (2)（三）实验要求 (2)（四）实验内容 (2)1.1 MATLAB基础知识 (2)1.1.1 MATLAB程序设计语言简介 (2)1.1.2 MA TLAB界面及帮助 (2)1.2 MA TLAB基本运算 (4)1.2.1 MA TLAB内部特殊变量和常数 (4)1.2.2 变量类型 (4)1.2.3 内存变量管理 (5)1.2.4 MA TLAB常用数学函数 (5)1.2.5 MA TLAB矩阵生成 (5)1.2.6 MA TLAB矩阵运算 (8)1.2.7 MA TLAB中的矩阵分析 (10)1.3 MA TLAB程序设计 (10)1.3.1 M文件 (10)1.3.2 程序控制结构 (12)实验五MALTAB基础知识（一）实验目的●了解MA TLAB 程序设计语言的基本特点，熟悉MA TLAB软件运行环境●掌握创建、保存、打开m文件及函数的方法●掌握变量等有关概念，具备初步的将一般数学问题转化为对应的计算机模型并进行处理的能力（二）实验设备计算机，Matlab软件（三）实验要求本实验属于验证实验，请根据（四）实验内容的步骤，运行相应的指令或例子，并将仿真结果截图至文档（请自己新建一个word文档，注意，并不一定所有指令或例子的实验结果都要截图，截图数目大于等于5个即可，自己选择性截图，答案不唯一，自由发挥）请在页眉处填写班级、学号、姓名，并将实验报告命名为“实验五_学号_姓名”，并通过FTP上传至指定文件夹。
（四）实验内容1.1 MATLAB基础知识1.1.1 MATLAB程序设计语言简介MA TLAB，Matrix Laboratory的缩写，是由MathWorks公司开发的一套用于科学工程计算的可视化高性能语言，具有强大的矩阵运算能力。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《数学实验》报告
实验名称 MATLAB绘图
学 院 东凌经济管理学院
专业班级 工管111
姓 名 李伊童
学 号 41171034
2013年 5月
2
一、 【实验目的】
插值与拟合
二、 【实验任务】
第六章8、10、12
三、 【实验程序】
8、x=[0.10,0.20,0.30,0.55,0.70,0.80,0.95];
y=[15,18,19,21,22.6,23.8,26];
p1=polyfit(x,y,1);
p3=polyfit(x,y,3);
p5=polyfit(x,y,5);
disp(‘一阶拟合函数'),f1=poly2str(p1,'x');
disp(‘三阶拟合函数'),f3=poly2str(p3,'x');
disp(‘五阶拟合函数'),f5=poly2str(p5,'x');
x1=0.10:0.01:0.95;
y1=polyval(p1,x1);
y3=polyval(p3,x1);
y5=polyval(p5,x1);
plot(x,y,'rp',x1,y1,'--',x1,y3,'k-.',x1,y5);
legend(‘拟合点’,'一次拟合','三次拟合','五次拟合');
10、x=[10,15,20,25,30];
y=[25.2,29.8,31.2,31.7,29.4];
3
xi=10:0.5:30;
yi1=interp1(x,y,xi,'*nearest');
yi2=interp1(x,y,xi,'*linear');
yi3=interp1(x,y,xi,'*spline');
yi4=interp1(x,y,xi,'*cubic');
plot(x,y,'ro',xi,yi1,'--',xi,yi2,'-',xi,yi3,'k.-',xi,yi4,'m:'), grid on;
legend('原始数据','最近值插点','线性插点','样条插点','立方插点');
12、[x,y]=meshgrid(-3:0.3:3);
>> z=x.^2/16-y.^2/9;
>> [x1,y1]=meshgrid(-3:0.1:3);
>> z1=x1.^2/16-y1.^2/9;
>> figure(1)
>> subplot(1,2,1),mesh(x,y,z),title('数据点')
>> subplot(1,2,2),mesh(x1,y1,z1),title('函数图象')
>> [xi,yi]=meshgrid(-3:0.05:3);
>> zi1=interp2(x,y,z,xi,yi,'*nearest');
>> zi2=interp2(x,y,z,xi,yi,'*linear');
>> zi3=interp2(x,y,z,xi,yi,'*spline');
>> zi4=interp2(x,y,z,xi,yi,'*cubic');
>> figure(2)
>> subplot(221),mesh(xi,yi,zi1),title('最近点插值')
>> subplot(222),mesh(xi,yi,zi2),title('线性插值')
4
>> subplot(223),mesh(xi,yi,zi3),title('样条插值')
>> subplot(224),mesh(xi,yi,zi4),title('立方插值')
四、 【实验结果】
8、
10、
5
12、
五、 【实验总结】
此次作业比较综合，与前面二维绘图和三维绘图相关。