Matlab编程训练题(一)
Matlab编程训练题（一）：找出10万以内的所有亲和数对
说明：什么是亲和数？
亲和数是这样一对正整数a和b，使得a的所有真因子的和等于b，而b的所有真因子的和等于a。亲和数问题最早由毕达哥拉斯学派发现和研究的。他们在研究数字的规律的时候发现有以下的性点的两个数：
1+2+4+5+10+11+20+22+44+55+110=284，
1+2+4+71+142=220，
就是220的真因子是1、2、4、5、10、11、20、22、44、55、110，它们的和是284；284的真因子是1、2、4、71、142，其和恰好是220。这是最早发现的一对亲和数，也是最小的一对亲和数。
考虑到1是每个整数的因子，把除去整数本身之外的所有因子叫做这个数的“真因子”。如果两个整数，其中每一个数的真因子的和都恰好等于另一个数，那么这两个数，就构成一对“亲和数”。
3对亲和数有兴趣的人大约在公元9世纪，杰出的阿拉伯数学家本·科拉建立了一个有名的亲和数公式：
设a=3·2x-1，b=3·2x-1-1，c=9·22x-1-1，这里x是大于1的自然数，如果a、b、c 全是素数的话，那么2x·ab与a x·c。便是一对亲和数。
例如，当x=2时，我们不难算出a=11，b=5，c=71，它们全都是素数，所以2x·ab=22·11·5=220；2x·c=22·71=284。
后来的人们对亲和数研究一直保持着极大的兴趣，特别是大数学家费尔马、笛卡尔和欧拉等都曾经研究过亲和数。1636年法国数学家费马发现了第二对亲和数，它们是17962与18416。1638年笛卡尔给出了第三对亲和数。第三对和第四对亲和数，即17926与18416及9363548与9437506。
要数对亲和数做过比较深入研究和为寻找亲和数花了很多功夫的人应当是瑞士的著名数学家欧拉。1747年大数学家欧拉一下子找出了30对，3年后，1750年欧拉向公众宣布了另外的30对亲和数，这样亲和数的数量又增加到了62对，并给出了一个有62对亲和数表。这样大的进展真的给人们一个大的惊喜。可是这样一来，人们反倒觉得既然大数学家欧拉都已经研究过亲和数了，而且他一个人就发现了60对亲和数。欧拉算出了长达几十位、天文数字般的亲和数，那么应该能够计算的数可能都被欧拉找出来了，肯定不会有什么遗漏。
但是，让人没有想到的是，除去最小的220与284之外，另一对亲和数1184与1210竟然被欧拉和另外几位数学大师都漏过了。这对亲和数是在一百多年之后，当“亲和数”的话题不那么热了，似乎已被世人淡忘的时候，1886年一个16岁的意大利男孩帕加尼尼发现这对亲和数，如果把亲和数按从小到大的顺序排列，那么这个少年发现的亲和数是排在第二位。这也可以说明一个现象，就是“百密一疏”，被漏过的恰恰是近在第一对亲和数身旁的第二对1184与1210，最容易的反倒是被人忽略了。
对于亲和数的性质我们知道的还不多，能否用一个公式求出所有的亲和数也不清楚。但是随着计算机的性能不断地提高，利用计算机计算亲和数要比过去容易得多，可以找出更多的亲和数。
目前已经知道的有1000多对亲和数，而10000以内的只有5对，在100000以内有13对，它们是：220和284、1184和1210、2620和2924、5020和5564、6232和6368、10744和10856、12285和14595、17296和18416、63020和76084、66928和66992、67095和71145、69615和87633、79750和88730。
《MATLAB程序设计与应用(刘卫国)》(第二版) 答案
《MATLAB程序设计与应用（刘卫国）》（第二版）实验一MATLAB运算基础 1.（1） z1=2*sin(pi*85/180)/(1+exp(2)) （2）x=[2,1+2i;-0.45,5]; z2=log(x+sqrt(1+x.^2))/2 （3）a=-3.0:0.1:3.0; z3=0.5*(exp(0.3*a)-exp(-0.3*a)).*sin(a+0.3)+log((0.3+a)/2) （4）t=0:0.5:2.5; z4=(t>=0&t<1).*t.^2+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^2-2*t+1) 2. A=[12,34,-4;34,7,87;3,65,7]; B=[1,3,-1;2,0,3;3,-2,7]; （1）a=A+6*B b=A-B+eye(size(A)) %I=eye(size(A)) （2）c=A*B d=A.*B （3）e=A^3 f=A.^3 （4）g=A/B h=B\A （5）m=[A,B] n=[A([1,3],:);B^2] 3. A=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,24,25]; B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11]; （1）C=A*B （2）D=C(3:end,2:end) 4.（1） a=100:999; b=rem(a,21)==0; c=find(b); d=length(c) （2）ch='Just as Bianhaiman said,Xiehong is ...'; e=find(ch>='A'&ch<='Z'); ch(e)=[] 实验二 MATLAB矩阵分析与处理 1.E=eye(3);
MATLAB 语言程序设计基础(2)
%第六章微分方程问题的解法 % 微分方程的解析解方法 % 常微分方程问题的数值解法 % 微分方程问题算法概述 % 四阶定步长Runge-Kutta算法及MATLAB 实现 % 一阶微分方程组的数值解 % 微分方程转换 % 特殊微分方程的数值解 % 边值问题的计算机求解 % 偏微分方程的解 % 6.1 微分方程的解析解方法 % y=dsolve(f1, f2, …, fm ,'x') % syms t; u=exp(-5*t)*cos(2*t+1)+5; % uu=5*diff(u,t,2)+4*diff(u,t)+2*u % syms t y; % y=dsolve(['D4y+10*D3y+35*D2y+50*Dy+24*y='... % '87*exp(-5*t)*cos(2*t+1)+92*exp(-5*t)*sin(2*t+1)+10'],'y(0)=3','Dy(0)=2','D2y(0)=0','D3y(0)=0') % [x,y]=dsolve('D2x+2*Dx=x+2*y-exp(-t)',... % 'Dy=4*x+3*y+4*exp(-t)') % syms t x; % x=dsolve('Dx=x*(1-x^2)+1') % Warning: Explicit solution could not be found; implicit solution returned. % > In D:\MATLAB6p5\toolbox\symbolic\dsolve.m at line 292 % x = % t-Int(1/(a-a^3+1),a=``..x)+C1=0 % 故只有部分非线性微分方程有解析解。 % 6.2 微分方程问题的数值解法 % 6.2.1 微方程问题算法概述 %Euler算法% %function [outx,outy]=MyEuler(fun,x0,xt,y0,PointNum) % fun 表示f(x,y); x0,xt:自变量的初值和终值; % y0:函数在x0处的值,其可以为向量形式; % PointNum表示自变量在[x0,xt]上取的点数 % if nargin<5 |PointNum<=0 %PointNum 默认值为100 % PointNum=100; % end % if nargin<4 % y0默认值为0
刘卫国版MATLAB程序设计与应用课后实验六八九
实验六 高层绘图操作 %第一题： 程序代码如下： x=linspace(0,2*pi,101); y=(0.5+3*sin(x)./(1+x.^2)).*cos(x); plot(x,y) 01234567 -1 -0.5 0.5 1 1.5 %第二题： %（1） 程序代码如下： x=linspace(-2*pi,2*pi,100); y1=x.^2; y2=cos(2*x); y3=y1.*y2; plot(x,y1,'b-',x,y2,'r:',x,y3,'y--'); text(4,16,'\leftarrow y1=x^2'); text(6*pi/4,-1,'\downarrow y2=cos(2*x)'); text(-1.5*pi,-2.25*pi*pi,'\uparrow y3=y1*y2');
-8 -6 -4 -2 2 4 6 8 -30-20 -10 10 20 30 40 %（2） 程序代码如下： x=linspace(-2*pi,2*pi,100); y1=x.^2; y2=cos(2*x); y3=y1.*y2; subplot(1,3,1);%分区 plot(x,y1); title('y1=x^2');%设置标题 subplot(1,3,2); plot(x,y2); title('y2=cos(2*x)'); subplot(1,3,3); plot(x,y3); title('y3=x^2*cos(2*x)');
-10 10 0510 15202530 35 40y1=x 2 -10 10 -1-0.8 -0.6 -0.4-0.200.20.4 0.6 0.8 1y2=cos(2*x) -10 10 -30-20 -10 10 20 30 40 y3=x 2*cos(2*x) %（3） 程序代码如下： x=linspace(-2*pi,2*pi,20); y1=x.^2; subplot(2,2,1);%分区 bar(x,y1); title('y1=x^2的条形图');%设置标题 subplot(2,2,2); stairs(x,y1); title('y1=x^2的阶梯图'); subplot(2,2,3); stem(x,y1); title('y1=x^2的杆图'); subplot(2,2,4); fill(x,y1,'r');%如果少了'r'则会出错 title('y1=x^2的填充图'); %其他的函数照样做。
Matlab语言程序设计
MATLAB 语言程序设计 1、应用题（20 分） 某工厂生产A和B两种产品，它们需要经过三种设备的加工，工时如表1所列。设备一、二、三每天可使用的时间分别不超过12、10和8小时。产品A和B的利润随市场的需求有所波动，如果预测未来某个时期内A和B的利润分别为4000元/吨和3000元/吨，问每天应安排产品A、B各多少吨，才能使工厂获利最大？ 表 1 生产产品工时表 产品设备一设备二设备三 A/（小时/吨） 3 3 4 B/（小时/吨） 4 3 2 设备每天最多可工作时数/小时12 10 8 请写出具体解决方案，并编写该解决方案的具体MATLAB程序代码。 假设每天应安排生产产品A和B分别为和吨，依题意可建数学模型如下： max z=4000x 1+ 3000x2 3x1+ 4x2≤12 3x1+ 3x2≤10 4x1+ 2x2≤8 x1≥0, x2≥0 首先将目标函数转换为标准形式:min z=-4000-3000 M文件如下： f = [-4000;-3000]; A = [3 4;3 3;4 2]; b = [12;10;8]; lb = zeros(2,1); [x,fval] = linprog(f,A,b,[],[],lb) 结果如下图1所示： 图1
2、结合MATLAB 图像的读写知识，请你写出具体MATLAB 图像读写方法，结合实例，给出图像读写步骤。（20分） 读写方法： imfinfo(FileNAme) %从FileNAme中获取或显示图像文件的特征数据 [X,cmap]=imread(FileNAme) %读取变址图像的数据矩阵与伴随色图，并分别送给X 和cmap X=imread(FileNAme) %从FileNAme中读取强度图像或真彩图像的数据矩阵并送给X imwrite(X,cmap,FileName,Parameter,Value) %将变址图像写入文件 imwrite（X ,FileName,Parameter,Value）%将强度图像或真彩图像写入文件 [例2-1] M文件如下： I=imread('m2.png'); imwrite(I,'m2.tif'); [X,cmap]=imread('m2.tif'); image(X); colormap(cmap); axis image off; 结果如下图2所示： 图2 3、结合课本第四章所提供的18种绘图方法，请实现其中8种绘图，提供程序代码及所绘制的图形。（30分）
matlab程序设计与应用第二版习题答案
matlab程序设计与应用第二版习题答案【篇一：matlab程序设计与应用(第二版)实验答案】 %实验一 matlab运算基础 %第1题 %（1） z1=2*sin(85*pi/180)/(1+exp(2)) %（2） x=[2,1+2i;-0.45,5]; z2=0.5*log(x+sqrt(1+x.^2)) %(3) a=-3.0:0.1:3.0; z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2) %(4) t=0:0.5:2.5; z4=t.^2.*(t=0t1)+(t.^2-1).*(t=1t2)+(t.^2-2*t+1).*(t=2t3) %第2题 a=[12 34 -4;34 7 87;3 65 7]; b=[1 3 -1;2 0 3;3 -2 7]; a+6*b a-b+eye(size(a)) a*b a.*b a^3 a.^3 a/b b\a [a,b] [a([1,3],:);b^2] %第3题 a=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25] b=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11] c=a*b f=size(c) d=c(f(1)-2:f(1),f(2)-1:f(2)) whos %第4题 %(1):
a=100:999; b=rem(a,21); c=length(find(b==0)) %(2): a=lsdhksdlkklsdkl; k=find(a=aa=z); a(k)=[] %实验二 matlab矩阵分析与处理 %第1题 e=eye(3); r=rand(3,2); o=zeros(2,3); s=diag([2,3]); a=[e,r;o,s]; a^2 b=[e,(r+r*s);o,s^2] %第2题 h=hilb(5) p=pascal(5) hh=det(h) hp=det(p) th=cond(h) tp=cond(p) %第3题 a=fix(10*rand(5)) h=det(a) trace=trace(a) rank=rank(a) norm=norm(a) %第4题 a=[-29,6,18;20,5,12;-8,8,5] [v,d]=eig(a) %数学意义略 %第5题方法一 %(1): a=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; b=[0.95,0.67,0.52]; x=inv(a)*b %(2):
05级通信工程、电子信息工程《MATLAB语言程序设计》考试卷
安徽农业大学2006―2007学年第二学期 《MATLAB 语言程序设计》试卷（A 卷） 考试形式: 闭卷笔试，2小时 适用专业： 05级通信工程、电子信息工程 一、填空题(本题满分30分，每空3分) 1． 设有程序 A=[1,2,3,4;4,3,2,1;1, -2 ,1, -2];B=[3,1,-1;0,4,2];A1=A(:,[1 3]); M=size(A1)+ length(B(:,2)) 将下列命令的运行结果填在横线上 M= ;A1+B' = . 2. A=[1,2,3,4;2,3,4,5;3,4,5,6]; A1=sum(A<4,2);A([2,3],:)=[],A2=A A1= ; A2= ; 3.P=[1,2,3,4;3,0,1,2];max(P)= ;mean(P,2)= 院： 专业班级： 姓名： 学号： 装 订 线
二、(本题满分12分)试编写计算程序 V; 与特征向量 U 的特征值 B 2 A 求 (4). E; 5 B A 2B BX 解矩阵 . (3) ); A B) (B(A R R 的秩 A (2). |; B) 2 B)(A (A | D (1). , 6 2 1 5 7 2 , 2 1 5 1 2 3 3 2 1 1 1 + + = + = + - = ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? = - - T T T B A 方程 设
三、(本题满分8分) 的程序组写出求解超定线性方程?????? ?=++-=++=--=++1 21212:321 321321 321x x x x x x x x x x x x 四、(本题满分16分)运用符号运算功能写出下列各题Matlab 程序 解求方程03.12=-x e x Adx dx A d x wx x x xe A x ??? ????=-1 022,,cos )sin(sin .2并计算 生成符号矩阵
matlab程序设计实例
MATLAB 程序设计方法及若干程序实例 樊双喜 (河南大学数学与 信息科学学院开封475004) 摘要本文通过对 MATLAB 程序设计中的若干典型问题做简要的分析和总结,并在此基础上着重讨论了有关算法设计、程序的调试与测试、算法与程序的优化以及循环控制等方面的问题.还通过对一些程序实例做具体解析,来方便读者进行编程训练并掌握一些有关MATLAB 程序设计方面的基本概念、基本方法以及某些问题的处理技巧等.此外,在文章的最后还给出了几个常用数学方法的算法程序, 供读者参考使用.希望能对初学者进行 MATLAB 编程训练提供一些可供参考的材料,并起到一定的指导和激励作用,进而为MATLAB 编程入门打下好的基础. 关键字算法设计；程序调试与测试；程序优化；循环控制 1 算法与程序 1.1 算法与程序的关系算法被称为程序的灵魂,因此在介绍程序之前应先了 解什么是算法.所谓算 法就是对特定问题求解步骤的一种描述.对于一个较复杂的计算或是数据处理的问题,通常是先设计出在理论上可行的算法,即程序的操作步骤,然后再按照算法逐步翻译成相应的程序语言,即计算机可识别的语言. 所谓程序设计,就是使用在计算机上可执行的程序代码来有效的描述用于解决特定问题算法的过程.简单来说,程序就是指令的集合.结构化程序设计由于采用了模块分化与功能分解,自顶向下,即分而治之的方法,因而可将一个较复杂的问题分解为若干子问题,逐步求精.算法是操作的过程,而程序结构和程序流程则是算法的具体体现. 1.2MATLAB 语言的特点 MATLAB 语言简洁紧凑,使用方便灵活,库函数极其丰富,其语法规则与科技人员的思维和书写习惯相近,便于操作.MATLAB 程序书写形式自由,利用其丰富
《MATLAB语言程序设计》练习题1(08通信)
《MATLAB语言程序设计》练习题1（2010.5） 一、填空题 1．设有程序 A=[1,2,3,4;6,-5,6,3];B=[2,4,5,7];A1=A(:,[1 3 4]);A2=A(1,[2,3]). 将下列命令的运行结果填在横线上 size(A)= ;A(1,: ).*B= ;length(B)= ; A1= ; A2= ;A+3= ;B.^2= 2. A=[1,2,3,4;6,-5,6,3];B=[-2,4,-6,7,9];A1=(A<3);A2=all(A(:,1)<5);A3=any(A(2,:)<10); A4=find(abs(B)<3| abs(B)>7); A5=find(abs(B)<7& abs(B)>3) A1= ; A2= ; A3= ; A4= ;A5= .
; ).4(; 2(3).);BA R(A R 2).A (|;)(|).1(,21 4154 101 ,32 1214 1531 1 1 V U A E B A B AX AB D B A ：、T T 与特征向量 的特征值 求解矩阵方程 的秩设题试编写程序计算下列各二+=+==???? ? ? ?--=????? ? ?-=--- ??? ?? ?? ??=+-=++-=++=-+=++???????=-++-=-++=-++=+++?? ?=-++=+++3 512332232.31 23322352:.2221 432:.13213213213213214321 43214321432143214321x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x x 、超定方程组 恰定方程组 欠定方程组 程组的程序 写出求解下列各线性方三 .)(,cos 3)sin()sin(cos 2sin )(.3;)(,,1 sin )ln(.2;02.1: ,4 321 2 2 2 2 3 变换的计算并计算生成符号矩阵 的根求方程下列各题运用符号运算功能计算四Laplace t A wt e wt t t te t t t te t t A dx x A dx A d x x e s x x A x x 、ut t t x ? ??? ? ??? ? ?=??? ? ??? ?+==+---? 的解 函数求解微分方程用1)0(')0(,sin '2".4===++x x t x x x dsolve （1） 若.2 c bt at y t y ++=的经验公式为 与试编写程序计算出上式中的a 、b 、c;
MATLAB 语言程序设计基础(1)
%第二章MATLAB 语言程序设计基础% % MATLAB 基本命令简介 % MATLAB 程序设计语言基础 % 基本数学运算 % MATLAB语言流程控制 % MATLAB 函数的编写 % 二维图形绘制 % 三维图形绘制 % 二维图形 % 基本平面 % plot % fplot%f(x) %fplot('tanh',[-2 2]) % loglog%双对数图形 % semilogx%x轴对数图形 % semilogy %zoom %meshgrid
% 特殊平面 %polar %bar %barh %compass %comet %errorbar %feather %hist 二维直方图%histc 直方图记数%rose 角度直方图%stairs %stem 柄形图 %stem3 %pie 饼形图 % 注释命令 % grid % gtext % text % legend % title
% xlabel,ylabel % 三维图形 % % 三维曲线，面填色命令 % comet3 三维彗星 % fill3 % 三维图形等高线 %clabel 二维等高线图中添加高度标签 % [x,y] = meshgrid(-2:.2:2); % z = x.*y.*exp(-x.^2-y.^2); % [C,h] = contour(x,y,z); % clabel(C,h); %contour %contourc 低级等高线图形计算命令。计算等高线矩阵c %contour3 三维空间等高线图 %contourf 填充二维等高线图 %contour，contour3和contourf %pie3
大学毕业设计-MATLAB语言程序设计
MATLAB语言程序设计 系别：电子电气工程系 班级：08级自动化班 姓名：XXX 学号：
1.编写一个求圆的面积的函数文件。 >> f(1) s = 3.1416 >> f(2) s = 12.5664 通过此题掌握了独立文件与函数文件的区别。 2.三次抛物线的方程为：y=a*x^3+b*x^2+c*x 试探讨参数a，b和c对其图形的影响。 >> clear all; >> subplot(1,3,1) >> fplot('(-2:2)*x.^3+x.^2+x',[-2 2]) %绘制变量a=-2,-1,0,1,2时的图形 >> grid,axis('equal'),axis([-2 2 -4 4]) %显示网格，纵横坐标轴保持一致图 >> %显示区间 >> subplot(1,3,2) >> fplot('x.^3+(-2:2)x.^2+x',[-2 2]) %绘制变量b=-2,-1,0,1,2时的图形 >> fplot('x.^3+(-2:2)*x.^2+x',[-2 2]) %绘制变量b=-2,-1,0,1,2时的图形 >> grid,axis('equal'),axis([-2 2 -4 4]) >> subplot(1,3,3) >> fplot('x.^3+x.^2+(-2:2)*x',[-2 2]) %绘制变量c=-2,-1,0,1,2时的图形 >> grid,axis('equal'),axis([-2 2 -4 4]) >> gtext('a=2'),gtext('a=-2'),gtext('a=0') %用gtext命令在其他两个子图上标注字符程序运行结果如下：
MATLAB程序设计及应用(第二版)课后实验答案
。 Matlab课后实验题答案 《 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 (1) 0 12 2sin85 1 z e = + (2) 21 ln( 2 z x =，其中 212 0.455 i x + ??=?? -?? (3) 0.30.3 3 0.3 sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22 a a e e a z a a - -+ =++=--
(4) 22 42011 122123t t z t t t t t ?≤
) 3. 设有矩阵A 和B 1234 53 166789101769,11 121314150 23416171819209 7021222324254 1311A B ???? ????-??? ?????==-??? ? ???????????? (1) 求它们的乘积C 。 (2) 将矩阵C 的右下角3×2子矩阵赋给D 。 (3) 查看MATLAB 工作空间的使用情况。 4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 解：(1) 结果：
MATLAB软件语言及程序设计报告(含源程序版
MATLAB软件语言及程序设计 姓名: 学号: 班级: 院系: 指导老师:张登峰 时间:2010-5-10
一、问答题(15分) 1. 在安装MATLAB软件时，哪些组件(模块)是必须选择的？否则，不能建立MATLAB工作环境。 Toolboxes、MATLAB、Complier 2. 写出MA TLAB中的6个预定义变量名，并说明各自的含义。 ans:计算结果的缺省变量名 pi:圆周率 i,j：虚数单位 inf：无穷大 eps：计算机的最小数 NaN：不定量，如0/0，00/00 nargin：函数的输入变量个数 nargout：函数的输出变量个数 realmin：最小正实数 realmax：最大正实数 lasterr：存放最新的错误信息 lastwarn：存放最新的警告信息 3. 在MA TLAB指令窗中运行指令(-243)^(1/5)后，会得到-3吗？ 的全部方根有几 个，请写出计算全部方根的M文件程序代码，要求对每条语句加注释说明。 不会，会得到2.4271 + 1.7634i，用nthroot(-243,5)可以得到一个根-3。 可以用solve解决，程序如下 solve('t^5+243') 得到结果为： ans = -3 -3/4*5^(1/2)+3/4-3/4*i*2^(1/2)*(5+5^(1/2))^(1/2) 3/4*5^(1/2)+3/4-3/4*i*2^(1/2)*(5-5^(1/2))^(1/2) 3/4*5^(1/2)+3/4+3/4*i*2^(1/2)*(5-5^(1/2))^(1/2) -3/4*5^(1/2)+3/4+3/4*i*2^(1/2)*(5+5^(1/2))^(1/2) 4. 指令clear，clf，clc各有什么用处，它们运行结果的区别是什么？ Clear ：清除内存中所有的或指定的变量和函数 Clf：擦除MATLAB当前工作窗口中的图形 Clc：擦除MATLAB工作窗口中所有显示的内容 5. 在MATLAB指令表达式中，标点符号“空格”和“逗号”各有什么作用？请举例说明，在哪些情况下，两者的作用相同？在哪些情况下，二者的作用不同？ 空格和逗号都能起到分隔符的作用，逗号可作为函数分割符，也用于区分行，显示运算结果，当然不加标点也显示运算结果，而在大多数情况下，MA TLAB对空格不予处理。
MATLAB程序设计与应用(第二版)
阅阅读读时时：：请请选选择择““视视图图||文文档档结结构构图图””，，弹弹出出文文档档中中的的标标题题链链接接。。数学软件 MATLAB 程序设计与应用
第1章MATLAB系统环境 1.1 MATLAB概貌 1.2 MATLAB环境的准备 1.3 MATLAB操作界面 1.4 MATLAB帮助系统 自上世纪80年代以来，出现了科学计算语言，亦称数学软件。 MATLAB Mathematica Mathcad Maple LINDO LINGO 1.1 MATLAB概貌p3 MATLAB 是MATrix LABoratory（矩阵实验室）的缩写。 1984年由Math Works 公司推出，现已成为国际公认的优秀的工程应用开发环境，是影响最大，流行最广的科学计算语言。 1.1.1 MATLAB的发展
1. 从MATLAB 4.2c开始，每个版本增加了一个建造编号； 2. 例如MATLAB7.6的建造编号是R2008a。说明MATLAB7.6与MATLAB2008a是等同的； 3. 对于建造编号，正规化以后，每年出两个版本。一般来说。a是测试版，b是正式版。a是前半 年出，b是后半年出。 教材采用MATLAB7.0（R14，2004） 实验室采用MATLAB 7.8（R2009a，2009.3，汉化） 1.1.2 MATLAB的主要功能p4 ◆数值计算和符号计算功能 ◆绘图功能 ◆语言体系 ◆MATLAB工具箱 (1) 数值计算和符号计算功能 MATLAB以矩阵作为数据操作的基本单位，还提供了十分丰富的数值计算函数。 MATLAB先后和著名的符号计算语言Maple与MuPAD（从MATLAB 2008b开始使用MuPAD）相结合，使得MATLAB具有符号计算功能。 (2) 绘图功能 可以绘制二维和三维图形。 MATLAB提供了两个层次的绘图操作： ●对图形句柄进行的低层绘图操作； ●建立在低层绘图操作之上的高层绘图操作。 (3) 语言体系 MATLAB具有程序结构控制、函数调用、数据结构、输入输出、面向对象等程序语言特征，而且简单易学、编程效率高。 MATLAB是解释性语言，不能脱离MATLAB环境而独立运行。 (4) MATLAB工具箱 MATLAB包含两部分内容：基本部分和各种可选的工具箱。 MATLAB工具箱分为两大类：功能性工具箱和学科性工具箱。 基本部分 构成MATLAB的核心内容，也是使用和构造工具箱的基础。 功能性工具箱 主要用来扩充其符号计算功能、可视建模仿真功能及文字处理功能等。 学科性工具箱 Control System Toolbox 控制系统工具箱 Signal Processing Toolbox 信号处理工具箱 Neural Network Toolbox 神经网络工具箱 Optimization Toolbox 最优化工具箱 Financial Toolbox 金融工具箱 Statistics Toolbox 统计学工具箱 开始→工具箱 MATLAB具备很强的开放性 除内部函数外，所有MATLAB基本文件和各工具箱文件都是可读、可改的源文件，用户可通过对源文件的修改或加入自己编写的文件去构成新的专用工具箱。
实验1 Matlab程序设计基础
实验1 Matlab 程序设计基础 一、实验名称：Matlab 程序设计基础. 二、实验目的：掌握Matlab 语言，会用Matlab 进行一些简单的程序设计. 三、实验要求： 1、对以下问题，编写M 文件： （1）用起泡法对10个数由小到大排序.即将相邻两个数比较，将小的调到前头. （2）有一个4×5矩阵，编程求其最大值及其所处的位置. （3）编程求∑=20 1!n n . （4）一球从100米高度自由落下，每次落地后反跳回原高度的一半，再落下.求它在第10次落地时，共经过多少米？第10次反弹有多高？ （5）有一函数f(x,y)=x^2+sin(x*y)+2*y ，写一程序，输入自变量的值，输出函数值. 2、在同一平面中的两个窗口分别画出心形线和马鞍面. 要求： （1）在图形上加格栅、图例和标注. （2）定制坐标. （3）以不同角度观察马鞍面. 3、请编写绘制以下图形的Matlab 命令，并展示绘得的图形： （1）4,12222=+=+y x y x 分别是椭圆14/2 2=+y x 的内切圆和外切圆； （2）指数函数和对数函数x y ln =的图像关于直线x y =对称. 四、实验内容： 1、1）源程序： function p=paixu(a) a=[1,2,5,6,8,7,9,4,3,0]; for i=9:-1:1 for j=1:i if a(j)>a(j+1) t=a(j); a(j)=a(j+1); a(j+1)=t; end end end a 调试结果：a = 0 1 2 3 4 5 6 7 8 9 2）源程序： function j=juzhen(a)
MATLAB程序设计与应用(第二版)实验参考答案
MATLAB程序设计与应用（第二版）实验参考答案 %实验一MATLAB运算基础 %第一题 %（1） z1=2*sin(85*pi/180)/(1+exp(2)) %（2） x=[2,1+2i;-0.45,5]; z2=0.5*log(x+sqrt(1+x.^2)) %(3) a=-3.0:0.1:3.0; z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2) %(4) t=0:0.5:2.5; z4=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3) %第二题 A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; A+6*B A-B+eye(size(A)) A*B A.*B A^3 A.^3 A/B B\A [A,B] [A([1,3],:);B^2] %第三题 A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25] B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11] C=A*B F=size(C) D=C(F(1)-2:F(1),F(2)-1:F(2)) whos %第四题 %(1): A=100:999; B=rem(A,21); C=length(find(B==0)) %(2): A='lsdhKSDLKklsdkl';
k=find(A>='A'&A<='Z'); A(k)=[] %实验二MATLAB矩阵分析与处理 %第一题 E=eye(3); R=rand(3,2); O=zeros(2,3); S=diag([2,3]); A=[E,R;O,S]; A^2 B=[E,(R+R*S);O,S^2] %第二题 H=hilb(5) P=pascal(5) Hh=det(H) Hp=det(P) Th=cond(H) Tp=cond(P) %第三题: A=fix(10*rand(5)) H=det(A) Trace=trace(A) Rank=rank(A) Norm=norm(A) %第四题: A=[-29,6,18;20,5,12;-8,8,5] [V,D]=eig(A) %数学意义略 %第五题方法一: %(1): A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; b=[0.95,0.67,0.52]'; x=inv(A)*b %(2): B=[0.95,0.67,0.53]'; x=inv(A)*B %(3): cond(A) %第五题方法二: A=hilb(4) A(:,1)=[] A(4,:)=[] B=[0.95,0.67,0.52]';
(完整word版)MATLAB程序设计及应用(第二版)课后实验答案
Matlab课后实验题答案 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 (1) 0 12 2sin85 1 z e = + (2) 21 ln( 2 z x =+，其中 212 0.455 i x + ??=?? -?? (3) 0.30.3 3 0.3 sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22 a a e e a z a a - -+ =++=--L (4) 2 2 4 2 01 112 2123 t t z t t t t t ?≤< ? =-≤< ? ?-+≤< ? ，其中t=0:0.5:2.5 2. 已知：
1234413134787,2033657327A 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 和B 1234 53 166789101769,11 121314150 23416171819209 7021222324254 1311A B ???? ????-??? ?????==-??? ? ???????????? (1) 求它们的乘积C 。 (2) 将矩阵C 的右下角3×2子矩阵赋给D 。 (3) 查看MATLAB 工作空间的使用情况。
4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 解：(1) 结果： (2). 建立一个字符串向量 例如： ch='ABC123d4e56Fg9';则要求结果是： 实验二 MATLAB 矩阵分析与处理 1. 设有分块矩阵33 322322E R A O S ?????? =? ??? ，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩
MATLAB程序设计复习题(最新整理)
《MATLAB程序设计》复习题 一填空题 1 已知A=[ 2 5 6;1 7 8];B=[1 4 3;4 7 0];写出下列各指令运行的结果。 A & B 的运行结果为 [1 1 1;1 1 0] 。 A== B的运行结果为 [0 0 0 ;0 1 0] 。 2 产生3阶单位阵的命令为 eye(3) ；产生4阶全1方阵的命令为ones(4)。 3 A=rand(5,8);b=size(A);c=length(A)；则b和c的值分别为[5,8]和8。 4 求矩阵的特征值的指令为 eig 。 5符号；可以使命令行不显示运算结果， % 用来表示该行为注释行。6在线形图型格式的设置中，字符r表示什么颜色（红色） 7fix函数四舍五入对数组[3.17，5.34，8.90，2.52]取整，结果为 [3 5 8 2] 。 8 P, Q分别是个多项式的系数矢量，求P对应的多项式的根，使用的命令是 roots(P) 9使2个plot的图形在同一个坐标显示，使用 hold on 命令进行图形保持，使用 grid on 命令为图形添加网格。 10 数据拟合采用___polyfit___函数；求不定积分采用___ int ___函数。 11 A=rand(2,5);b=size(A);c=length(A)；则b和c的值分别为_2 5和5__ 12已知A是矩阵，求A的对角矩阵函数是___diag(A)____，求A的下三角矩阵函数是__tril(A)____。 13条形图由函数__bar___来实现。 14已知s=‘显示“hey”’，则s的元素个数是9 二解答题 1、A=[71，3，-8; 2，-9，8; 0，4，5]
MATLAB课程设计报告(绝对完整)
课程设计任务书 学生姓名：董航专业班级：电信1006班 指导教师：阙大顺，李景松工作单位：信息工程学院 课程设计名称：Matlab应用课程设计 课程设计题目：Matlab运算与应用设计5 初始条件： 以上版本软件； 课程设计辅导资料：“Matlab语言基础及使用入门”、“Matlab及在电子信息课程中的应用”、线性代数及相关书籍等； 先修课程：高等数学、线性代数、电路、Matlab应用实践及信号处理类相关课程等。 要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） 课程设计内容：根据指导老师给定的7套题目，按规定选择其中1套完成； 本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析，针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结。具体设计要求包括： 初步了解Matlab、熟悉Matlab界面、进行简单操作； MATLAB的数值计算：创建矩阵矩阵运算、多项式运算、线性方程组、数值统计； 基本绘图函数：plot, plot3, mesh, surf等，要求掌握以上绘图函数的用法、简单图形标注、简单颜色设定等； 使用文本编辑器编辑m文件，函数调用； 能进行简单的信号处理Matlab编程； 按要求参加课程设计实验演示和答辩等。 课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括： 目录； 与设计题目相关的理论分析、归纳和总结； 与设计内容相关的原理分析、建模、推导、可行性分析； 程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结； 课程设计的心得体会（至少500字）； 参考文献（不少于5篇）； 其它必要内容等。 时间安排：周（分散进行） 参考文献： (美)穆尔,高会生,刘童娜,李聪聪．MATLAB实用教程(第二版) . 电子工业出版社,2010. 王正林,刘明．精通MATLAB(升级版) .电子工业出版社,2011. 陈杰. MATLAB宝典(第3版) . 电子工业出版社,2011. 刘保柱,苏彦华,张宏林. MATLAB 从入门到精通(修订版) . 人民邮电出版社,2010. 指导教师签名：年月日 系主任（或责任教师）签名：年月日
MATLAB程序设计与应用课后习题答案
西安科技大学MATLAB程序设计 专业：信息与计算科学 班级：1001班 学号：1008060129 姓名：刘仲能 2012年6月27日
实验一 2.已知： ?? ??? ?????-=76538773443412A ，???? ? ?????--=723302131B 求下列表达式的值： （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 和B ????? ??? ??? ?? ???=252423222120191817161514 13121110987654321A ，??????? ? ????????--=11134079423096171603B （1） 求它们的乘积C 。 （2） 将矩阵C 的右下角3×2子矩阵赋给D 。 （3） 查看MATLAB 工作空间的使用情况 (1) (2) （3）
4.完成下列操作 （1）求[100,999]之间能被21整除的数的个数。（2）建立一个字符串向量，删除其中的大写字母。 (1)(2)
实验二 3.建立一个5×5矩阵，求它的行列式值、迹、秩和范数。 运行截图： A 矩阵的行列式值、迹、秩分别如下： 范数如下： 4.已知 ?? ?? ? ?????--=5881252018629A 求A 的特征值及特征向量，并分析其数学意义。 运行截图：
5.下面是一个线性方程组：???? ? ?????=????????????????????52.067.095.06/15/14/15/14/13/14/13/12/1321x x x （1） 求方程的解； （2） 将方程右边向量元素 改为0.53，在求解，并比较 的变化和解的相对 变化； （3） 计算系数矩阵A 的条件数并分析结论。 （2） 变大，其解中，相对未变化前的 的解：x1变大，x2变小，x3变大。 （3） 由于A 矩阵的条件数很大，故当线性方程组中的b 变大时，x 也将发生很大的变
《MATLAB7.x》程序设计语言(第二版)课后题答案
《MATLAB7.x》程序设计语言（第二版）部分课后习题答案（楼顺天，姚若玉，沈俊霞编著） 说明：所有答案均是本人在备考过程中亲自整理的，收录了一部分题目的答案。答案可能存在不足甚至谬误，很多算法也可能不是最优的，仅供参考。本人尽力整理的一点点心血，希望对大家有所帮助。 第二章 8、 a=rand(5,5); [i,j]=find(a>0.5); for u=1:length(i); b(u)=a(i(u),j(u)); end b' 10、aa=any(a'); i=find(aa==0) a(i,:)=[]; 第四章 1、function flag=isprime(m) %m为素数，flag=1 %m不为素数，flag=0 k=sqrt(m); flag=1; for i=2:k; if rem(m,i)==0; flag=0; break end end m=input('input an integer\n'); flag=isprime(m); if flag==1 disp([num2str(m),'是素数']) else disp([num2str(m),'不是素数']) end 4、function [x1,x2]=jfc(a,b,c) d=b^2-4*a*c; if d>0; x1=(-b-sqrt(d))/(2*a); x2=(-b+sqrt(d))/(2*a); elseif d==0; x1=-b/(2*a); x2=x1; end
a=input('a='); b=input('b='); c=input('c='); [x1,x2]=jfc(a,b,c); if x1~=x2; disp('原方程存在两个不同的根') disp(['x1=',num2str(x1)]) disp(['x2=',num2str(x2)]); elseif x1==x2; disp('原方程存在两个相同的根') disp(['x1=x2=',num2str(x1)]) else disp('原方程的根不存在') end %MATLAB上机作业，184页习题4，求二次方程的实根 function [s1,s2]=solve(a,b,c) d=b^2-4*a*c; if d>0 s1=(-b-sqrt(d))/(2*a) s1=(-b-sqrt(d))/(2*a) %disp(['原方程有两个不同的根'num2str(s1)'和'nums2tr(s2)]) elseif d==0 s1=(-b-sqrt(d))/(2*a); s1=(-b-sqrt(d))/(2*a); %disp(['原方程有两个相同的根'num2str(s1)) else disp('原方程无实根') end end 5、x=-3:0.01:3; if x>=-3 & x<-1; y=(-x.^2-4*x-3)/2; elseif x>=-1 & x<1; y=-x.^2+1; else y=(-x.^2+4*x-3)/2; end plot(x,y) 6、%MATLAB上机作业，184页习题6，点不同按钮产生不同分布的数 s=menu('请选择随机数类型','U[-10,10]','U[-5,5]','U[-1,1]','N(0,1)'); switch s case 1,n=(rand(1)-0.5)*20; case 2,n=(rand(1)-0.5)*10; case 3,n=(rand(1)-0.5)*2; case 4,n=randn(1); otherwise disp('请选择！')