matlab作业(7)

合集下载

matlab应用解析作业

MATLAB应用解析第2章1．创建double的变量，并进行计算。

(1) a=87，b=190，计算 a+b 、a-b 、a*b 。

(2) 创建 uint8 类型的变量，数值与(1)中相同，进行相同的计算。

(1)>> a=87a =87>> b=190b =190>> a+bans =277>> a-bans =-103>> a*bans =16530(2)>> c=uint8(87)c =87>> d=uint8(190)d =190>> c+dans =255>> c-dans =>> c*dans =2552．计算：(1) ()sin 60(2) e 3 (3) 3cos 4⎛⎫π ⎪⎝⎭（1）>> sind(60)ans =0.8660（2）>> exp(3)ans =20.0855（3）>> cos(3*pi/4)ans =-0.70713．设2u =，3v =，计算： (1) 4log uv v(2)()22e u v v u +-(3)(1) >> u=2;>> v=3;>> 4*u*v/log(v)ans =21.8457(2)>> (exp(u)+v)^2/(v^2-u)ans =15.4189(3)>> sqrt(u-3*v)/(u*v)ans =0 + 0.4410i4．计算如下表达式：(1) ()()3542i i -+(2) ()sin 28i -(1)>> (3-5*i)*(4+2*i)ans =22.0000 -14.0000i(2)>> sin(2-8*i)ans =1.3553e+03 + 6.2026e+02i5．判断下面语句的运算结果。

MATLAB练习作业(已做)2012

说明:如在测试时请将程序中的中文字符改为西文(部分程序为手写,并未经matlab验证)2012年春《MATLAB基础》第一次上机练习与作业练习１：脚本文件建立一个文件名为helloWorld的脚本文件，运行该文件时显示出如下文字：Hello World!I am going to learn MATLAB•提示: 用disp显示字符串，将要显示的字符号串用单引号引起来。

如'This is a string'练习２：变量获取并保存当前的日期与时间●用函数clock生成一个变量，变量名为start●用size查看start的维数，它是一个行向量还是一个列向量？●start包含什么内容？用help clock查看●用函数datestr将向量start转换成字符串，得到新的变量,名为startString●将start与startString保存为mat文件，文件名为startTime●在练习１建立的脚本文件helloWorld.m文件中，用load函数导入变量startTime, 并显示如下文字：I started learning MATLAB on * start date and time*练习3：标量你将要以指数增长的速度来学习MATLAB, 将如下内容添加到helloWorld.m文件中●假设你的学习时间是一个常量,为1.5 days，将此时间用秒表示，赋给变量tau●假设课程持续时间为5 days. 将这个时间单位转换为秒，保存在变量endofClass中？●将学到的知识描述为t的函数，函数方程为：●用函数datestr将向量start转换成字符串，得到新的变量startString●在课程结束时间endofClass，你将学到多少知识？用变量knowledgeAtEnd表示(指数函数exp)●用变量knowledgeAtEnd,显示如下语句：At the end of the MATLAB class, I will know X% of MATLAB提示: 将数转换成字符串，用num2str练习４：向量运算计算从课程开始到现在经过的时间，用秒表示.●在helloWorld.m中,建立一系列变量,变量名分别为:secPerMin, secPerHour,secPerDay,secPerMonth(假设每个月30.5天),以及secPerYear.●将变量按次序secPerYear,secPerMonth, secPerDay,secPerHour, secPerMin,1 排列,构成一个行向量,命名为secConversion？●用时钟函数clock生成一个向量currentTime●计算经过的时间elapsedTime,用currentTime与start相减.●通过向量secondConversion与elapsedTime的数量积运算计算时间t,●显示当前的时间水平,用变量currentKnowledge表示. ()●显示如下语句：At this time, I know X% of MATLAB练习5：向量函数计算你的学习轨迹.●在helloWorld.m中,建立线性时间向量tVec,包含从0到endOfClass的10000个样本点.●计算在每个时间点处对应的知识值,仍然用函数练习6：向量寻址什么时候你将学到50%的MATLAB知识?.●在向量knowledgeVec中,查找最接近0.5的元素所在的位置.●用halfTime保存对应的时间●显示如下语句:I will know half of MATLAB after X days要将halfTime用secPerDay转换成天数练习7：绘图画出学习的轨迹图?.●在helloWorld.m中,打开一个新的图形窗口(figure)●用向量tVec与knowledgeVec画出知识轨迹,画图时,将时间单位转换为天●用图形窗口中的zoom in 查看halfTime,与前面计算结果相比较.clear,clc;disp('Hello World!');disp('I am going to learn MA TLAB');start=clock;[startx,starty]=size(start);if startx>startydisp('一个行向量')elsedisp('一个列向量');endstartString=datestr(start);save startTime.mat start startString;waittime=input('为使程序达到较好的执行效果，请输入程序需要暂停时间（秒）：'); state = 0;h = waitbar(0, '请等待...', 'Name', '进度条', 'CreateCancelBtn', ...'state = 1; delete(h); clear h');h1 = findall(h, 'style', 'pushbutton');set(h1, 'string', '取消', 'fontsize', 10)for i = 1 : 100waitbar(i/100, h, ['已暂停时间百分比' num2str(i) '%']);pause(waittime/100);if statebreakendendif exist('h') == 1delete(h);endload startTimedisp('I started learning MATLAB on * start date and time*');tau=1.5*24*3600;endofClass=5*24*3600;knowledgeAtEnd=1-exp(-tau/endofClass);disp(['At the end of the MATLAB class, I will know ',num2str(knowledgeAtEnd*100),' % of MATLAB']);secPerMin=60;secPerHour=3600;secPerDay=24*secPerHour;secPerMonth=30.5*secPerDay;secPerYear=12*secPerMonth;secondConversion=[secPerYear,secPerMonth, secPerDay,secPerHour, secPerMin,1]; currentTime=clock;elapsedTime=currentTime-start;elapsedTime=elapsedTime.*secondConversion;currentknowledge=1-exp(-elapsedTime/endofClass);disp(['At this time, I know ',num2str(sum(currentknowledge)*100),'% of MATLAB']);tVec=linspace(0,endofClass,1000);knowledgeVec=1-exp(-tVec./endofClass);halfTimeindex=find(abs(knowledgeVec-0.5)<1e-3);halfTime=sum(tVec(halfTimeindex))/length(halfTimeindex)/secPerDay;disp(['I will know half of MATLAB after ',num2str(halfTime),' days']);plot(tVec/secPerDay,knowledgeVec);xlabel('Day');2012年春《MATLAB基础》第二次上机练习与作业练习１：函数文件建立一个文件名为plotSin的函数文件，要求：●函数声明为:functon plotSin(f1)●该函数用于画出频率为的正弦波在区间上的图像.●每个周期含有16个数据点function plotsin(f1)xnum=16*f1;x=linspace(0,2*pi,xnum);plot(x,sin(f1*x));练习２：条件语句●修改函数plotSin(f1)使之有两个输入变量,形式为plotSin(f1,f2)●如果输入的变量个数为1,执行以前写的绘图命令.否则显示文字'Two inputs were given'●提示:输入变量的数目用内置变量名:narginfunction plotsin(f1,f2)if nargin==1xnum=16*f1;x=linspace(0,2*pi,xnum);plot(x,sin(f1*x));elseerror('Two inputs were given');end练习3：绘图命令●修改PlotSin中的绘图命令,将数据点类型改为正方形,线形改为红色虚线,线宽为2.设置数据点填充色为黑色(属性名为:LineWidth, MarkerFaceColor)●如果有2个输入变量,新开一个绘图窗口,上下排列,并且激活上面图形窗口plotSin(6)plotSin(2,3)function plotsin(f1,f2)if nargin==1h1=figure(1);xnum=16*f1;x=linspace(0,2*pi,xnum);plot(x,sin(f1*x),'--rs','MarkerFaceColor','k','LineWidth',2);elseh2=figure(2);h1=figure(1);xnum=16*f1;x=linspace(0,2*pi,xnum);plot(x,sin(f1*x),'--rs','MarkerFaceColor','k','LineWidth',2);end练习４：半对数图在过去5年里,某课程学生人数分别为:15,25,55,115,144.班级规模似乎呈指数增长趋势,用半对数图来验证.要求绘制的图形效果如下●相关函数:semilogy,xlabel,ylabel,axis,title●数据点为玫红,正方形,没有折线相连.点的大小为10,点边线的宽度为4.x=1:5;y=[15,25,55,115,144];h=semilogy(x,y,' sr','MarkerSize',10,'Linewidth',4);xlabel('x');ylabel('y');title('My Semilogy plot');axis([0,6,10,300])set(h,'Color',[1 0.3 0.6]);练习5：条形图生成一个有5个随机数构成的向量,画出一幅条形图,条形块为红色.效果如下:bar(1:5,rand(5,1),'r')练习6：圆的绘制●建立一个画圆的函数文件.function [x y]=getCircle(center,r)运行程序得到的是圆上点的坐标●调用程序画出同心圆及五环图,效果如下col=['b','k','r','y','g'];subplot 121r1=1:5;center1=[0,0];for i=1:5[x,y]=getCircle(center1,r1(i));plot(x,y,col(i),'Linewidth',2);hold onendaxis equalaxis([-5 5 -5 5])subplot 122r2=0.5;center2=[-1 0;0 0;1 0;-0.5 -0.5;0.5 -0.5];for i=1:5[x,y]=getCircle(center2(i,:),r2);plot(x,y,col(i),'Linewidth',2);hold onendaxis equalaxis([-1.5 1.5 -1.5 1])-505-101-1.5-1-0.50.512012年春《MATLAB 基础》第三次上机练习与作业练习１：演示文件命令窗口中，键入：demo，出现MATLAB的demo演示窗口，然后选择：MATLAB下的Matrices，依次选取其中的Basic Matrix operation 、Matrix manipulation两个示例，进行MATLAB矩阵运算的学习。

matlab编程题整理

作业一1输出x,y两个中值较大的一个值x=input(‘x’);y=input(‘y’);if x>yxelseyend2输入x，计算y的值。

计算函数的值y=x+1,x<0,y=2x-1,x≧0x=input错误！未指定书签。

(‘x);if x<0y=x+1elsey=2*x-1end3输入一学生成绩，评定其等级，方法是：90~100分为“优秀”，80~89分为“良好”，70~79分为“中等”，60~69分为“及格”，60分为“不合格”x=input(‘x’)if x>100|x<0y=’输入错误’elseif x>=90y=’优秀’elseif x>=80y=’良好’elseif x>=70y=’中等’elseif x>=60y=’及格’elsey=’不合格’emd4某超市节日期间举办购物打折的促销活动，优惠办法是：每位顾客当天一次性购物在100元以上者，按九五折优惠；在200元以上者，按九折优惠；在300元以上者，按八五折优惠；在500元以上者，按八折优惠。

x=input(‘x’);if x>=500y=x*0.8elseif x>=300y=x*0.85elseif x>=200y=x*0.9elseif x>=100y=x*0.95else y=xend 5编程计算：s=1+2+3+…+100sum=0;for i=1:100sum=sum+i;endsum引申1!-2!+3!-4!+5!- (99)sum=0;for i=1:99pdr=1;for k=1:ipdr=pdr*k;endsum=sum+pdr*(-1)^(i-1);endsum引申1*2*3*4*……*100sum=1；for i=1:100sum=sum*iendsum6计算1~100的奇数和sum=0;for i=1:2:100sum=sum+i;endsum7百元买百鸡问题。

MATLAB作业

MATLAB作业⼀、必答题：1. MATLAB系统由那些部分组成？答：MATLAB系统主要由开发环境、MATLAB语⾔、MATLAB数学函数库、图形功能和应⽤程序接⼝五个部分组成。

2. 如何启动M⽂件编辑/调试器？答：在操作界⾯上选择“建⽴新⽂件”或“打开⽂件”操作时，M⽂件编辑/调试器将被启动。

在命令窗⼝中键⼊“edit”命令也可以启动M⽂件编辑/调试器。

3. 存储在⼯作空间中的数组能编辑吗？如何操作？答：存储在⼯作空间的数组可以通过数组编辑器进⾏编辑：在⼯作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输⼊修改内容即可。

4. 在MATLAB中有⼏种获得帮助的途径？答：在MATLAB中有多种获得帮助的途径：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器；（2）help命令：在命令窗⼝键⼊“help” 命令可以列出帮助主题，键⼊“help 函数名”可以得到指定函数的在线帮助信息；（3）lookfor命令：在命令窗⼝键⼊“lookfor 关键词”可以搜索出⼀系列与给定关键词相关的命令和函数（4）模糊查询：输⼊命令的前⼏个字母，然后按Tab键，就可以列出所有以这⼏个字母开始的命令和函数。

5. 有⼏种建⽴矩阵的⽅法？各有什么优点？答：（1）以直接列出元素的形式输⼊；（2）通过语句和函数产⽣；（3）.在m⽂件中创建矩阵；（4）从外部的数据⽂件中装⼊。

6. 命令⽂件与函数⽂件的主要区别是什么？答：命令⽂件: M⽂件中最简单的⼀种，不需输出输⼊参数，⽤M ⽂件可以控制⼯作空间的所有数据。

运⾏过程中产⽣的变量都是全局变量。

运⾏⼀个命令⽂件等价于从命令窗⼝中顺序运⾏⽂件⾥的命令，程序不需要预先定义，只要依次将命令编辑在命令⽂件中即可。

函数⽂件：如果M⽂件的第⼀个可执⾏⾏以function开始，便是函数⽂件，每⼀个函数⽂件定义⼀个函数。

基于simulink的Matlab仿真作业(电气工程专业)7

基于MATLAB/SIMULINK 直流电动机调速系统仿真建模张三（陕西西安西安科技大学 710054）摘要：论述了一种基于MA TLAB 语言的直流电机调速控制仿真系统，通过MA TLAB 语言中的SIMULINK 模块和模糊控制箱实现模糊控制仿真关键词：MA TLAB. SIMULINK,仿真系统。

0 引言MATLAB 是矩阵实验室（Matrix Laboratory ）的简称，是一种数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB 和Simulink 两大部分。

Simulink 提供一个动态系统建模、仿真和综合分析的集成环境。

在该环境中，无需大量书写程序，而只需要通过简单直观的鼠标操作，就可构造出复杂的系统。

Simulink 具有适应面广、结构和流程清晰及仿真精细、贴近实际、效率高、灵活等优点，并基于以上优点Simulink 已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。

同时有大量的第三方软件和硬件可应用于或被要求应用于Simulink 。

纵观实验研究发展历程，仿真起着重要的作用，其经历了物理模型仿真，模似计算机仿真和数字计算机仿真，仿真给我们研究系统带来巨大方便，但在进行程序设计时．尤其是在计算矩阵运算或画图时，利用FORTRAN 、c 语言编程调试工作效率较低，很不方便。

在1984年，MATLAB 的推出，为研究者打开了一个新局面。

在系统仿真应用中很快得到了发展。

一台直流并励电动机，铭牌额定参数为kW 17=N P ，V 220=N U ，A 9.88=NI ，r/min 3000=N n ，电枢回路总电阻Ω=087.0a R ，励磁回路总电阻Ω=5.181f R ，电动机转动惯量2m K 76.0⋅=g J 。

要求仿真该电动机的直接起动的过程。

1计算电动机参数励磁电流为 I=220/181.5=1.21A 励磁电感在恒定磁场控制时可取零点数电阻R a =0.087Ω，电枢电感估算L=0.0032H式中，p 为极对数；C 为计算系数，对于无补偿电机C=0.1，补偿电机C=0.4。

MATLAB与控制系统仿真大作业

>> G1=tf(num,den);
>> G1=zpk(G1)
G1 =
5 (s+0.6) ----------------(s+3) (s+2) (s+1)
>> num=[2 1]; den=[1 2.9 1]; G2=tf(num,den); >> G2=zpk(G2)
G2 =
2 (s+0.5)
>> xlabel('x') ylabel('exp(-2.*x)')
三、simulink 建模（每题 10 分，共 20 分）
1、已知单位负反馈的开环传递函数为 G(s) =
2 s2 + 4s
，试利用
simulink
建立系在单位阶跃输入作用下的模型。
要求答案包括：（1）simulink 建模结构图；
（2）在同一个坐标中的阶跃信号和响应曲线图。
2、已知系统的开环传递函数为
G1(s)
=
2s2 + 5s + 6 s2 + 2s + 3
、G2 (s)
=
s2
s+6 + 7s +1
，H
(s)
=
5(s + 2) s +10
求：建立 Simulink 仿真模型，并求出其系统在单位阶跃响应；
要求答案包括：（1）simulink 建模结构图；
x=
-2.9709
0.5491
3.6000
0.0509
2. 已知下列矩阵
2 3 1
−1 3 5
A
=

MATLAB实验练习题计算机南邮MATLAB数学实验大作业答案

“M A T L A B ”练习题要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。

1、求230x e x -=的所有根。

（先画图后求解）（要求贴图）>>solve('exp(x)-3*x^2',0)ans=-2*lambertw(-1/6*3^(1/2))-2*lambertw(-1,-1/6*3^(1/2))-2*lambertw(1/6*3^(1/2))2、求下列方程的根。

1)5510x x ++=a=solve('x^5+5*x+1',0);a=vpa(a,6)a=1.10447+1.05983*i-1.00450+1.06095*i-.-1.00450-1.06095*i1.10447-1.05983*i2）1sin 02x x -=至少三个根 >>fzero('x*sin(x)-1/2',3)ans=2.9726>>fzero('x*sin(x)-1/2',-3)ans=-2.9726>>fzero('x*sin(x)-1/2',0)ans=-0.74083）2sin cos 0x x x -=所有根>>fzero('sin(x)*cos(x)-x^2',0)ans=>>fzero('sin(x)*cos(x)-x^2',0.6)ans=0.70223、求解下列各题：1）30sin lim x x x x ->- >>symx;ans=1/62）(10)cos ,x y e x y =求 >>symx;>>diff(exp(x)*cos(x),10)ans=(-32)*exp(x)*sin(x)3）21/20(17x e dx ⎰精确到位有效数字）>>symx;>>vpa((int(exp(x^2),x,0,1/2)),17)ans=0.4）42254x dx x +⎰ >>symx;>>int(x^4/(25+x^2),x)ans=125*atan(x/5)-25*x+x^3/35）求由参数方程arctan x y t⎧⎪=⎨=⎪⎩所确定的函数的一阶导数dy dx 与二阶导数22d y dx 。

matlab大作业

Matlab语言及应用周一晚作业廖育州07010559%第一组程序%绘制选用潍柴发动机WD615.50时汽车行驶性能曲线图，Ft--Ua图clear，clc；%清除内存变量，清屏%n为发动机的转速n=[1000,1100,1200,1300,1400,1500,1600,1700,1800,1900,2000,2100,2200];%Te为发动机在对应以上转速n下的外特性输出转矩Te=[1100,1125,1140,1145,1150,1180,1140,1125,1090,1060,1010,975,940];%ig为变速器不同档下的速比ig=[11.40;7.94;5.63;4.06;2.81;1.96;1.39;1.00];%nt为传动效率，是变速器、辅助变速器、传动轴、万向节、主减速器等四者的效率的乘积nt=0.849;Rr=0.5145; %Rr是车轮实际滚动半径i0=4.625; %i0是所选用的主减速器的速比ig1=1./(i0.*ig); %ig1只是为计算方便而设立的中间值，无物理意义Ua=0.377*Rr.*ig1*n; %Ua是在不同档位下，发动机不同扭矩下的车速Ft=ig*Te*i0*nt/(Rr*1000); %Ft是在不同档位下，发动机不同扭矩下的驱动力for n=1:8 % 8个档位对应八条曲线，八次循环p1=polyfit(Ua(n,:),Ft(n,:),5); %把先前计算所得的同一档位下的速度和驱动力拟合成曲线x2=Ua(n,1):0.1:Ua(n,13); % 设定曲线的X轴的取值范围和精度y1=polyval(p1,x2); % 计算曲线的Y轴的取值plot(x2,y1,'r-'); % 用红线绘制曲线hold on; % 把持当前的曲线，准备绘制下一条曲线end % 循环到此结束set(gca,'ygrid','on'); % 绘制水平线title('选用潍柴发动机WD615.50时汽车行驶性能曲线，Ft--Ua图'); %设定图片的标题名ylabel('Ft/(kN)'),xlabel('Ua/(km/h)'); %设定图片的x、y轴名称%******************************************************************************* %第二组程序绘制一个900*900的黑白格图片，并输出clear，clc；%清除内存中的变量，清屏B=ones(150,150); %定义一个150*150的一矩阵，作为白格C=zeros(150,150); %定义一个150*150的零矩阵，作为黑格A=[B,C,B;C,B,C;B,C,B]; %定义黑白格位置关系imshow(A)%******************************************************************************* %第三组程序计算曲线拟合系数clear，clc; %清除内存中的变量，清屏%n是发动机转速n=[1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 2100 2200];n1=n/(1000); %n1无物理意义，只是作为方便计算的中间变量%ge为与转速n对应的油耗ge=[202 200 195 194 192 194 196 199 202 207 212 217 223 ];p1=polyfit(n1,ge,5); %用5阶拟合n1和gep1 %输出p1，p1为比油耗特性拟合多项式的系数向量%Te为发动机在不同转速下的转矩Te=[1100 1125 1140 1145 1150 1180 1140 1125 1090 1060 1010 975 940];%转矩p2=polyfit(n1,Te,5); %用5阶拟合n1和Tep2 %输出p2，p2为转矩特性拟合多项式的系数向量%******************************************************************************* %第四组程序求符号表达式的和与差syms x fx gx %定义符号变量等于符号表达式fx=2*x^2+3*x-5; %定义fxgx=x^2+7; %定义gxfx+gx %求解fx+gx的符号表达式%******************************************************************************* %第五组程序绘制三维立体峰图[x,y]=meshgrid(-3:1/10:3,-3:1/10:3); %定义x，y的取值范围及网格划分精度z=peaks(x,y); %计算z值surfc(x,y,z); %绘图%******************************************************************************* %第六组程序计算导数clc，clear %清屏，清除内存中的变量syms a b t x y z; %定义a b t x y z等为变量符合f=sqrt(1+exp(x)); %定义原始式子fA=diff(f) %对f求一阶导数B=diff(f,x,3) %对f求三阶导数%****************************************************************************** %第七组程序求方程的所有解clear，clc；%清除内存中的变量，清屏p=[5,10,3,-71,9,13]; %定义方程x=roots(p) %求解方程的根%******************************************************************************* %第八组程序%找出一串奇数个数字串中的中间那个数字，或者偶数个数字串中的中间那两个数字的平均值clc , clear; %清屏,清除内存中的变量,A=input('请输入一串数据，并以回车键结束输入','s'); % 输入一串数字字符n=length(A); %计算字符串A的长度i=mod(n,2); %求余if i==1 %如果奇数个字符m=(n-1)/2+1; %找出中间那个数字middle=A(m) %输出中间那个数字else %如果是偶数个字符m=n/2;z=A(m)+A(m+1)-96; %把字符的ASCII码转换成数字s=z/2 %求出中间那两个数字的平均值end %程序到此结束。

MATLAB平时作业(图文版)

MATLAB 平时作业第一章习题16. 以下两种说法对吗？（1）“MATLAB 的数值表达精度与其指令窗中的数据显示精度相同。

”答：此种说法错误。

MATLAB 提供了控制数据显示格式的控制指format ，该指令并不改变MATLAB 内存中变量的精度，只是改变其显示精度。

（2）“MATLAB 指令窗中显示的数值有效位数不超过7位。

”答：此种说法错误。

当变量小于1000时，使用format 或format short 后，或者默认情况下，变量的显示精度最多不超过7位，但显示精度不等于变量的精度。

7. 想要在MATLAB 中产生二维数组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321S ，下面哪些指令能实现目的？ S=[1,2,3;4,5,6;7,8;9]S=[1 2 3;4 5 6;7 8 9]S=[1，2，3；4，5，6；7，8，9] %整个指令在中文状态下输入答：操作如图：第1、2条指令可以实现，第3条指令不可实现。

第三章习题31．在MATLAB 中，先运行指令A=magic(3), B=[1,2,1;3,4,3;5,6,7], C=reshape(1:6,3,2)生成阵列33⨯A ，23⨯B ，23⨯C ，然后根据运行结果回答以下问题：运行结果如图：（1）计算A*B, B*A ，这两个乘积相同吗？计算结果如图：答：不同。

（2）计算A\B, B/A ，左除、右除结果相同吗？计算结果如图：答：不同。

（3）计算B( : ,[1,2]).*C和C.*B( : , [1,2])，这两个乘积相同吗？计算结果如图答：相同。

（4）计算A\A和A.\A，这两个计算结果相同吗？计算结果如图：答：相同。

（5）计算A\eye(3)和inv(A)，这两个计算结果相同吗？计算结果如图：答：不同。

（提示：根据对计算结果的目测回答问题）2．在MATLAB中，先运行A=[1, 2; 3, 4],b=0.5,C=[4, 2; 1, 0.5], 然后根据计算结果回答以下问题：创建数据步骤略（1）计算A^b和A.^b, 这两个计算结果相同吗？答：不同。

MATLAB作业

>>demoThe element type"name"must be terminated by the matching end-tag"</name>".org.xml.sax.SAXParseException:The element type"name"must be terminated by the matching end-tag"</name>".at org.apache.xerces.parsers.DOMParser.parse(DOMParser.java:235)at org.apache.xerces.jaxp.DocumentBuilderImpl.parse(DocumentBuilderImpl.java:201)at com.mathworks.mde.help.Demos.createNodes(Demos.java:416)at com.mathworks.mde.help.Demos.access$900(Demos.java:46)at com.mathworks.mde.help.Demos$LoadDemosThread.run(Demos.java:276)at com.mathworks.jmi.NativeMatlab.dispatchMTRequests(NativeMatlab.java:250)Could not parse G:\toolbox\ccslink\ccslink\demos.xmldemo：演示作用>>A=magic(3)A=816357492>>[sum(A),sum(A'),trace(A),trace(rot90(A))]ans=1515151515151515>>t=[0:0.1:2*pi];y=t.*sin(2.*t);plot(t,y)>>[x,y]=meshgrid(-5:0.1:5);>>z=(x.^2+y).*exp(-x.^2-(y.^2)/2); >>surf(x,y,z),colorbar>>whoYour variables are:A ans t x y z>>whosName Size Bytes ClassA3x372double array ans1x864double array t1x63504double array x101x10181608double array y101x10181608double array z101x10181608double array Grand total is30683elements using245464byteswho和whos：显示当前工作空间的变量情况>>which first'first'not foundwhich：显示出某个MATLAB文件所在目录>>workspaceworkspace：打开工作空间的浏览器窗口quit和exit：退出程序>>pathMATLABPATHG:\toolbox\matlab\generalG:\toolbox\matlab\opsG:\toolbox\matlab\langG:\toolbox\matlab\elmatG:\toolbox\matlab\elfunG:\toolbox\matlab\specfunG:\toolbox\matlab\datafunG:\toolbox\matlab\polyfunG:\toolbox\matlab\funfunG:\toolbox\matlab\sparfunG:\toolbox\matlab\scribeG:\toolbox\matlab\graph2dG:\toolbox\matlab\graph3dG:\toolbox\matlab\specgraphG:\toolbox\matlab\graphicsG:\toolbox\matlab\uitoolsG:\toolbox\matlab\strfunG:\toolbox\matlab\imagesciG:\toolbox\matlab\iofunG:\toolbox\matlab\audiovideoG:\toolbox\matlab\timefunG:\toolbox\matlab\datatypesG:\toolbox\matlab\verctrlG:\toolbox\matlab\codetoolsG:\toolbox\matlab\helptoolsG:\toolbox\matlab\winfunG:\toolbox\matlab\demosG:\toolbox\localG:\toolbox\matlabxl\matlabxlG:\toolbox\simulink\simulinkG:\toolbox\simulink\blocksG:\toolbox\simulink\componentsG:\toolbox\simulink\fixedandfloatG:\toolbox\simulink\fixedandfloat\fxpdemos G:\toolbox\simulink\fixedandfloat\obsoleteG:\toolbox\simulink\simdemosG:\toolbox\simulink\simdemos\aerospaceG:\toolbox\simulink\simdemos\automotiveG:\toolbox\simulink\simdemos\simfeaturesG:\toolbox\simulink\simdemos\simfeatures\mdlref G:\toolbox\simulink\simdemos\simgeneralG:\toolbox\simulink\simdemos\simnewG:\toolbox\simulink\deeG:\toolbox\shared\dastudioG:\toolbox\stateflow\stateflowG:\toolbox\rtw\rtwG:\toolbox\shared\hdsG:\toolbox\shared\timeseriesG:\toolbox\stateflow\sfdemosG:\toolbox\rtw\rtwdemosG:\toolbox\rtw\rtwdemos\rsimdemosG:\toolbox\rtw\targets\asap2\asap2G:\toolbox\rtw\targets\asap2\asap2\userG:\toolbox\rtw\targets\rtwin\rtwinG:\toolbox\simulink\acceleratorG:\toolbox\rtw\accelG:\toolbox\aeroblks\aeroblksG:\toolbox\aeroblks\aerodemosG:\toolbox\aeroblks\aerodemos\textureG:\toolbox\bioinfo\bioinfoG:\toolbox\bioinfo\microarrayG:\toolbox\bioinfo\proteinsG:\toolbox\bioinfo\biomatricesG:\toolbox\bioinfo\biodemosG:\toolbox\rtw\targets\common\can\blocksG:\toolbox\rtw\targets\common\configuration\resource G:\toolbox\rtw\targets\common\tgtcommonG:\toolbox\rtw\targets\c166\c166G:\toolbox\rtw\targets\c166\blocksG:\toolbox\rtw\targets\c166\c166demosG:\toolbox\ccslink\ccslinkG:\toolbox\ccslink\ccsblksG:\toolbox\ccslink\ccsdemosG:\toolbox\cdma\cdmaG:\toolbox\cdma\cdmamasksG:\toolbox\cdma\cdmamexG:\toolbox\cdma\cdmademosG:\toolbox\combuilder\combuilderG:\toolbox\comm\commG:\toolbox\comm\commdemosG:\toolbox\comm\commdemos\commdocdemosG:\toolbox\comm\commobsoleteG:\toolbox\commblks\commblksG:\toolbox\commblks\commmasksG:\toolbox\commblks\commmexG:\toolbox\commblks\commblksdemosG:\toolbox\commblks\commblksobsolete\v2p5G:\toolbox\commblks\commblksobsolete\v2G:\toolbox\commblks\commblksobsolete\v1p5G:\toolbox\compilerG:\toolbox\control\controlG:\toolbox\control\ctrlguisG:\toolbox\control\ctrlobsoleteG:\toolbox\control\ctrlutilG:\toolbox\control\ctrldemosG:\toolbox\shared\controllibG:\toolbox\curvefit\curvefitG:\toolbox\curvefit\cftoolguiG:\toolbox\shared\optimlibG:\toolbox\daq\daqG:\toolbox\daq\daqguisG:\toolbox\daq\daqdemosG:\toolbox\database\databaseG:\toolbox\database\dbdemosG:\toolbox\database\vqbG:\toolbox\datafeed\datafeedG:\toolbox\datafeed\dfguiG:\toolbox\dialsG:\toolbox\dspblks\dspblksG:\toolbox\dspblks\dspmasksG:\toolbox\dspblks\dspmexG:\toolbox\dspblks\dspdemosG:\toolbox\rtw\targets\ecoderG:\toolbox\rtw\targets\ecoder\ecoderdemos G:\toolbox\rtw\targets\mptG:\toolbox\rtw\targets\mpt\mptG:\toolbox\rtw\targets\mpt\user_specific G:\toolbox\exlinkG:\toolbox\symbolic\extendedG:\toolbox\filterdesign\filterdesignG:\toolbox\filterdesign\quantizationG:\toolbox\filterdesign\filtdesdemosG:\toolbox\finance\financeG:\toolbox\finance\calendarG:\toolbox\finance\findemosG:\toolbox\finance\finsupportG:\toolbox\finderiv\finderivG:\toolbox\finfixed\finfixedG:\toolbox\fixedpoint\fixedpointG:\toolbox\fixedpoint\fidemosG:\toolbox\fixedpoint\fimexG:\toolbox\fixpointG:\toolbox\ftseries\ftseriesG:\toolbox\ftseries\ftsdemosG:\toolbox\ftseries\ftsdataG:\toolbox\ftseries\ftstutorialsG:\toolbox\fuzzy\fuzzyG:\toolbox\fuzzy\fuzdemosG:\toolbox\gadsG:\toolbox\gads\gadsG:\toolbox\gads\gadsdemosG:\toolbox\garch\garchG:\toolbox\garch\garchdemosG:\toolbox\rtw\targets\hc12\hc12G:\toolbox\rtw\targets\hc12\blocks G:\toolbox\rtw\targets\hc12\codewarrior G:\toolbox\rtw\targets\hc12\hc12demos G:\toolbox\hdlfilter\hdlfilterG:\toolbox\hdlfilter\hdlfiltdemosG:\toolbox\shared\hdlsharedG:\toolbox\ident\identG:\toolbox\ident\idobsoleteG:\toolbox\ident\idguisG:\toolbox\ident\idutilsG:\toolbox\ident\iddemosG:\toolbox\ident\idhelpG:\toolbox\images\imagesG:\toolbox\images\imdemosG:\toolbox\imaq\imaqG:\toolbox\imaq\imaqdemosG:\toolbox\instrument\instrumentG:\toolbox\instrument\instrumentdemos G:\toolbox\lmi\lmictrlG:\toolbox\lmi\lmilabG:\toolbox\map\mapG:\toolbox\map\mapdemosG:\toolbox\map\mapdispG:\toolbox\map\mapformatsG:\toolbox\map\mapprojG:\toolbox\mbc\mbcG:\toolbox\mbc\mbcdataG:\toolbox\mbc\mbcdesignG:\toolbox\mbc\mbcexprG:\toolbox\mbc\mbcguitoolsG:\toolbox\mbc\mbclayoutsG:\toolbox\mbc\mbcmodelsG:\toolbox\mbc\mbcsimulinkG:\toolbox\mbc\mbctoolsG:\toolbox\mbc\mbcviewG:\toolbox\physmod\mech\mechG:\toolbox\physmod\mech\mechdemosG:\toolbox\physmod\pmimport\pmimportG:\toolbox\modelsim\modelsimG:\toolbox\modelsim\modelsimdemosG:\toolbox\mpc\mpcG:\toolbox\mpc\mpcdemosG:\toolbox\mpc\mpcguisG:\toolbox\mpc\mpcobsoleteG:\toolbox\mpc\mpcutilsG:\toolbox\shared\slcontrollibG:\toolbox\rtw\targets\mpc555dkG:\toolbox\rtw\targets\mpc555dk\common\configuration G:\toolbox\rtw\targets\mpc555dk\mpc555demosG:\toolbox\rtw\targets\mpc555dk\mpc555dkG:\toolbox\rtw\targets\mpc555dk\pilG:\toolbox\rtw\targets\mpc555dk\rt\blockset\mfiles G:\toolbox\rtw\targets\mpc555dk\rt\blocksetG:\toolbox\mutools\commandsG:\toolbox\mutools\subsG:\toolbox\nnet\nnetG:\toolbox\nnet\nnutilsG:\toolbox\nnet\nncontrolG:\toolbox\nnet\nndemosG:\toolbox\nnet\nnobsoleteG:\toolbox\opc\opcG:\toolbox\opc\opcguiG:\toolbox\opc\opcdemosG:\toolbox\optimG:\toolbox\rtw\targets\osek\osekG:\toolbox\rtw\targets\osek\osekdemosG:\toolbox\rtw\targets\osek\blocksG:\toolbox\rtw\targets\osek\osekworksG:\toolbox\rtw\targets\osek\proosekG:\toolbox\pdeG:\toolbox\physmod\powersys\powersysG:\toolbox\physmod\powersys\powerdemoG:\toolbox\simulink\reqmgtG:\toolbox\simulink\reqmgt\rmidemosG:\toolbox\rf\rfG:\toolbox\rf\rfdemosG:\toolbox\rf\rftoolG:\toolbox\rfblks\rfblksG:\toolbox\rfblks\rfblksmasksG:\toolbox\rfblks\rfblksmexG:\toolbox\rfblks\rfblksdemosG:\toolbox\robustG:\toolbox\rptgen\rptgenG:\toolbox\rptgen\rptgendemosG:\toolbox\rptgenext\rptgenextG:\toolbox\rptgenext\rptgenextdemosG:\toolbox\signal\signalG:\toolbox\signal\sigtoolsG:\toolbox\signal\sptoolguiG:\toolbox\signal\sigdemosG:\toolbox\slcontrol\slcontrolG:\toolbox\slcontrol\slctrlguisG:\toolbox\slcontrol\slctrlutilG:\toolbox\slcontrol\slctrldemosG:\toolbox\slestim\slestimG:\toolbox\slestim\slestmexG:\toolbox\slestim\slestdemosG:\toolbox\slestim\slestguisG:\toolbox\sloptim\sloptimG:\toolbox\sloptim\sloptguisG:\toolbox\sloptim\sloptdemosG:\toolbox\sloptim\sloptobsoleteG:\toolbox\simulink\slvnvG:\toolbox\simulink\simcoverageG:\toolbox\simulink\simcoverage\simcovdemos G:\toolbox\splinesG:\toolbox\statsG:\toolbox\symbolicG:\toolbox\rtw\targets\tic2000\tic2000G:\toolbox\rtw\targets\tic2000\tic2000blksG:\toolbox\rtw\targets\tic2000\tic2000demos G:\toolbox\shared\etargets\etargetsG:\toolbox\rtw\targets\tic6000\tic6000G:\toolbox\rtw\targets\tic6000\tic6000blksG:\toolbox\rtw\targets\tic6000\tic6000demos G:\toolbox\shared\etargets\rtdxblksG:\toolbox\vr\vrG:\toolbox\vr\vrdemosG:\toolbox\wavelet\waveletG:\toolbox\wavelet\wavedemoG:\toolbox\webserver\webserverG:\toolbox\webserver\wsdemosG:\toolbox\rtw\targets\xpc\xpcG:\toolbox\rtw\targets\xpc\target\build\xpcblocksG:\toolbox\rtw\targets\xpc\xpcdemosG:\toolbox\rtw\targets\xpc\xpc\xpcmngrG:\toolbox\rtw\targets\xpc\target\kernel\embeddedG:\workpath：显示MATLAB包含的路径>>cdG:\workcd：显示当前路径>>addpath a>In path at115In addpath at95addpath：将一个新的字符串表示的子目录名添加到MATLAB包含的路径下rmpath：从当前的MATLAB包含的某个路径下清除某个字符串表示的子目录>>format long>>pians=3.14159265358979format long：小数点后精确到7位（单精度）或14位（双精度）有效数字>>format short e>>pians=3.1416e+000format short e：5位科学记数表示>>format long e>>pians=3.141592653589793e+000format long e：16位科学记数表示>>format short g>>pians=3.1416format short g：从format short和format short e中自动选择最佳计数方式>>format long g>>pians=3.14159265358979format long g：从format long和format long e中自动选择最佳计数方式>>format hex>>pians=400921fb54442d18format hex：十六进制表示>>format rat>>pians=355/113format rat：近似有理数表示clc：清除命令窗口cla：清除当前坐标系下的内容clf：清除图形窗口close：关闭所有的命令窗口copyfile：复制文件>>pwdans=G:\workpwd：显示当前的工作目录名称>>mkdir a>>dir...amkdir：建立目录dir：获取当前目录并保存在列表中>>delete aWarning:'G:\work\a'is a e rmdir to delete directories. delete：删除文件。

matlab综合大作业(附详细答案)

m a t l a b综合大作业(附详细答案)-标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII《MATLAB语言及应用》期末大作业报告1．数组的创建和访问（20分，每小题2分）：1)利用randn函数生成均值为1，方差为4的5*5矩阵A；实验程序：A=1+sqrt(4)*randn(5)实验结果：A =0.1349 3.3818 0.6266 1.2279 1.5888-2.3312 3.3783 2.4516 3.1335 -1.67241.2507 0.9247 -0.1766 1.11862.42861.5754 1.6546 5.3664 0.8087 4.2471-1.2929 1.3493 0.7272 -0.6647 -0.38362)将矩阵A按列拉长得到矩阵B；实验程序：B=A(:)实验结果：B =0.1349-2.33121.25071.5754-1.29293.38183.37830.92471.65461.34930.62662.4516-0.17665.36640.72721.22793.13351.11860.8087-0.66471.5888-1.67242.42864.2471-0.38363)提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的矩阵C；实验程序：C=[A(2,2),A(2,4);A(3,2),A(3,4)]实验结果：C =3.3783 3.13350.9247 1.11864)寻找矩阵A中大于0的元素；]实验程序：G=A(find(A>0))实验结果：G =0.13491.25071.57543.38183.37830.92471.65461.34930.62662.45165.36640.72721.22793.13351.11860.80871.58882.42864.24715)求矩阵A的转置矩阵D；实验程序：D=A'实验结果：D =0.1349 -2.3312 1.2507 1.5754 -1.29293.3818 3.3783 0.9247 1.6546 1.34930.6266 2.4516 -0.1766 5.3664 0.72721.2279 3.1335 1.1186 0.8087 -0.66471.5888 -1.67242.4286 4.2471 -0.38366)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E；实验程序：E=flipud(fliplr(A))实验结果：E =-0.3836 -0.6647 0.7272 1.3493 -1.29294.2471 0.80875.3664 1.6546 1.57542.4286 1.1186 -0.1766 0.9247 1.2507-1.6724 3.1335 2.4516 3.3783 -2.33121.5888 1.2279 0.6266 3.3818 0.13497)删除矩阵A的第2列和第4列得到矩阵F；实验程序：F=A;F(:,[2,4])=[]实验结果：F =0.1349 0.6266 1.5888-2.3312 2.4516 -1.67241.2507 -0.17662.42861.5754 5.3664 4.2471-1.2929 0.7272 -0.38368)求矩阵A的特征值和特征向量；实验程序：[Av,Ad]=eig(A)实验结果：特征向量Av =-0.4777 0.1090 + 0.3829i 0.1090 - 0.3829i -0.7900 -0.2579 -0.5651 -0.5944 -0.5944 -0.3439 -0.1272-0.2862 0.2779 + 0.0196i 0.2779 - 0.0196i -0.0612 -0.5682 -0.6087 0.5042 - 0.2283i 0.5042 + 0.2283i 0.0343 0.6786 0.0080 -0.1028 + 0.3059i -0.1028 - 0.3059i 0.5026 0.3660 特征值Ad =6.0481 0 0 0 00 -0.2877 + 3.4850i 0 0 00 0 -0.2877 - 3.4850i 0 00 0 0 0.5915 00 0 0 0 -2.30249)求矩阵A的每一列的和值；实验程序：lieSUM=sum(A)实验结果：lieSUM =-0.6632 10.6888 8.9951 5.6240 6.208710)求矩阵A的每一列的平均值；实验程序：average=mean(A)实验结果：average =-0.1326 2.1378 1.7990 1.1248 1.24172．符号计算（10分，每小题5分）：1)求方程组20,0++=++=关于,y z的解；uy vz w y z w实验程序：S = solve('u*y^2 + v*z+w=0', 'y+z+w=0','y,z');y= S. y, z=S. z实验结果：y =[ -1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))-w] [ -1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))-w] z =[ 1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))] [ 1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))]2)利用dsolve 求解偏微分方程,dx dyy x dt dt==-的解；实验程序：[x,y]=dsolve('Dx=y','Dy=-x')实验结果：x =-C1*cos(t)+C2*sin(t)y = C1*sin(t)+C2*cos(t)3．数据和函数的可视化（20分，每小题5分）：1)二维图形绘制：绘制方程2222125x y a a +=-表示的一组椭圆，其中0.5:0.5:4.5a =；实验程序：t=0:0.01*pi:2*pi; for a=0.5:0.5:4.5; x=a*cos(t); y=sqrt(25-a^2)*sin(t); plot(x,y) hold on end实验结果：2) 利用plotyy 指令在同一张图上绘制sin y x =和10x y =在[0,4]x ∈上的曲线；实验程序：x=0:0.1:4; y1=sin(x); y2=10.^x;[ax,h1,h2]=plotyy(x,y1,x,y2); set(h1,'LineStyle','.','color','r'); set(h2,'LineStyle','-','color','g'); legend([h1,h2],{'y=sinx';'y=10^x'});实验结果：3)用曲面图表示函数22z x y =+；实验程序：x=-3:0.1:3; y=-3:0.1:3; [X,Y]=meshgrid(x,y); Z=X.^2+Y.^2; surf(X,Y,Z)实验结果：4)用stem 函数绘制对函数cos 4y t π=的采样序列；实验程序：t=-8:0.1:8;y=cos(pi.*t/4); stem(y)实验结果：4. 设采样频率为Fs = 1000 Hz ，已知原始信号为)150π2sin(2)80π2sin(t t x ⨯+⨯=，由于某一原因，原始信号被白噪声污染，实际获得的信号为))((ˆt size randn x x+=，要求设计出一个FIR 滤波器恢复出原始信号。

MATLAB大作业

M A T L A B大作业作业要求：（1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。

（2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。

打印稿不要求MATLAB程序，但电子稿要包含MATLAB程序。

（3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。

90问题五：利用MATLAB软件绘制一朵鲜花，实现一定的仿真效果。

提示：二维/三维绘图，对花瓣、花蕊、叶片、花杆等的形状和颜色进行详细设置。

第二类：插值与拟合。

（B级）问题一：有人对汽车进行了一次实验，具体过程是，在行驶过程中先加速，然后再保持匀速行驶一段时间，接着再加速，然后再保持匀速，如此交替。

注意，整个实验过程中从未减速。

在一组时间段50个时间点的速度。

（2）绘制插值图形并标注样本点。

问题二：估算矩形平板各个位置的温度。

已知平板长为5m，宽为3m，平板上3×5栅格点上的温度值为44，25，20，24，30；42，21，20，23，38；25，23，19，27，40。

（1）分别使用最近点插值、线性插值和三次样条插值进行计算。

（2）用杆图标注样本点。

（3）绘制平板温度分布图。

对a,b,c,d的值。

提示：曲线拟合并绘图分析第三类：定积分问题。

（B级）问题一：地球密度随着离中心（r=0）距离的变化而变化，不同半径处的密度如表所示，试估问题二：河道平均流量Q（m3/s）可使用速度和深度的乘积的积分来计算（河道横截面不规则），公式如下。

其中V(x)是离岸x（m）距离处的水速（m/s）,H(x)是离岸x距离处的水深（m）。

根据收集到过5（1（2（3（Q，单位是m（1（2（1（2（3）将节点1的力改为方向向上，计算这种改变对H2和V2的影响。

（4）将节点1的力撤销，而在节点1和2处施加1500N的水平外力，求节点3处垂直反作用力（V3）。

MATLAB第一章作业答案

第一章M A T L A B概况与基本操作1.选择题：(1)最初的MATLAB核心程序是采用A语言编写的。

(2)即将于2011年9月发布的MATLAB新版本的编号为D。

2011Ra 2011Rb R2011a R2011b(3)在默认设置中，MATLAB中的注释语句显示的颜色是D。

A.黑色B.蓝色C.红色D.绿色(4)如果要以科学计数法显示15位有效数字，使用的命令是B。

long long e long g long d(5)在命令窗口新建变量a、b，如果只查看变量a的详细信息，使用的命令为B。

a a(6)如果要清除工作空间的所有变量，使用的命令为 C 。

all C.两者都可 D.两者都不可(7)在创建变量时，如果不想立即在命令窗口中输出结果，可以在命令后加上D。

A.冒号B.逗号C.空格D.分号(8)如果要重新执行以前输入的命令，可以使用B键。

A.下箭头↓B.上箭头↑C.左箭头←D.右箭头→(9)如果要查询函数inv的相关信息，并显示在命令窗口，应使用命令A。

D.三者均可(10)如果要启动Notebook文档，下列D操作是可行的。

A.在命令窗口输入notebook命令B.在命令窗口输入notebook filename命令C.在Word中启动M-book文档D.三者均可2.填空题：(1)MATLAB是MATrix和LABoratory两个单词前三个字母的组合，意为“矩阵实验室”，它的创始人是Cleve Moler和Jack Little。

(2)在MATLAB的默认设置中，关键字显示的字体为蓝色，命令、表达式、计算结果显示的字体为黑色，字符串显示的字体为紫色，注释显示的字体为绿色，错误信息显示的字体为红色。

(3)在命令窗口中，输出结果显示为各行之间添加空行的命令为format compact，各行之间不添加空行的命令为format compact。

备注：本题布置给大家时有一点小错误，现在予以更正。

MATLAB基础训练作业(含答案)

实验一 MATLAB 工作环境熟悉及简单命令的执行一、实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。

二、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单（Menu ）、工具条（Toolbar ）、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。

三、实验步骤：1、启动MATLAB ，熟悉MATLAB 的桌面。

2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。

（1）（365-522-70）3（2）>>area=pi*^2（3）已知x=3，y=4，在MATLAB 中求z ：()232y x y x z -= （4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。

m1=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡11514412679810115133216 执行以下命令>>m1( 2 , 3 )>>m1( 11 )>>m1( : , 3 )>>m1( 2 : 3 , 1 : 3 )>>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1)（5）执行命令>>help abs查看函数abs 的用法及用途，计算abs( 3 + 4i )（6）执行命令>>x=0::6*pi;>>y=5*sin(x);>>plot(x,y)（6）运行MATLAB 的演示程序，>>demo ，以便对MATLAB 有一个总体了解。

四、思考题1、以下变量名是否合法为什么（1）x2（2）3col（3）_row（4）for2、求以下变量的值，并在MATLAB 中验证。

MATLAB第三章和第四章作业

习题三1.选择题(1)已知a=4，b='4'，下面说法错误的为C 。

A. 变量a比b占用的存储空间大B. 变量a和b可以进行加、减、乘、除运算C. 变量a和b的数据类型相同D. 变量b可以用eval命令执行(2)已知s='显示"hello"'，则s的元素个数是B 。

A. 12B. 9C. 7D. 18(3)运行命令>>a=sym('pi','d')则对于变量a的描述D是正确的。

A. a是符号变量B. a显示为10位的数值C. a显示为32位的数值D. a不存在(4)对于符号表达式g=sym('sin(a*z)+cos(w*v)')，如果运行命令diff(g)，则求导自变量是D。

A. aB. zC. wD. v(5)运行命令>>a=double(sym('sin(pi/2)'))，则变量a是 C。

A. 符号变量B. 字符串'1'C. double型的数值1D. 出错(6)运行命令>>y=dsolve('x*D2y-3*Dy=x^2','t')求解微分方程，则D。

A. Dy是指dy/dxB. 得出的y是通解有一个常数符C1C. D2y是指d2y/dxD. 得出的y是通解有两个常数符C1和C2(7)有一个2行2列的元胞数组c，则c(2)是指C 。

A. 第1行第2列的元素内容B. 第2行第1列的元素内容C. 第1行第2列的元素D. 第2行第1列的元素(8)对于一个2行2列的元胞数组c，如果要删除第二列的元素，应该使用的命令是B 。

A.c{:,2}=[]B.c(:,2)=[]C.A、B两种方式均可D.A、B两种方式均不可(9)有一个3行10列的构架数组student，每个构架有name（姓名）、scores两个字段，其中scores是以1×5数组表示的五门课的成绩，那么要查看第4个学生的第2门课成绩，采用的命令中达不到要求的是。

MATLAB大作业

M A T L A B大作业(总15页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--MATLAB大作业作业要求：（1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。

（2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。

打印稿不要求MATLAB程序，但电子稿要包含MATLAB程序。

（3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。

不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。

（4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。

（5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。

第一类：绘制图形。

（B级）问题一：斐波那契（Fibonacci）螺旋线，也称黄金螺旋线（Golden spiral），是根据斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典黄金比例。

斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契螺旋线，如图所示。

问题二：绘制谢尔宾斯基三角形（Sierpinskitriangle）是一种分形，由波兰数学家谢尔宾斯基在1915年提出,它是一种典型的自相似集。

其生成过程为：取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三角形，然后去掉中间的那一个小三角形。

接下来对其余三个小三角形重复上述操作，如图所示。

问题三：其他分形曲线或图形。

分形曲线还有很多，教材介绍了科赫曲线，其他还有皮亚诺曲线、分形树、康托（G. Cantor）三分集、Julia集、曼德布罗集合（Mandelbrot set），等等。

matlab作业练习及答案(有程序和截图)

1. 给出如下方程组：16323313915111081697612534141543s u p w s u p w s u p w s u p w +++=+++=+++=+++=求s ，u ，p ，w 的值，并求出系数行列式。

[答案：s ＝-0.1258，u ＝-8.7133，p ＝11.2875，w ＝-0.0500。

行列式=7680。

]>> A=[16 32 33 13;5 11 10 8;9 7 6 12;34 14 15 1];>> B=[91 16 5 43]';>> C=A\B2. 求矩阵H1()H X X X X -''=其中1731565419289121110X ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦〔答案：H ＝〔0.7294，0.9041，0.4477，0.9188〕'〕>> X=[17 31 5;6 5 4;19 28 9;12 11 10]; >> H=X/(X'*X)*X'>> X=[17 31 5;6 5 4;19 28 9;12 11 10]; >> H=X*inv((X'*X))*X'3．用MA TLAB语言实现下面的分段函数,()/,||,h x Dy f x h Dx x Dh x D>⎧⎪==≤⎨⎪-<-⎩。

y=h*(x>D)+h/D*x.x>=-D&x<=D)-h*(x<-D)4．636263 0S=2124822 ii==++++++∑①试不采用循环的形式，用数值方法求出上式的解。

②由于数值方法采用double形式进行计算，难以保证有效数字，试采用符号运算的方法求该式的精确解。

并给出保留16位有效数字的结果。

>> sum(2.^[1:63])>> sum(sym(2).^[1:63])5．编写一个矩阵相加函数mat_add（），使其具体的调用格式为，要求该函数能接受任意多个矩阵进行加法运算。