matlab判断素数程序设计题

matlab程序设计题

用一个MATLAB语言编写一个程序：输入一个自然数，判断它是否是素数，如果是，输出“It is one prime”，如果不是，输出“It is not one prime.”。要求通过调用子函数实现。最好能具有如下功能：①设计较好的人机对话界面，程序中

含有提示性的输入输出语句。②能实现循环操作，由操作者输入相关命令来控制是否继续进行素数的判断。如果操作者希望停止这种判断，则可以退出程序。③

如果所输入的自然数是一个合数，除了给出其不是素数的结论外，还应给出至少一种其因数分解形式。例：输入6，因为6不是素数。则程序中除了有“It is not one prime”的结论外，还应有：“6=2*3”的说明。

【程序】

a=input('请输入一个整数(end in 0):');

while(a~=0)

if a==1

disp('1 is not a prime or a composite number.')

elseif isprime(a)==1

fprintf('%d is a prime.\n',a)

elseif isprime(a)==0

f=factor(a);m=length(f);

fprintf('%d is not a prime,%d=%d',a,a,f(1))

for i=2:m

fprintf('*%d',f(i));

end

fprintf('\n');

end

a=input('请输入一个整数(end in 0):');

end

实验二 Matlab程序设计基本方法1

实验二Matlab程序设计基本方法 覃照乘自092 电气工程学院 一、实验目的： 1、熟悉MATLAB 程序编辑与设计环境 2、掌握各种编程语句语法规则及程序设计方法 3、函数文件的编写和设计 4、了解和熟悉跨空间变量传递和赋值 二、实验基本知识： ◆for循环结构 语法：for i＝初值：增量：终值 语句1 …… 语句n end 说明：1．i＝初值：终值，则增量为1。 2．初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑。 ◆while 循环结构 语法：while 逻辑表达式 循环体语句 end 说明：1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。若表达式的值为真，执行循环体语句一次、在反复执行时，每次都要进行判断。若表达 式的值为假，则程序执行end之后的语句。 2、为了避免因逻辑上的失误，而陷入死循环，建议在循环体语句的适当位 置加break语句、以便程序能正常执行。（执行循环体的次数不确定； 每一次执行循环体后，一定会改变while后面所跟关系式的值。） 3、while循环也可以嵌套、其结构如下：

while逻辑表达式1 循环体语句1 while逻辑表达式2 循环体语句2 end 循环体语句3 end ◆if-else-end分支结构 if 表达式1 语句1 else if 表达式2(可选) 语句2 else(可选) 语句3 end end 说明：1.if结构是一个条件分支语句，若满足表达式的条件，则往下执行；若不满足，则跳出if结构。 2．else if表达式2与else为可选项，这两条语句可依据具体情况取舍。 3．注意：每一个if都对应一个end，即有几个if，记就应有几个end。 ◆switch-case结构 语法：switch表达式 case常量表达式1 语句组1 case常量表达式2 语句组2 …… otherwise 语句组n end

《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程序设计实践-牛顿法解非线性方程

中南大学MATLAB程序设计实践学长有爱奉献，下载填上信息即可上交，没有下载券的自行百度。所需m文件照本文档做即可，即新建（FILE）→脚本（NEW-Sscript）→复制本文档代码→运行（会跳出保存界面，文件名默认不要修改，保存）→结果。第一题需要把数据文本文档和m文件放在一起。全部测试无误，放心使用。本文档针对做牛顿法求非线性函数题目的同学，当然第一题都一样，所有人都可以用。←记得删掉这段话 班级： ？ 学号： 姓名：

一、《MATLAB程序设计实践》Matlab基础 表示多晶体材料织构的三维取向分布函数（f＝f（φ1，φ，φ2））是一个非常复杂的函数，难以精确的用解析函数表达，通常采用离散 空间函数值来表示取向分布函数，是三维取向分布函数的一个实例。 由于数据量非常大，不便于分析，需要借助图形来分析。请你编写一 个matlab程序画出如下的几种图形来分析其取向分布特征： （1）用Slice函数给出其整体分布特征； " ~ （2）用pcolor或contour函数分别给出(φ2＝0, 5, 10, 15, 20, 25, 30, 35 … 90)切面上f分布情况(需要用到subplot函数)；

(3) 用plot函数给出沿α取向线(φ1=0~90，φ＝45，φ2＝0)的f分布情况。 (

备注：数据格式说明 解： （1）（ （2）将文件内的数据按照要求读取到矩阵f(phi1,phi,phi2)中，代码如 下： fid=fopen(''); for i=1:18 tline=fgetl(fid); end phi1=1;phi=1;phi2=1;line=0; f=zeros(19,19,19); [ while ~feof(fid) tline=fgetl(fid); data=str2num(tline); line=line+1;数据说明部分，与 作图无关此方向表示f随着 φ1从0,5,10,15, 20 …到90的变化而 变化 此方向表示f随着φ 从0,5,10,15, 20 … 到90的变化而变化 表示以下数据为φ2＝0的数据，即f（φ1，φ，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分）

实验3MATLAB程序设计

1，编写M 函数实现求一个数是否为素数，再编写一主程序（脚本文件），要求通过键盘输入一个整数，然后调用判断素数函数，从而确定它是否素数。 x=input('请输入一个整数x:'); if myprime(x) disp('您输入的整数x是一个素数。') else disp('您输入的数x不是一个素数。') end function y=myprime(x) y=1; for i=2:fix(sqrt(x)) if mod(x,i)==0 y=0; end end 2，编写M 函数统计一数值中零的个数，然后编写脚本文件，实现统计从1—2007 中零的总个数。 function num=number0(a) %统计十进制数值中0的个数 sa=num2str(a);%将数值装化为字符串 num=length(find(sa=='0'));% ));%求取字符串中'0’的个数 y=0;

for a=1:2006 num=number0(a); y=num+y; end disp(y) 504 3，编写程序计算x∈[-3,3],字长0.01：并画出曲线x = -3:0.01:3; y=zeros(size(x)); for i = 1:length(x) if -3<= x(i)& x(i)<=-1 y(i)=(-x(i).^2-4*x(i)-3)/ 2; elseif -1<= x(i) & x(i)<=1 y(i)=-x(i).^2+1; elseif 1<=x(:,i)<=3 y(i)=(-x(i).^2+4*x(i)-3)/2; end end plot(x,y) -3-2-10123

1实验一 Matlab程序设计与M文件(1)

实验一 Matlab 程序设计与M 文件 一、实验目的 1. 掌握Matlab 程序设计常用命令，如，循环、选择、暂停、显示输出、输入变量值等。 2. 掌握Matlab 的M 文件的创建于使用，包括脚本文件和函数文件。 3. 掌握常用的编程技巧。 二、实验内容 1. 输入一个百分制成绩，要求输出成绩等级A 、B 、C 、D 、E 。其中90分~100分为A ，80分~89分为B ，79分~79分为C ，60分~69分为D ，60分以下为E 。 要求： (1) 分别用if 语句和switch 语句实现。 (2) 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 2. 硅谷公司员工的工资计算方法如下： (1) 工作时数超过120小时者，超过部分加发15%。 (2) 工作时数低于60小时者，扣发700元。 (3) 其余按每小时84元计发。 试编程按输入的工号和该号员工的工时数，计算应发工资。 3. 根据2 222211116123n π=++++，求π的近似值。当n 分别取100、1000、10000时，结果是多少？ 要求：分别用循环结构和向量运算（使用sum 函数）来实现。 4. 考虑以下迭代公式： 1n n a x b x +=+ 其中a 、b 为正的学数。 (1) 编写程序求迭代的结果，迭代的终止条件为|x n+1-x n |≤10-5，迭代初值x 0=1.0，迭代次数不超过500次。 (2) 如果迭代过程收敛于r ，那么r 的准确值是，当(a,b)的值 取(1,1)、(8,3)、(10,0.1)时，分别对迭代结果和准确值进行比较。

5. 已知 12312311021 323 n n n n f n f n f n f f f f n ---==??==??==??=-+>? 求f 1~f 100中： (1) 最大值、最小值、各数之和。 (2) 正数、零、负数的个数。 6. 若两个连续自然数的乘积减1是素数，则称这两个边疆自然数是亲密数对，该素数是亲密素数。例如，2×3-1=5，由于5是素数，所以2和3是亲密数，5是亲密素数。求[2,50]区间内： (1) 亲密数对的对数。 (2) 与上述亲密数对对应的所有亲密素数之和。 7. 设2411()(2)0.1(3)0.01 f x x x =+-+-+，编写一个MATLAB 函数文件fx.m ，使得调用f(x)时，x 可用矩阵代入，得出的f(x)为同阶矩阵。 8. 一物理系统可用下列方程组来表示： 111211 12 220cos sin 0sin 0cos 000sin 000cos 1a m m a m g m N m N m g θθθθθθ--??????????????????=??????-??????-?????? 从键盘输入m 1、m 2和θ的值，求a 1、a 2、N 1和N 2的值。其中g 取9.8，输入θ时以角度为单位。 要求：定义一个求解线性方程组AX=B 的函数文件，然后在命令文件中调用该函数文件。

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 程序设计实验报告 一、实验目的 1. 通过实验熟悉MATLAB 仿真软件的使用方法； 2. 掌握用MATLAB 对连续信号时域分析、频域分析和s 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征； 3. 掌握用MATLAB 对离散信号时域分析、频域分析和z 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征； 4. 通过绘制信号运算结果的波形，了解这些信号运算对信号所起的作用。 二、实验设备 1. 计算机 ： 2. MATLAB R2007a 仿真软件 三、实验原理 对系统的时域分析 信号的时域运算包括信号的相加、相乘，信号的时域变换包括信号的平移、反折、倒相及信号的尺度变换。 （1）信号的相加和相乘：已知信号)(1t f 和)(2t f ，信号相加和相乘记为 )()(1t f t f =)(2t f +；)()(1 t f t f =)(2t f *。 （2）信号的微分和积分：对于连续时间信号，其微分运算是用diff 函数来完成的，其语句格式为：diff(function,’variable’,n)，其中function 表示需要进行求导运算的信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；n 为求导的阶数，默认值为求一阶导数。连续信号的积分运算用int 函数来完成，语句格式为：diff(function,’variable’,a,b)，其中function 表示需要进行被积信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；a,b 为积分上、下限，a 和b 省略时为求不定积分。 （3）信号的平移、翻转和尺度变换 信号的平移包含信号的左移与右移，信号的翻转包含信号的倒相与折叠，平移和翻转信号不会改变信号)(t f 的面积和能量。信号的尺度变换是对信号)(t f 在时间轴上的变化，可使信号压缩或扩展。)(at f 将原波形压缩a 倍，)/(a t f 将原波形扩大a 倍。 ￥ 对系统频率特性的分析

matlab程序设计作业

Matlab程序设计作业 姓名: 学号: 专业:

? MATLAB 程序设计》作业 1、考虑如下x-y 一组实验数据: x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y 二[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot 的原始数据、一次拟合曲线和三次拟合曲线，给出 原始曲线 MATLAB 代码和运行结果。 7 6 5 4 3 2 2 3 4 5 6 7 8 9 10

7 6.5 6 5.5 5 4.5 4 3.5 3 2.5 10 一次拟合 三次拟合

x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]; y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2]; figure; plot(x,y) p1=polyfit(x,y,1); y1=polyval(p1,x); figure; plot(x,y1) p2=polyfit(x,y,3); y2=polyval(p2,x); figure; plot(x,y2) 2、在［0, 3n区间，绘制y二Sin(x)曲线(要求消去负半波，即(n 2n)区间内的函数值置零)，求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 x=0:pi/1000:3*pi; y=Sin(x); y1=(y>=0).*y; %消去负半波figure(1); plot(x,y1, 'b' ); a=mean(y1) %求出y1 的平均值 b=max(y1) %求出y1 的最大值b, 以及最大值在矩阵中的位置; d=x(find(y1==b)) >> ex1 a = 0.4243 b = 1 d = 1.5708 7.8540 >>

Matlab程序设计实验报告

实验七Matlab 程序设计 实验目的： 1、掌握建立和执行M 文件的方法； 2、掌握实现选择结构的方法； 3、掌握实现循环结构的方法。 实验内容： 1. 编写用 5 次多项式拟合函数y=sin(x), x [0, 2 ]的脚本M 文件，要求绘图观察拟合的效果。 function shiyan1 x=0:0.5:2*pi y=sin(x) p=polyfit(x,y,5) x1=0:0.2:2*pi y1=polyval(p,x1) plot(x,y, 'b' ,x1,y1, '*r' x =

Columns 1 through 9 0 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000 Columns 10 through 13 4.5000 5.0000 5.5000 6.0000 y = Columns 1 through 9 0 0.4794 0.8415 0.9975 0.9093 0.5985 0.1411 -0.3508 -0.7568 Columns 10 through 13 -0.9775 -0.9589 -0.7055 -0.2794 p = -0.0056 0.0881 -0.3967 0.2671 0.8902 0.0029 x1 = Columns 1 through 10 0 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 1.4000 1.6000 1.8000 Columns 11 through 20

2. 2.2000 2.4000 2.6000 2.8000 3.0000 3.2000 3.4000 3.6000 1.8001 Columns 21 through 30 4.0 4.2000 4.4000 4.6000 4.8000 5.0000 5.2000 5.4000 5.6000 5.8000 Columns 31 through 32 6.0 6.2000 y1 = Columns 1 through 10 0.29 0.1886 0.3786 0.5585 0.7172 0.8461 0.9391 0.9926 1.0048 0.9761 Columns 11 through 20 0.9083 0.8048 0.6701 0.5098 0.3301 0.1381 -0.0590 -0.2538 -0.4389 -0.6073 Columns 21 through 30 -0.7524 -0.8685 -0.9505 -0.9949 -0.9991 -0.9626 -0.8863 -0.7732 -0.6288 -0.4606 Columns 31 through 32

matlab程序设计与应用(第二版)第三章部分课后答案

第三章1. (1)A=eye(3) (2)A=100+100*rand(5,6) (3)A=1+sqrt(0.2)*randn(10,50) (4)B=ones(size(A)) (5)A+30*eye(size(A)) (6)B=diag(diag(A)) 2. B=rot90(A) C=rot90(A,-1) 3. B=inv(A) ;A的逆矩阵 C=det(A) ;A的行列式的值 D=A*B E=B*A D=E 因此A与A-1是互逆的。 4. A=[4 2 -1;3 -1 2;12 3 0]; b=[2;10;8]; x=inv(A)*b x = -6.0000 26.6667 27.3333 5. (1) diag(A) ；主对角线元素 ans = 1 1 5 9 triu(A) ；上三角阵

ans = 1 -1 2 3 0 1 -4 2 0 0 5 2 0 0 0 9 tril(A) ；下三角阵 ans = 1 0 0 0 5 1 0 0 3 0 5 0 11 15 0 9 rank(A) ；秩 ans = 4 norm(A) ；范数 ans = 21.3005 cond(A) ；条件数 ans = 11.1739 trace(A) ；迹 ans = 16 （2）略 6. A=[1 1 0.5;1 1 0.25;0.5 0.25 2] A = 1.0000 1.0000 0.5000 1.0000 1.0000 0.2500 0.5000 0.2500 2.0000

[V,D]=eig(A) V = 0.7212 0.4443 0.5315 -0.6863 0.5621 0.4615 -0.0937 -0.6976 0.7103 D = -0.0166 0 0 0 1.4801 0 0 0 2.5365

matlab程序设计与应用实验指导书

附件2 《matlab程序设计与应用》实验指导书 山东建筑大学信息与电气工程学院

前言 一、实验目的 本课程是电气工程及其自动化、自动化、电力工程与管理专业本科生的 学科基础选修课，它在线性代数、信号分析和处理、控制系统设计和仿真等 方面有着广泛的应用。主要是学习MATLAB的语法规则、基本命令和使用环境，使学生掌握MATLAB的基本命令和基本程序设计方法，提高使用该语言 的应用能力，具有使用MATLAB语言编程和调试的能力，以便为后续多门课 程使用该语言奠定必要的基础。 通过上机实验，使学生掌握MATLAB在线帮助功能的使用、熟悉MATLAB运行环境和MATLAB语言的主要特点，掌握MA TLAB语言的基本 语法规则及基本操作命令的使用，学会M文件的建立和使用方法以及应用MATLAB实现二维和三维图形的绘制方法，具有使用MATLAB语言编程和 调试的能力。 二、实验前预习 每次实验前，学生须仔细阅读本实验指导书的相关内容，明确实验目的、要求；明确实验步骤、测试数据及需观察的现象；复习与实验内容有关的理论知识；预习仪器设备的使用方法、操作规程及注意事项；做好预习要求中提出的其它事项。 三、实验注意事项 1．实验开始前，应先检查本组电脑设备是否工作正常，matlab软件工作是否正常。 2．实验时每位同学应独立完成实验任务，避免抄袭。 3．实验后应及时将实验数据进行记录与存盘，避免因电脑故障或其它原因造成实验数据的丢失。 4．实验中严格遵循电脑操作规程，如电脑发生异常现象，应立即切断电源，

报告指导教师检查处理。 5．测量数据或观察现象要认真细致，实事求是。 6．未经许可，不得对电脑进行软件及硬件的更改操作。 7．实验结束后，实验记录交指导教师查看并认为无误后，方可关掉电脑离开。 8．爱护公物，发生仪器设备等损坏事故时，应及时报告指导教师，按有关实验管理规定处理。 9．自觉遵守学校和实验室管理的其它有关规定。 四、实验总结 每次实验后，应对实验进行总结，即实验数据进行整理，绘制波形和图表，分析实验现象，撰写实验报告。实验报告除写明实验名称、日期、实验者姓名、同组实验者姓名外，还包括： 1．实验目的； 2．实验内容 3．实验流程图； 4．程序命令清单； 5．运行结果； 6．实验的收获与体会； 7．回答每项实验的有关问答题。

matlab程序设计实例

MATLAB 程序设计方法及若干程序实例 樊双喜 (河南大学数学与 信息科学学院开封475004) 摘要本文通过对 MATLAB 程序设计中的若干典型问题做简要的分析和总结,并在此基础上着重讨论了有关算法设计、程序的调试与测试、算法与程序的优化以及循环控制等方面的问题.还通过对一些程序实例做具体解析,来方便读者进行编程训练并掌握一些有关MATLAB 程序设计方面的基本概念、基本方法以及某些问题的处理技巧等.此外,在文章的最后还给出了几个常用数学方法的算法程序, 供读者参考使用.希望能对初学者进行 MATLAB 编程训练提供一些可供参考的材料,并起到一定的指导和激励作用,进而为MATLAB 编程入门打下好的基础. 关键字算法设计；程序调试与测试；程序优化；循环控制 1 算法与程序 1.1 算法与程序的关系算法被称为程序的灵魂,因此在介绍程序之前应先了 解什么是算法.所谓算 法就是对特定问题求解步骤的一种描述.对于一个较复杂的计算或是数据处理的问题,通常是先设计出在理论上可行的算法,即程序的操作步骤,然后再按照算法逐步翻译成相应的程序语言,即计算机可识别的语言. 所谓程序设计,就是使用在计算机上可执行的程序代码来有效的描述用于解决特定问题算法的过程.简单来说,程序就是指令的集合.结构化程序设计由于采用了模块分化与功能分解,自顶向下,即分而治之的方法,因而可将一个较复杂的问题分解为若干子问题,逐步求精.算法是操作的过程,而程序结构和程序流程则是算法的具体体现. 1.2MATLAB 语言的特点 MATLAB 语言简洁紧凑,使用方便灵活,库函数极其丰富,其语法规则与科技人员的思维和书写习惯相近,便于操作.MATLAB 程序书写形式自由,利用其丰富

D实验五 M文件和MATLAB程序设计

实验五M文件和MATLAB程序设计 一、实验目的 matlab作为一种高级计算机语言，不仅可以命令行方式完成操作，也具有数据结构、控制流、输入输出等能力，本次实验通过熟悉和掌握m文件的建立与使用方法，以及函数与控制程序流程语句的使用，使学生具备一定的编程和程序调试能力。 1．掌握M文件的使用方法。 2．掌握if语句和switch语句的使用 3. 掌握循环语句的使用 4. 通过练习理解MATLAB编程方法。 二、实验原理 1．m文件 用matlab语言编写的程序，称为m文件。M文件根据调用方式的不同分为两类，命令文件（Script file）和函数文件（Function file）。区别? 2．程序控制结构 1）顺序结构 2）选择结构 （1）if语句a) 单分支if语句b) 双分支if语句c) 多分支if语句 （2）switch 语句 （3）try语句 3）循环结构 （1）for 语句 （2）while语句 （3）break语句、continue语句、return使用,区别? 3．函数文件 function 输出形参表＝函数名（输入形参表） 注释说明部分 函数体语句 注意事项？ 三、实验要求 1．首先上机练习PPT中各种流程控制语句的有关实例。 2．然后上机练习下面的实验习题。 四、实验习题

1．数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2→1 3→10→5→16→8→4→2→1 6→3→10→5→16→8→4→2→1 运行下面的程序，按程序提示输入n=1,2,3,5,7，8,9等数来验证这一结论。 %classic "3n+1" problem from number theory. while 1 n=input('Enter n,negative quits:'); if n<=0 break end a=n; while n>1 if rem(n,2)==0 n=n/2; else n=3*n+1; end a=[a,n]; end a end Enter n,negative quits: 2． 编程求满足∑=>m i i 11000020的最小m 值。 a=0; i=1; while (a<100000) a=a+pow2(i); i=i+1; end m=i-1 3． 编写一个函数，计算下面函数的值，给出x 的值，调用该函数后，返回y 的值。 function [y]=myfun1(x)

《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程序设计实训

实习报告 实习名称Matlab程序设计实训 专业班级**** 姓名*** 学号*** 成绩评定 电气与信息工程学院 和谐勤奋求是创新

实习考核和成绩评定办法 1．实习成绩的考核由指导教师根据实习表现、实习报告、实习成果、现场操作、设计、口试或笔试等几个方面，给出各项权重，综合评定。该实习考核教研室主任审核，主管院长审批备案。 2．成绩评定采用五级分制，即优、良、中、及格、不及格。 3．参加本次实习时间不足三分之二或旷课四天以上者，不得参加本次考核，按不及格处理。 4．实习结束一周内，指导教师提交实习成绩和实习总结。 5．实习过程考核和实习成绩在教师手册中有记载。 实习报告内容 实习报告内容、格式各专业根据实习类别（技能实习、认识实习、生产实习、毕业实习等）统一规范，经教研室主任审核、主管院长审批备案。 注：1. 实习任务书和实习指导书在实习前发给学生，实习任务书放置在实习报告封面后和正文目录前。 2. 为了节省纸张，保护环境，便于保管实习报告，统一采用A4纸，实习报告建议双面打印（正文采 用宋体五号字）或手写，右侧装订。

2016/2017学年第二学期 《Matlab程序设计实训》计划 指导教师：****** 班级：自动化1541、2班 实习地点：1教8楼机房 一、实习目的 实习的目的：了解并掌握MATLAB软件的各种应用开发环境；了解并掌握MATLAB软件强大的科学运算功能；了解并掌握利用MATLAB语言程序设计流程进行程序设计的方法；了解并掌握MATLAB的图形处理功能；了解MATLAB软件的交互式仿真功能；了解并掌握MATLAB软件在控制系统建模、分析以及设计中的应用。 实习的任务：通过实习可以使学生将MATLAB软件与控制系统理论及仿真相结合，加深对控制系统理论知识的理解，验证理论知识的正确性，提高软件的实际操作能力，可以有效地培养学生分析及解决实际问题的能力，不断提高自学能力，养成良好的科学态度以及实事求是、严谨踏实的学习风气，成为真正的应用型本科人才。通过实训达到边学边练、理论与实际操作相结合的学习目的。 二、实习内容 掌握MATLAB编程环境、基础知识；掌握矩阵生成及运算、数组运算；矩阵函数；掌握常用绘图命令及方法；掌握基本图形处理、二维图形处理方法；掌握可视化M文件程序设计及调试；掌握Simulink仿真工具的建立和基本仿真方法。 了解并掌握基于MATLAB的控制系统数学模型及分析方法；了解并掌握基于MATLAB的控制系统PID控制器的基本设计方法。综合应用MATLAB程序设计语言，根据设计任务要求，完成软件程序的编制并调试运行。 四、实习要求 实习期间，严格遵守课堂纪律，不许迟到、早退；实习期间，学习态度端正，认真上机练习；熟练掌握MATLAB软件的应用方法；完成实习报告一份，内容充实，写出实习的体会与收获；实习报告打印成文。