哈工大_Matlab__2014年秋季学期《MATLAB语言及应用》试题
- 格式:doc
- 大小:184.00 KB
- 文档页数:11
第5章习题与答案5.1用矩阵三角分解方法解方程组123123123214453186920x x x x x x x x x +-=⎧⎪-+=⎨⎪+-=⎩ 解答:>>A=[2 1 -1;4 -1 3;6 9 -1] A =2 1 -1 4 -13 6 9 -1 >>b=[14 18 20]; b =14 18 20 >> [L, U, P]=lu(A) L =1.0000 0 0 0.6667 1.0000 0 0.3333 0.2857 1.0000 U =6.0000 9.0000 -1.0000 0 -7.0000 3.6667 0 0 -1.7143 P =0 0 1 0 1 0 1 0 0 >> y=backsub(L,P*b’) y =20.0000 4.6667 6.0000 >> x=backsub(U,y) x =6.5000 -2.5000 -3.5000 5.2 Cholesky 分解方法解方程组123121332352233127x x x x x x x ++=⎧⎪+=⎨⎪+=⎩ 解答:>> A=[3 2 3;2 2 0;3 0 12] A =3 2 32 2 03 0 12>> b=[5;3;7]b =537>> L=chol(A)L =1.7321 1.1547 1.73210 0.8165 -2.44950 0 1.7321>> y=backsub(L,b)y =-11.6871 15.7986 4.0415>> x=backsub(L',y)x =-6.7475 28.8917 49.93995.3解答:观察数据点图形>> x=0:0.5:2.5x =0 0.5000 1.0000 1.5000 2.0000 2.5000 >> y=[2.0 1.1 0.9 0.6 0.4 0.3]y =2.0000 1.1000 0.9000 0.6000 0.4000 0.3000 >> plot(x,y)图5.1 离散点分布示意图从图5.1观察数据点分布,用二次曲线拟合。
matlab考试题目及答案1. 题目:编写一个MATLAB函数,实现计算并返回一个向量中所有元素的平方和。
答案:以下是一个简单的MATLAB函数,用于计算并返回一个向量中所有元素的平方和。
```matlabfunction sumOfSquares = calculateSumOfSquares(vector)sumOfSquares = sum(vector.^2);end```2. 题目:给定一个3x3的矩阵A,使用MATLAB编写代码,求出矩阵A 的转置。
答案:可以通过简单的转置操作来求得矩阵A的转置。
```matlabA = [1 2 3; 4 5 6; 7 8 9];A_transpose = A';```3. 题目:编写一个MATLAB脚本,实现对一个二维数组进行排序,并返回排序后的数组。
答案:以下是一个MATLAB脚本,用于对一个二维数组进行排序,并返回排序后的数组。
```matlabfunction sortedArray = sort2DArray(array)sortedArray = sort(array(:));end```4. 题目:给定一个向量x,使用MATLAB编写代码,计算并返回向量x的元素个数。
答案:可以通过内置函数`numel`来计算向量x的元素个数。
```matlabx = [1, 2, 3, 4, 5];numElements = numel(x);```5. 题目:编写一个MATLAB函数,实现计算并返回两个向量元素的点积。
答案:以下是一个简单的MATLAB函数,用于计算两个向量的点积。
```matlabfunction dotProduct = calculateDotProduct(vector1, vector2)dotProduct = dot(vector1, vector2);end```6. 题目:给定一个矩阵B,使用MATLAB编写代码,求出矩阵B的行列式。
《MATLAB语言及应用》期末大作业题目1.数组的创建和访问(20分,每小题2分):1)利用randn函数生成均值为1,方差为4的5*5矩阵A;2)将矩阵A按列拉长得到矩阵B;3)提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的矩阵C;4)寻找矩阵A中大于0的元素;]5)求矩阵A的转置矩阵D;6)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E;7)删除矩阵A的第2列和第4列得到矩阵F;8)求矩阵A的特政值和特征向量;9)求矩阵A的每一列的和值;10)求矩阵A的每一列的平均值;程序代码:clear;clc;A=1+sqrt(4)*randn(5) %生成均值为1,方差为4的5*5矩阵A;B=A(:) %将矩阵A按列拉长得到矩阵B;C=A([2 3],[2 4]) %提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的矩阵C;n=find(A>0) %寻找矩阵A中大于0的元素;x=A(n)D=A' %求矩阵A的转置矩阵D;E1=flipud(A); %对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E;E=fliplr(E1)F=A(:,[1 3 5]) %删除矩阵A的第2列和第4列得到矩阵F;[Av,Ad]=eig(A) %求矩阵A的特征值和特征向量;S=sum(A,1) %求矩阵A的每一列的和值;Avg=S/5 %求矩阵A的每一列的平均值;运行结果:A =2.3333 2.1171 0.8568 2.1971 -0.7526-1.7853 0.4453 -3.8292 1.2944 0.4690 -1.6011 -1.5874 -0.3887 0.7971 0.3448 -0.2100 -0.7769 -1.7828 -4.2700 -1.3165 -1.9771 -0.9730 1.6593 1.0561 2.1601B =2.3333-1.7853-1.6011-0.2100-1.97712.11710.4453-1.5874-0.7769-0.97300.8568-3.8292-0.3887-1.78281.65932.19711.29440.7971-4.27001.0561-0.75260.46900.3448-1.31652.1601C =0.4453 1.2944-1.5874 0.7971n =167111516171820222325x =2.33332.11710.44530.85681.65932.19711.29440.79711.05610.46900.34482.1601D =2.3333 -1.7853 -1.6011 -0.2100 -1.97712.1171 0.4453 -1.5874 -0.7769 -0.97300.8568 -3.8292 -0.3887 -1.7828 1.65932.1971 1.2944 0.7971 -4.2700 1.0561-0.7526 0.4690 0.3448 -1.3165 2.1601E =2.1601 1.0561 1.6593 -0.9730 -1.9771-1.3165 -4.2700 -1.7828 -0.7769 -0.21000.3448 0.7971 -0.3887 -1.5874 -1.60110.4690 1.2944 -3.8292 0.4453 -1.7853-0.7526 2.1971 0.8568 2.1171 2.3333F =2.3333 0.8568 -0.7526-1.7853 -3.8292 0.4690-1.6011 -0.3887 0.3448-0.2100 -1.7828 -1.3165-1.9771 1.6593 2.1601Av =Columns 1 through 40.1004 + 0.2832i 0.1004 - 0.2832i0.6302 -0.5216-0.5969 -0.5969 -0.4811 0.0856 -0.4405 + 0.0006i -0.4405 - 0.0006i -0.3078 0.21200.2732 - 0.4899i 0.2732 + 0.4899i0.0244 -0.1780-0.0617 + 0.2024i -0.0617 - 0.2024i 0.5254 0.8025 Column 50.3903-0.49590.02180.1929-0.7511Ad =Columns 1 through 4-2.6239 + 1.7544i 0 0 00 -2.6239 - 1.7544i 0 00 0 -0.2434 00 0 0 3.54550 0 0 0Column 50 0 0 02.2257S =-3.2403 -0.7749 -3.4846 1.0747 0.9049Avg =-0.6481 -0.1550 -0.6969 0.2149 0.18102.符号计算(10分,每小题5分):1) 求方程组20,0uy vz w y z w ++=++=关于,y z 的解;2) 利用dsolve 求解偏微分方程,dx dyy x dt dt==-的解;程序代码:clc[u,v,w] = solve('u*y^2 + v*z + w = 0','y + z + w = 0','u,v,w')[x y]=dsolve('Dx=y','Dy=-x')运行结果:u =(-v*z+y+z)/y^2v =vw = -y-zx =-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:4a =;程序代码:clccleara=0.5:0.5:4.5;t=-2*pi:0.1:2*pi;N=length(a); for i=1:1:Nx=a(i)*cos(t);y=sqrt(25-a(i).^2).*sin(t);plot(x,y)hold onend运行结果:-5-4-3-2-1123452) 利用plotyy 指令在同一张图上绘制sin y x =和10x y =在[0,4]x ∈上的曲线;程序代码:clcclearx=0:0.01:4;y1=sin(x); y2=10.^x;plotyy(x,y1,x,y2) %用双y 轴绘制二维图形运行结果:3) 用曲面图表示函数22z x y =+;程序代码:clcclear[X,Y] = meshgrid(-2:0.05:2); %产生xy 平面上的网格数据Z = X.^2 + Y.^2;surf(X,Y,Z) %绘着色曲面图hold off运行结果:4) 用stem 函数绘制对函数cos 4y t π=的采样序列;程序代码:clcclearfs=25;Ts=1/fs;n=1:1:200;yn=cos(pi*n*Ts/4);stem(n,yn) %绘离散数据的火柴杆图运行结果:4. 设采样频率为Fs = 1000 Hz ,已知原始信号为)150π2sin(2)80π2sin(t t x ⨯+⨯=,由于某一原因,原始信号被白噪声污染,实际获得的信号为))((ˆt size randn x x +=,要求设计出一个FIR 滤波器恢复出原始信号。
计算机实践1 Matlab基础考试题目1、考试为上机操作。
2、每个考试同学考试题号由程序随机现场产生,为使每个考生之间题目难易程度相当,考生每次抽取的题目都是从同类型题目中抽取。
每位同学应做的题号在本文件的最后两页3、操作过程,请注意保存结果。
请先建立以自己“姓名拼音+学号”命名的M文件,例如:lishaohong1340502113。
在该M文件中编辑答案,每题以%第**题开始。
例如,第57题,以“%第57题”开始。
4、本测试并不确保每一个题目数学上正确,每个同学自己去判断,自己认为是怎样的结果,就写怎样的结果。
5、操作完成,将文件以附件发送至10307086@并告知老师已发送可以看书,可以看笔记,可以上网查资料。
禁止交流、相互传阅资料。
如发现,视为作弊。
一、简单计算值的计算机实现:1. i i e i i 34)3sec(7ln 26))92()2sin()8945)(233((43++-+-+π2. )94)14sec(27ln 259)(8945)(233()92()2sin(63i i i i e ++-+--π 3. i e i i 379)59sec(34ln 1116))2sin()879165)(533((43+++-+π 4.ie i i 3)59csc(34log 1611))2cos())()3/5()5^3((sin(2113+⨯++π5.i e i 3)2cot(34log 1619)2)4cos())3/1()6^2((tan(612+⨯⋅++π二、求极限、微分和积分6. 求7. 求函数)5(/)3sin(x tg x y =在0=x 处的极限8. 求函数)()sin(π-=x x y 在π=x 处的极限9. 求函数 10. 求xx )tan(,当0→x 时的极限11. 求2/1))sin((x e x y x +=的3阶导数 12. 求)3cos(2/12x e y x -=的4阶导数 13. 求33/12+-=x x y 的50阶导数 14. 求的2阶导数 15. 求的3阶导数 16. 求b x x ++)sin(在]8,0[=x 的定积分17. 求)1/(1)(2x x g +=在],1[∞=x 的定积分 18. 求2/12)(a x +在]2,2[-=x 的定积分?)()cos(1lim )(22220022=++-+→→y x y xey x y x ?)(lim 3221=++→-→y x xy y x y x xe x x x y -=1sin )(x x x x y ln 1cos )(-=19. 求 20. 求?ln 10=⎰dt t x21. 由表达式()2ln f x y =+,求/f x ∂∂,/f y ∂∂,2/f x y ∂∂∂22. 由表达式)ln(y x xy f +=,求/f x ∂∂,/f y ∂∂,2/f x y ∂∂∂23. 由表达式3y x ey +=,求/f x ∂∂,/f y ∂∂,2/f x y ∂∂∂24. 由表达式y)cos()sin(x f =,求/f x ∂∂,/f y ∂∂,2/f x y ∂∂∂25. 由表达式)cos(/)sin(x y f =,求/f x ∂∂,/f y ∂∂,2/f x y ∂∂∂26. 计算x x f +=)sin(表达式在/6x π=处的值,并将结果设置为以下精度:小数点之后 1 位、5 位和20位有效数字27. 计算1)cos(+=x f 表达式在/6x π=处的值,并将结果设置为以下精度:小数点之后 3 位、6 位和14位有效数字28. 计算)sin()tan(x x f =表达式在/6x π=处的值,并将结果设置为以下精度:小数点之后 2位、7 位和12位有效数字29. 计算)cos(/)sin(x x f =表达式在/6x π=处的值,并将结果设置为以下精度:小数点之后 1 位、8位和22位有效数字30. 计算)sin()cos(x x f +=表达式在/6x π=处的值,并将结果设置为以下精度:小数点之后 1 位、4 位和10位有效数字 三、向量、矩阵运算31. 已知A=[1 2 3];B=[3 1 4];C=[9 -1 4]求A 和B 的点积,B 和C 的叉积 32. 计算向量a=(4, 2, 8)、b=(4, 9, 6)和c=(-2, 6, -1) 的混合积)(c b a ⨯∙?cos 0==⎰∞dx xx I33. 计算向量a=(6, 1, 8)、b=(3, 9, 5)和c=(-1, 3, -5) 的混合积c b a ⨯∙)( 34. 计算向量a=(5, 2, 1)、b=(2, 4, 3)和c=(1, 2,3) 的混合积c b a ∙⨯)( 35. 计算向量a=(3, 2,6)、b=(3, 5, 7)和c=(5, 3,1) 的混合积c b a ⨯∙)(36. 求矩阵⎪⎪⎪⎭⎫ ⎝⎛--=314020112A 的特征值和特征向量 37. 求解4阶随机矩阵的特征值和特征向量38. 设A 为6阶魔方矩阵,分别对A 求逆,求行列式的值,求迹 39. 设B 为4阶魔方矩阵,分别对B 求逆,求行列式的值,求秩40. 求矩阵⎪⎪⎪⎭⎫ ⎝⎛--=623212113A 的特征值、特征向量和迹 四、求解方程组41. ⎪⎪⎩⎪⎪⎨⎧-=++=++=++=+-550754322232132132121x x x x x x x x x x x42. ⎪⎪⎩⎪⎪⎨⎧=-+-=++-=+--=-++0422123432143214324321x x x x x x x x x x x x x x x43. ⎪⎩⎪⎨⎧-=+---=+-=-+-42344324214324321x x x x x x x x x x44. ⎪⎩⎪⎨⎧=++=++=++113212523321321321x x x x x x x x x45.⎪⎩⎪⎨⎧=++=+=++1041023433232121321xxxxxxxx五、数据统计处理46.已知A=34*rand(12,10) 分别求矩阵A中各列和各行元素中的最大值,并求整个矩阵的最大值和最小值。
M a t l a b期末复习题(答案)-CAL-FENGHAI.-(YICAI)-Company One1一、 选择题(有10个题,共20分)1.有一个4×4矩阵A,表示取得矩阵A 的第一行和第二行的第二列至第四列所有元素的命令是( A )(A) A(1:2,2:4) (B) A(2:4,1:2) (C) A(1:2:4) (D) A(4:2:1)2.将一个图形窗口分成3×4个子区并表示第一个子区的命令是( C )A . subplot(12,1)B .subplot(1,12)C .subplot(3,4,1)D .subplot(1,3,4)3.利用Matlab 的符号计算功能计算微分方程初值问题⎪⎩⎪⎨⎧=+-=2)1(42y x x y dx dy 在区间 [1,2]上的解析解时所用到的函数是( D )(A) ode23 (B) ode45 (C) solve (D) dsolve4.以下有关Matlab 子函数的性质中不正确的是( D )(A) 子函数只能被处于同一个文件的主函数或其他子函数调用;(B) 在M 函数文件中,任何指令通过“名字”对函数进行调用时,子函数的优先级仅次于内装函数;(C) 在M 函数文件内,主函数的位置不可改变,但子函数的排列次序可以任意改变;(D) help, lookfor 等帮助指令都能提供关于子函数的任何帮助信息。
5.以下语句中生成按钮控件的语句是( B )(A) h_1=uimenu(gcf,'Label','&Blue');(B) h_1=uicontrol(gcf,'style','push', 'string','grid off','callback','grid off');(C) h_1=uicontrol(gcf,'style','text', 'horizontal','left','string',{'输入'});(D) h_1=axes('unit','normalized','position',[0,0,1,1],'visible','off');6.以下有关Matlab 函数及其句柄的说法中正确的是( A,C )(A) 在该函数的目录下可定义该函数的句柄,句柄定义后,可用feval 函数和句柄在任何目录下调用该函数;(B) 函数不在Matlab 的当前工作目录下不可执行,但可以用feval 函数和句柄调用该函数;(C) 函数不在Matlab 的当前工作目录下时,可以定义函数的句柄;(D) 函数不在Matlab 的当前工作目录下时,不可以定义函数的句柄。
第一套华东交通大学2015-2016学年第一学期Matlab期末考试一、填空题(4*6=24分)1、在MATLAB^令窗口中的“ >>”标志为MATLAB勺_提示符,标志为提示符。
2、符号表达式sin(2*a+t)+m 中独立的符号变量为__ 。
3、在通常情况下,左除x=a\b是—的解,右除x=b/a是—的解,一般情况下,_______________ 。
4、为了使两个plot的图形在同一个坐标显示,可以使用 _______ 令进行图形保持;可以使用—命令为图形添加网格。
5、倘若要是对x进行赋值,从5到25,间隔是0.1,为_________ ;倘若要是对x进行赋值,从2到20,中间间隔100个点,为 ______6、A=[1,2;3,4];B=[1,0;0,1]; A*B= _ , A.*B=_ _ ;二、选择题(4*6=24分)7、如果x=1: 2 : 8,则x(1)和x⑷分别是()A. 1,8B. 1, 7C. 2, 8D. 2, 78、运行如下程序后,输入9回车,命令窗口( comma nd win dows)显示的结果为()c=input('请输入一个字符','s');if c>='A' & c<='Z'disp(setstr(abs(c)+abs('a')-abs('A')));elseif c>='a'& c<='z'disp(setstr(abs(c)- abs('a')+abs('A')));elseif c>='0' & c<='9'disp(abs(c)-abs('0'));else disp(c); endA. 9B. 8C. 7D. 69、MATLA表达式2*2A3A2的结果是()A. 128B. 4096C. 262144D. 25610、在循环结构中跳出循环,执行循环后面代码的命令为()(A) return (B) break (C) con ti nue (D) keyboard11、在图形指定位置加标注命令是()A. title(x,y, ^y=sin(x) ?B. xlabel(x,y, ?y=s in(x)C. text(x,y, ?y=sin(x) ?);D. lege nd(x,y, ?y=s in(x)12、下列哪个变量的定义是不合法的()。
《Matlab语言及其应用》形考作业1-3答案形考作业1一、判断题(8小题,每小题2分,共16分)1.Matlab提供了两种创建符号变量和表达式的函数:sym和syms。
(A)判断题(2 分) 2分A.对B.错2.xy12_是合法的Matlab变量名。
(A)判断题(2 分) 2分A.对B.错3.Matlab中的变量必须遵循先定义后使用的原则。
(B)判断题(2 分) 2分A.对B.错4.Matlab的含义为矩阵实验室(Matrix Laboratory)。
(A)判断题(2 分) 2分A.对B.错5.Matlab中注释行的开头是#。
(B)判断题(2 分) 2分A.对B.错6.plot( )是Matlab最基本的绘图函数。
(A)判断题(2 分) 2分A.对B.错7.双纵坐标图形的绘制函数是plotyy()。
(A)判断题(2 分) 2分B.错8.在Matlab中,二维饼图由函数pie( )来实现。
(A)判断题(2 分) 2分A.对B.错二、单选题(9小题,每小题3分,共27分)9.清空Matlab工作空间内所有变量的指令是(C)。
单选题(3 分) 3分A.clcB.clsC.clearD.clf10.在循环结构中跳出循环,执行循环后面代码的命令为(B)。
单选题(3 分) 3分A.returnB.breakC.continueD.keyboard11.图形窗口的分割函数subplot(2,2,[1 3])可以把一个窗口分割为(B)个子窗口。
单选题(3 分) 3分A.2B.3C.4D.无法确定12.一般情况下,绘图命令每执行一次就刷新当前图形窗口,图形窗口原有图形将不复存在。
若希望在已绘图形上叠加新的图形,可以使用(A)命令。
单选题(3 分) 3分A.hold onB.hold offD.plot13.Matlab中求极限的函数为(B)。
单选题(3 分) 3分A.diff( )B.limit( )C.int( )D.solve( )14.Matlab中求导数的函数为(A)。
2013年春季学期《MATLAB语言及应用》课程试卷姓名:学号:学院:专业:必答题1.常用的matlab界面由哪些窗口组成,各有什么主要作用?(4分)(1)菜单和工具栏功能:包括【File】【Edit】【Debug】等对文件,程序等进行操作(2)命令窗口功能:用于输入命令并显示除图形以外的所有执行结果(3)历史命令窗口功能:主要用于记录所有执行过的命令(4)当前工作目录窗口功能:对文件和目录进行操作(5)工作空间窗口功能:查看、载入和保存变量2.如何设置当前目录和搜索路径,在当前目录上的文件和在搜索路径上的文件有什么区别?(2分)方法一:在MATLAB命令窗口中输入editpath或pathtool命令或通过【File】/|【SetPath】菜单,进入“设置搜索路径”对话框,通过该对话框编辑搜索路径。
方法二:在命令窗口执行“path(path,…D:\Study ‟)”,然后通过“设置搜索路径”对话查看“D:\Study”是否在搜索路径中。
方法三:在命令窗口执行“addpath D:\Study- end”,将新的目录加到整个搜索路径的末尾。
如果将end改为begin,可以将新的目录加到整个搜索路径的开始。
区别:当前文件目录是正在运行的文件的目录,显示文件及文件夹的详细信息,且只有将文件设置为当前目录才能直接调用。
搜索路径中的文件可以来自多个不同目录,在调用时不用将其都设置为当前目录,为同时调用多个文件提供方便。
3.有几种建立矩阵的方法?各有什么优点?(4分)1.在命令窗口中直接输入优点:适合输入不规则和较小的矩阵2.通过语句和函数生成矩阵优点:适合输入规则且较大的矩阵3.通过M文件来建立矩阵,从外部数据文件中导入矩阵优点:方便创建和导入大型矩阵4.说明break语句、continue语句和return语句的用法。
(3分)1. break语句用于终止for循环和while循环的执行。
当遇到break时,则退出循环体继续执行循环体外的下一个语句。
在嵌套循环中,break往往存在于内层的循环中。
2.continue命令用于结束本次for或while循环,与break命令不同的是continue指结束本次循环而继续进行下次循环。
3. return语句用来终止被调用函数的运行,后面的程序代码将不再执行,直接返回到上一级调用函数。
5.GUI开发环境中常用的控件包括哪些,各有何作用?(4分)①按钮(Push Buttons):执行某种预定的功能或操作;②开关按钮(Toggle Button):产生一个动作并指示一个二进制状态(开或关),当鼠点击它时按钮将下陷,并执行callback(回调函数)中指定的内容,再次点击,按钮复原,并再次执行callback 中的内容;③单选框(Radio Button):单个的单选框用来在两种状态之间切换,多个单选框组成一个单选框组时,用户只能在一组状态中选择单一的状态,或称为单选项;④复选框(Check Boxes):单个的复选框用来在两种状态之间切换,多个复选框组成一个复选框组时,可使用户在一组状态中作组合式的选择,或称为多选项;⑤文本编辑器(Editable Texts):用来使用键盘输入字符串的值,可以对编辑框中的内容进行编辑、删除和替换等操作;⑥静态文本框(Static Texts):仅用于显示单行的说明文字;⑦滚动条(Slider):可输入指定范围的数量值;⑧边框(Frames):在图形窗口圈出一块区域;⑨列表框(List Boxes):在其中定义一系列可供选择的字符串;⑩弹出式菜单(Popup Menus):让用户从一列菜单项中选择一项作为参数输入;⑩坐标轴(Axes):用于显示图形和图象。
6.进行simulink建模仿真的主要步骤包括哪些?(4分)(1)根据要仿真的系统框图,在SIMULINK 窗口的仿真平台上构建仿真模型。
(2)设置模块参数。
(3)设置仿真参数。
(4)启动仿真。
(5) 观测仿真结果。
7.a=[0 0.2 1 0 -5 9]在进行逻辑运算时,a相当于什么样的逻辑量。
(2分)解:a=[0,1,1,0,1,1],a是一个逻辑向量。
8.生成4阶魔方矩阵A,再分别计算A的数组平方和矩阵平方,给出执行语句以及计算结果。
(3分)解:>> A = magic(4)B = A^2SSR = sum>> A = magic(4)B = A^2SSR = sum(B(:))MSR = A^2A =16 2 3 135 11 10 89 7 6 124 14 15 1B =256 4 9 16925 121 100 6481 49 36 14416 196 225 1SSR =1496MSR =345 257 281 273257 313 305 281281 305 313 257273 281 257 345 (B(:))MSR = A^2A =16 2 3 135 11 10 89 7 6 124 14 15 1B =256 4 9 16925 121 100 6481 49 36 14416 196 225 1SSR =1496MSR =345 257 281 273257 313 305 281281 305 313 257273 281 257 3459.生成一个8阶的魔方矩阵A,再采用相应的MATLAB命令,将其全部偶数行提取出来,赋给矩阵B,给出命令及结果(4分)M = magic(8);>> B = M(~mod(M,2))A =64 2 3 61 60 6 7 579 55 54 12 13 51 50 1617 47 46 20 21 43 42 24 40 26 27 37 36 30 31 33 32 34 35 29 28 38 39 25 41 23 22 44 45 19 18 48 49 15 14 52 53 11 10 56 8 58 59 5 4 62 63 1 B =9 55 54 12 13 51 50 16 40 26 27 37 36 30 31 33 41 23 22 44 45 19 18 48 8 58 59 5 4 62 63 110. 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=986810141248a ,计算a 的行列式和逆矩阵。
(4分) >> a=[8 4 -12;14 11 5 ;6 8 9];>> ad=det(a) ad =-728>> ai=inv(a) ai =-0.0357 0.1813 -0.2088 0.1071 -0.1978 0.3187 -0.0714 0.0549 -0.033011. 对于B AX =,如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753487294A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=183657B ,求解X 。
(3分) >> A=[4 9 2;7 18 4;3 5 7];>> B=[57 ;36 ;18]; >> X=A\B X = -3.0180 8.1018 -1.922212. 试求出如下极限:(4分)(1)xxx x 1)93(lim +∞→; (2)11lim0-+→→xy xy y x(1)解:>> clear >> syms x;>> limit(((x^3+x^6)^(1/x)),x,inf) ans = 9(2)解:>> syms x y; >> t=4*x*yt = 4*x*y>> limit(limit(t/(sqrt(t+1)-1),x,0),y,0)ans = 213. 求多项式43810x x +-乘以、除以多项式223x x -+的结果。
(4分) 求乘积>> p=[1,8,0,0,-10]; >> q=[0,0,2,-1,3]; >> k=conv(p,q) k =0 0 2 15 -5 24 -20 10 -30 结果为2x^6+15x^5-5x^4+24x^3-20x^2+10x-30 相除 >> q=[2,-1,3]; >> [k,r]=deconv(p,q) k =0.5000 4.2500 1.3750 r =0 0 0 -11.3750 -14.125043810x x +-=(223x x -+)*(0.500x^2+4.250x+1.375)+(-11.375x-14.125)14. 根据x= [66 49 71 56 38 46 90]绘制饼图,并将第五个切块分离出来。
(3分) >> x=[23 71 12 59 41 36 16 55]; >> L=[0 0 0 0 1]; >> pie(x,L);编写一个函数,使其能够产生如下的分段函数:⎪⎩⎪⎨⎧≤≤<-≤=x x x x x x f 65.0620.251.525.0)(,,,,并调用此函数,绘制在x=[0,2]范围内的()(2)f x f x ∙+曲线。
(4分) M 函数 function y=f(x)y=(x<=2).*(0.5*x)+(x>=2&x<=6).*(1.5-0.25*x)+(x>3).*(0.5); 函数图像如下 >> x=0:0.1:2; >> y=f(x).*f(x+2);>> plot(x,y),grid on15. 假定某天的气温变化记录如下表,试作图描述这一天的气温变化规律,再用二次多图像如下二次多项式拟合 >> pf=polyfit(t,w,2) pf =-0.1081 3.0257 5.7159 气温w=-0.1081t^2+3.0257t+5.715916. 计算z=yx 2+3y 2x+2y 3的z y ∂∂和2z y x∂∂∂。
(3分)>> clear;>> syms x y;>> z=y*x^2+3*y^2*x+2*y^3; >> dt=diff(z,y); >> dxy=diff(dt,x); >> dt,dxy dt =x^2 + 6*x*y + 6*y^2 dxy = 2*x + 6*yz y ∂∂= x^2 + 6*x*y + 6*y^2 2z y x∂∂∂=2*x + 6*y17. 用符号计算验证三角等式:sin(ϕ1)cos(ϕ2)-cos(ϕ1)sin(ϕ2) =sin(ϕ1-ϕ2)。