MATLAB练习作业(已做)2012
matlab2012实验1参考答案
MATLAB 实验一 MATLAB 数值计算试验报告说明：1 做试验前请先预习，并独立完成试验和试验报告。
2 报告解答方式：将MATLAB 执行命令和最后运行结果从命令窗口拷贝到每题的题目下面，请将报告解答部分的底纹设置为灰色，以便于批阅。
3 在页眉上写清报告名称，学生姓名，学号，专业以及班级。
3 报告以Word 文档书写。
文档命名方式: 学号+姓名+_（下划线）+试验几.doc 如：110400220张三_试验1.doc 4 试验报告doc 文档以附件形式发送到maya_email@ 。
凡文档命名不符合规范，或者发送方式不正确，不予登记。
5 每次试验报告的最后提交期限：下次试验课之前。
一 目的和要求1 熟练掌握MATLAB 变量的使用2 熟练掌握矩阵的创建3 熟练掌握MATLAB 的矩阵和数组的运算4 使用元胞数组和结构数组二 试验内容1 创建矩阵（必做）1.1使用直接输入，from:step:to ，linspace ，logspace 等方式创建矩阵。
1.2 输入矩阵1234246836912a ⎛⎫⎪= ⎪ ⎪⎝⎭1.2-1）分别使用全下标和单下标达方式取出元素“8”>>a=[1 2 3 4;2 4 6 8;3 6 9 12]>> a(2,4) %全下标方式 >> a(11) % 单下标方式1.2-2）分别用不同的方式从矩阵a 中取出子矩阵⎥⎦⎤⎢⎣⎡1286432%方法一：全下标方式 a([2,3],[1 2 4])%方法二：单下标方式 a([2 5 11;3 6 12])% 方法三：利用逻辑向量 l1=logical([0 1 1])l2=logical([1 1 0 1]) a(l1,l2)1.2-3）分别写出以下的执行结果：a ( :, 3 ) a (:,end) a (:, 2 : 3) a ([1 3], : )a (:) a (2 : 3) a([2 3 4])本题解答略，希望大家熟练掌握数组运算中 ： end 的用法以及 全下标，单下标等方式2 矩阵的运算（必做）2.1使用format 命令，将 pi 分别用short 和long 格式显示其值。
(完整版)matlab基础练习题(带答案)
Matlab 基础练习题常量、变量、表达式1、 MATLAB 中，下面哪些变量名是合法的？（ ）（A ）_num （B ）num_ （C ）num- （D ）-num 2、 在MA TLAB 中，要给出一个复数z 的模，应该使用（ ）函数。
（A ）mod(z) （B ）abs(z) （C ）double(z) （D ）angle(z) 3、 下面属于MATLAB 的预定义特殊变量的是？（ ）（A ）eps （B ）none （C ）zero （D ）exp4、 判断：在MA TLAB 的内存工作区中，存放一个英文字符 'a' 需要占用1个字节，存放一个中文字符‘啊’需要占用2个字节。
（ 错，都是2个字节 ）5、 判断：MA TLAB 中，i 和j （ 对 ）6、 判断：MA TLAB 中，pi 代表圆周率，它等于3.14。
（ 错，后面还有很多位小数 ）7、 在MA TLAB 中，若想计算的51)3.0sin(21+=πy 值，那么应该在MA TLAB 的指令窗中输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。
8、 在MA TLAB 中，a = 1，b = i ，则a 占_8__个字节，b 占_16_个字节，c 占________字节。
9、 在MA TLAB 中，inf 的含义是__无穷大__，nan 的含义是__非数（结果不定）___。
数组1、 在MA TLAB 中，X 是一个一维数值数组，现在要把数组X 中的所有元素按原来次序的逆序排列输出，应该使用下面的（ ）指令。
（A ）X[end:1] （B ）X[end:-1:1] （C ）X (end:-1:1) （D ）X(end:1) 2、 在MA TLAB 中，A 是一个字二维数组，要获取A 的行数和列数，应该使用的MATLAB的命令是（ ）。
（A ）class(A) （B ）sizeof(A) （C ）size(A) （D ）isa(A)3、 在MATLAB 中，用指令x=1:9生成数组x 。
#数值分析Matlab作业
数值分析编程作业2012年12月第二章14.考虑梯形电阻电路的设计，电路如下：电路中的各个电流{i1,i2,…,i8}须满足下列线性方程组：121232343454565676787822/252025202520252025202520250i i V R i i i i i i i i i i i i i i i i i i i i -=-+-=-+-=-+-=-+-=-+-=-+-=-+=这是一个三对角方程组。
设V=220V ，R=27Ω，运用追赶法，求各段电路的电流量。
Matlab 程序如下：function chase () %追赶法求梯形电路中各段的电流量 a=input('请输入下主对角线向量a='); b=input('请输入主对角线向量b='); c=input('请输入上主对角线向量c='); d=input('请输入右端向量d='); n=input('请输入系数矩阵维数n='); u(1)=b(1); for i=2:nl(i)=a(i)/u(i-1); u(i)=b(i)-c(i-1)*l(i); endy(1)=d(1); for i=2:ny(i)=d(i)-l(i)*y(i-1); endx(n)=y(n)/u(n); i=n-1; while i>0x(i)=(y(i)-c(i)*x(i+1))/u(i); i=i-1;end x输入如下: >> chase请输入下主对角线向量a=[0,-2,-2,-2,-2,-2,-2,-2]; 请输入主对角线向量b=[2,5,5,5,5,5,5,5];请输入上主对角线向量c=[-2,-2,-2,-2,-2,-2,-2,0]; 请输入方程组右端向量d=[220/27,0,0,0,0,0,0,0]; 请输入系数矩阵阶数n=8 运行结果如下：x = 8.1478 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194 0.0477第三章14.试分别用(1)Jacobi 迭代法；(2)Gauss-Seidel 迭代法解线性方程组1234510123412191232721735143231211743511512x x x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥---⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥---⎣⎦⎣⎦⎣⎦ 迭代初始向量(0)(0,0,0,0,0)T x =。
2012Matlab试题及解答
1. 要在命令窗口中查看对数函数log 的帮助信息，所用命令是help log 。
2. 删除工作空间中保存的变量A 的命令是clear A 。
3. 在Matlab 命令窗口输入（2*sin(85*pi/180))/(1+exp(2))可以计算表达式22sin 851e +的值。
4. 在同一图形窗口中开多个子窗口的Matlab 函数名是subplot 。
5. 设a 是一维数组，a 的倒数第3个元素表示为 a （end-2）。
6 求多项式43810xx --的根的命令为roots([1,8,0,0,-10])。
二、阅读程序并回答问题（本大题共3小题，每小题5分，共15分） 7. A=zeros(2,4); A(:)=1:8;s=[2 3 5];Sa=[10 20 30]'; A(s)=Sa ;AA =1 20 30 7 10 4 6 8..............................5分 8. 写出下列程序的输出结果。
B{1,1}=1; B{1,2}='Breden'; B{2,1}=reshape(1:9,3,3); B{2,2}={12,34,2;54,21,3;4,23,67}; B(3)=[]; B B =[1] [3x3 double] {3x3 cell} 9. 写出下列程序的输出结果： s=0；a=[12,13,14;15,16,17;18,19,20;21,22,23]; for k=a for j=1:4if rem(k(j),2)~=0s=s+k(j);endendendss =108..............................5分三、编写程序解答下列各题（本大题共4小题，每小题8分，共32分）10.从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们的和。
sum=0;n=0;x=input('Enter a number (end in 0):');while (x~=0)sum=sum+x;n=n+1;x=input('Enter a number (end in 0):');endif (n>0)summean=sum/nend11. 一个3位数各位数字的立方和等于该数本身则称该数为水仙花数。
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开始，便是函数⽂件，每⼀个函数⽂件定义⼀个函数。
(完整版)matlabR2012a课后习题答案第三章
第3章数值数组及其运算习题3及解答1 要求在闭区间]2,0[ 上产生具有10个等距采样点的一维数组。
试用两种不同的指令实现。
〖目的〗●数值计算中产生自变量采样点的两个常用指令的异同。
〖解答〗%方法一t1=linspace(0,2*pi,10)%方法二t2=0:2*pi/9:2*pi %要注意采样间距的选择，如这里的2*pi/9.t1 =Columns 1 through 70 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 Columns 8 through 104.88695.58516.2832t2 =Columns 1 through 70 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 Columns 8 through 104.88695.58516.28322 由指令rng('default'),A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5的元素的位置，分别求出它们的“全下标”和“单下标”。
〖目的〗●数组下标的不同描述：全下标和单下标。
●sub2ind, int2str, disp的使用。
●随机发生器的状态控制：保证随机数的可复现性。
〖解答〗rng('default')A=rand(3,5)[ri,cj]=find(A>0.5);id=sub2ind(size(A),ri,cj);ri=ri';cj=cj';disp(' ')disp('大于0.5的元素的全下标')disp(['行号 ',int2str(ri)])disp(['列号 ',int2str(cj)])disp(' ')disp('大于0.5的元素的单下标')disp(id')A =0.8147 0.9134 0.2785 0.9649 0.9572 0.9058 0.6324 0.5469 0.1576 0.4854 0.1270 0.0975 0.9575 0.9706 0.8003大于0.5的元素的全下标行号 1 2 1 2 2 3 1 3 1 3 列号 1 1 2 2 3 3 4 4 5 5大于0.5的元素的单下标1 2 4 5 8 9 10 12 13 153 采用默认全局随机流，写出产生长度为1000的“等概率双位（即取-1，+1）取值的随机码”程序指令，并给出 -1码的数目。
Matlab基本操作练习
实验一Matlab基本操作练习一、实验目的（1）了解和学会使用matlab平台的命令窗口、帮助窗口、编辑窗口、工作空间。
（2）学会使用matlab的.m文件、.mat文件（3）学会基本的矩阵赋值、运算和查找等操作（4）了解for循环、while循环的执行过程（5）掌握函数的定义和调用的基本方法二、实验内容1、了解save命令和load命令的使用。
操作步骤：（1）在命令窗口输入如下命令：clear ;savefile = 'TestFile.mat';p = rand(1, 4);q = ones(10);r = ones(10,2);save(savefile, 'p', 'q', 'r')（2）观察并记录下当前工作空间中的变量；（3）输入如下命令：clear;（4）观察当前工作空间中的变量；（5）在命令窗口输入如下命令：load TestFile（6）观察当前工作空间中的变量；2、练习基本的数学运算（1）逐句输入并运行以下命令，观察每一步的运行结果，弄明白每一条语句的含义。
（2）计算A*inv(C)，将结果与divrAC比较计算A*inv(C)，将结果与divrAC比较3、练习矩阵定位函数find的使用逐句输入并运行以下命令，观察每一步的运行结果，弄明白每一条语句的含义。
clear;X = [1 0 4 -3 0 0 0 8 6]Indices1 = find(X)indices2=find(X>=6)clear;A=[3 2 0; -5 0 1]Indices3=find(A)[row,col,value]=find(A)4、练习for循环和while循环的使用操作步骤：（1）在命令窗口输入如下命令clear ;clc;（2）在命令窗口输入如下命令，并观察s1的值：a=1:100;s1=0;for k=1:100s1=s1+k;end（3）在命令窗口输入如下命令，并观察s2的值：s2=0;k=1;while k<=100s2=s2+k;k=k+1;end（4）在命令窗口输入如下命令，并观察s3的值：s3=sum(a);[s1 s2 s3]5、练习函数的定义和调用；定义和调用一个连续求和的函数SumOfN，能实现功能：给一个参数n，求出从1到n连续求和的结果。
MATLAB Char03-数据可视化与绘图作业20120412
【作业1】绘制函数 y xe 和 y2=exsin(x) 在0 x 1 时的曲线。两条曲线分别用红色实 线和绿色虚线表示,并用legend函数进行标识。
x
【作业2】使用axis equal语句画出球面正弦波：
a 10.0, b 1.0, c 0.3, 0 t 2
绘出下述级数的图形。除非特别指明，所有级数仅求 2n n 1,3,5,...
1 1 2 2
4
【作业6】请绘制如下分段函数的图形，其中x的取值范围为 [-10,10]，并为图形的横坐标加上标签“x”，为图形的纵坐 标加上标签“f(x)”。
x 1 3 x 2, f x x, 1 x 1 2 x 3, x 1
【作业7】标准条件下钢管中流动的空气压降可由下式给出：
x cos(t ) b2 c 2 cos 2 (at ) y sin(t ) b2 c 2 cos 2 (at ) z c cos(at )
【作业3】用surf函数画出正螺旋面：
(c 1/ 2 , 0.5 u 0.5, 2 v 2 )
x u cos(v ) y u sin(v) z cv
【作业4】给定如下8点的x和y坐标向量，请由此7 点拟合成一条4次曲线方程，并在一图形窗口绘出 点和拟合曲线的图形。 x=[2 2.5 3 4 5 6 7 10] y=[4 3 5 6 8 5.6 9.8 12]
【作业5】从球面坐标到迪卡尔坐标的转换公式为: x=b sin a cos c y=b sin a sin c z=b cos a a 在[0 90]之间的范围内分10等份，c在[0 360]的范围内分24等 份。令b=2 ,计算x,y,z，并利用mesh(x,y,z)绘出这个半球形。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
说明:如在测试时请将程序中的中文字符改为西文(部分程序为手写,并未经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矩阵运算的学习。