习题集六(Matlab数值计算)课后习题集
MATLAB复习题
习
题
1、Simulink的仿真模型主要由 输入模块 、 状态模块 、输出模 块 三个部分构成。
2、建立仿真模型的过程主要由 模型结构分析、模块的复制、模块 连接、参数设置（模块参数及仿真参数）、运行仿真 五个步骤完 成。
习
1、提示从键盘输入变量的命令为：
题
input
2、在m文件中，用于暂停程序执行的命令为:
【 A 】
10、在MATLAB语言中，下列常用操作命令解释错误的是 A、dir为显示当前目录下的文件 C、path为显示搜索目录 E、clear为清除命令窗口中的内容 11、若只知要搜索函数的前两个字母，可用 B、clc为清除工作区中的变量 D、clf为清空当前图形窗口 F、path为显示搜索路径
【 BE
pause
函数文件
3、如果M文件的可执行语句以function 开始，则该文件为 定义全局变量 while if 和 和 for switch
5、在M 文件中，对流程控制实现循环结构的语句有 6、在M 文件中，对流程控制实现选择结构的语句有
【 B
】
7、在MATLAB语言中，提示其后语句为输入指令的行提示符为
C、“，” D、“>>”
【 D
】
8、在MATLAB语言中，在指令窗口运行过的指令再次调出运行的按键为 B、“↑” 或“↓” B、“：” C、“，” D、“>>”
【 B
】
9、在MATLAB语言中，若希望结果不被显示，可在语句之后加上符号 C、“，” D、“>>”
】
lookfor 命令搜索此关键字的帮助信息。
12、在MATLAB语言中，可采用
Tab
键进行模糊查询。
matlab基础测习题(带答案)
精心整理Matlab基础练习题常量、变量、表达式1、MATLAB中，下面哪些变量名是合法的？（）2、3、4、5、）6、小数7、8、在MATLAB中，a=1，b=i，则a占_8__个字节，b占_16_个字节，c占________字节。
9、在MATLAB中，inf的含义是__无穷大__，nan的含义是__非数（结果不定）___。
数组1、在MATLAB中，X是一个一维数值数组，现在要把数组X中的所有元素按原来次序的逆序排列输出，应该使用下面的（）指令。
（A）X[end:1] （B）X[end:-1:1] （C）X(end:-1:1)（D）X(end:1) 2、在MATLAB中，A是一个字二维数组，要获取A的行数和列数，应该使用的MATLAB的命令是（）。
（A）class(A) （B）sizeof(A) （C）size(A)（D）isa(A)3、4、这5、6、___linspace(0,5,50)___7、在MATLAB中，A=[0:1/2:2]*pi，那么sin(A)=___[010-10]_____。
8、在MATLAB中，A=[1,2,3;4,5,6;7,8,0]，B=[2,1,6;8,5,2;14,2,1]。
写出下面MATLAB语句执行的结果：(为节省篇幅，把矩阵写成mat2str的形式)（1）A==B [000;010;000]（2）A.*B [2218;322512;98160]（3）A(:)' [147258360]（4）A(1,:)*B(:,3) 139、在MATLAB中，写出下面MATLAB语句执行的结果：（1）c lear,A=ones(2,6)A=（5）B=A（:,end:-1:1)B=211713951231915117310、请编写一段matlab程序，完成以下功能：（1） 生成一个100行，200列的二维随机数组；A=rand(100,200);（2） 找出数组A 中所有大于0.49且小于0.51的元素的单下标；Idx=find(A(:)>0.49&A(:)<0.51)（3） 数组A 中满足（2）中的条件的元素有多少个？11、12、 22cos 0.7sin (1)x y x x ⎡⎤=+⎢⎥+⎣⎦，编写一段matlab 程序，要求如下（1） 在[0,2]π区间，每隔0.01取一x 数值，计算出相应的y 的函数值；x=0:0.01:2*pi;y=(0.7+2*cos(x)./(1+x.^2)).*sin(x)（2）根据MATLAB计算出的数据，找出在[0,2]π内该函数的极小值的坐标。
MATLAB上机答案
X
1
2
3
4
5
6
7
8
9
10
Y
16
32
70
142 260 436 682 1010 1432 1960
>> x=1:10;y=[16 32 70 142 260 436 682 1010 1432 1960];
>> p1=polyfit(x,y,1)
>>
p1 =
p2=polyfit(x,y,2),y2=polyval(p2,9.5)
0.01
0.005
0
-0.005
-0.01
-0.015
0
1
2
3
4
5
6
7
（3）大气压强 p 随高度 x 变化的理论公式为
，为验证这一公式，
测得某地大气压强随高度变化的一组数据如表所示。试用插值法和拟合法进行计算并绘图，
看那种方法较为合理，且总误差最小。
高度/m
0
300
600
1000
1500
2000
压强/Pa
equally spaced points between X1 and X2. 以 X1 为首元素，X2 为末元素平均生成 100 个元素的行向量。
LINSPACE(X1, X2, N) generates N points between X1 and X2. For N < 2, LINSPACE returns X2.
ans =
pi =
ans =
5
0
3.1416
答：3 次执行的结果不一样。exist（）函数是返回变量搜索顺序的一个函数。在第一次
matlab练习题
10.写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值 x=[0,1,0,2,0,3,0,4]; for k=1:8 if x(k)==0 x(k)=k;
else x(k)=2*k+1; end end disp(x); 1 5 3 9 5 13 7 17
11.写出执行以下代码后C，D，E的值 A=[1,2,3;4:6;7:9]; C=[A;[10,11,12]] D=C(1:3,[2 3]) E=C(2,[1 2]) C = 1 4 7 10 2 5 8 11 3 6 9 12
D = 2 5 8 3 6 9
E = 4 >> 12、输入 x,y 的值，并将它们的值互换后输出。 x=input('Input x please:'); y=input('Input y please:'); z=x; x=y; y=z; disp(x); disp(y); 5
13.当 n=100 时，求 y = ∑
15.一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全 部水仙花数。
function function3()
x=[]; for i=100:999 n1=fix(i/100); n2=fix((i-n1*100)/10); n3=i-n1*100-n2*10; if (n1^3+n2^3+n3^3)==i x=[x i]; end end
i =1
n
1 的值。 2i − 1
y=0; n=100; for i=1:n y=y+1/(2*i-1); end y 14.某商场对顾客所购买的商品实行打折销售，标准如下（商品价格用 price 来 表示） ： price<200 没有折扣 200≤price<500 3%折扣 500≤price<1000 5%折扣 1000≤price<2500 8%折扣 2500≤price<5000 10%折扣 5000≤price 14%折扣 输入所售商品的价格，求其实际销售价格。 function function2()
第六讲 MATLAB数值计算
5.2.4 矩阵的秩
求矩阵秩的函数: rank(A)。
例如，求例 5.7 中方程组系数矩阵 D 的秩，命令是：
D=[2, 2, -1, 1; 4, 3, -1, 2; 8, 5, -3, 4; 3, 3, -2, 2]; r=rank(D) r= 4
说明D是一个满秩矩阵。
5.2.5 向量和矩阵的范数
第五讲 MATLAB数值计算
5.1 5.2 5.3 5.4 5.5 5.6 5.7 5.8 5.9 特殊矩阵 矩阵分析 矩阵分解与线性方程组求解 数据处理与多项式计算 傅立叶分析 数值微积分 常微分方程的数值求解 非线性方程的数值求解 稀疏矩阵
5.1 特殊矩阵
5.1.1 对角阵与三角阵
1. 矩阵的对角元素 (1)提取矩阵的对角线元素 设 A 为 m×n 矩阵， diag(A) 函数用于提取矩阵 A主对 角线元素产生一个具有min(m,n)个元素的列向量。 (2)构造对角矩阵
2. 矩阵的伪逆
对奇异方阵和长方阵，求矩阵伪逆的函数是 pinv(A)。
例5.5 求A的伪逆，并将结果送B。 A=[3,1,1,1; 1,3,1,1; 1,1,3,1]; B=pinv(A) 例5.6 求矩阵A的伪逆。 A=[0,0,0; 0,1,0; 0,0,1]; pinv(A)
5.2.3 方阵的行列式
函数A = pascal(n)生成一个n阶的帕斯卡矩阵。
1 1 3 4 6 10 10 20 15 35
1 5 15 35 70
5.2 矩阵分析
5.2.1 矩阵结构变换
1. 矩阵的转置
转置运算符是单撇号(')。
2. 矩阵的旋转
矩阵的旋转利用函数 rot90(A,k) ，功能是将矩阵 A旋 转90º 的k倍，当k为1时可省略。
MATLAB补充习题及答案
2.5 计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。
>> a=[5 3 5;3 7 4;7 9 8];>> b=[2 4 2;6 7 9;8 3 6]; >> a+b ans =7 7 7 9 14 13 15 12 14 2.7 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8]; >> a.*b ans =12 36 3 8 42 402.9 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’; >> X=A\B X =-0.5118 4.0427 1.33182.10 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
>> a=[1 2 3;4 5 6;7 8 9];>> a.^2 ans =1 4 9 16 25 36 49 64 81 >> a^2 ans =30 36 42 66 81 96 102 126 1502.14 角度[]604530=x ，求x 的正弦、余弦、正切和余切。
>> x=[30 45 60];>> x1=x/180*pi; >> sin(x1) ans =0.5000 0.7071 0.8660 >> cos(x1) ans =0.8660 0.7071 0.5000 >> tan(x1) ans =0.5774 1.0000 1.7321 >> cot(x1) ans =1.7321 1.0000 0.57742.17 将矩阵⎥⎦⎤⎢⎣⎡=5724a 、⎥⎦⎤⎢⎣⎡=3817b 和⎥⎦⎤⎢⎣⎡=2695c 组合成两个新矩阵： （1）组合成一个4⨯3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵元素，第三列为按列顺序排列的c 矩阵元素，即 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡237912685574 （2）按照a 、b 、c 的列顺序组合成一个行矢量，即 []296531877254>> a=[4 2;5 7]; >> b=[7 1;8 3]; >> c=[5 9;6 2];% (1)>> d=[a(:) b(:) c(:)] d =4 75 5 86 2 1 97 3 2 % (2)>> e=[a(:);b(:);c(:)]' e =4 5 2 7 7 8 1 3 5 6 9 2 或利用（1）中产生的d >> e=reshape(d,1,12) ans =4 5 2 7 7 8 1 3 5 6 9 2 3.3 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。
MATLAB复习题
MATLAB复习题简答题：0.变量的命名？（1）变量名必须以字母开头，且只能由字母、数字或者下划线3类符号组成，不能含空格和标点符号：（2）变量名区分字母的⼤⼩写。
（3）变量名不能超过63个字符。
（4）关键字不能作为变量名。
1.结构数组的创建？两种（1）对域赋值创建；（2）利⽤函数struct创建。
0.简述使⽤M⽂件与在matlab命令窗⼝中直接输⼊命令有何异同？有何优缺点？在命令窗⼝中输⼊，⼀次只能执⾏⼀⾏命令，脚本⽂件是命令是集中，只需输⼊脚本⽂件名，就可以按顺序执⾏所有命令。
2.简述脚本形式的M⽂件与函数形式的M⽂件的异同？脚本⽂件是命令的集中，只需输⼊脚本⽂件名，就可以⼀次执⾏所有命令。
没有输⼊参数，没有输出参数，⽂件名可以任意取，⼀般需要由其他语句调⽤，⾸⾏语句必须以function开始，有函数名，⽂件名必须与函数名⼀致，函数空间⽣成的变量是局部变量。
3.与其他计算机语⾔相⽐较，MATLAB语⾔突出的特点是什么？MATLAB具有功能强⼤、使⽤⽅便、输⼊简捷、库函数丰富、开放性强等特点。
4.MATLAB系统由那些部分组成？MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语⾔、图形功能和应⽤程序接⼝五个部分组成。
5.安装MATLAB时，在选择组件窗⼝中哪些部分必须勾选，没有勾选的部分以后如何补安装？在安装MATLAB时，安装内容由选择组件窗⼝中个复选框是否被勾选来决定，可以根据⾃⼰的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。
第⼀次安装没有选择的内容在补安装时只需按照安装的过程进⾏，只是在选择组件时只勾选要补装的组件或⼯具箱即可。
6.MATLAB操作桌⾯有⼏个窗⼝？如何使某个窗⼝脱离桌⾯成为独⽴窗⼝？⼜如何将脱离出去的窗⼝重新放置到桌⾯上？在MATLAB操作桌⾯上有五个窗⼝，在每个窗⼝的右上⾓有两个⼩按钮，⼀个是关闭窗⼝的Close按钮，⼀个是可以使窗⼝成为独⽴窗⼝的Undock按钮，点击Undock按钮就可以使该窗⼝脱离桌⾯成为独⽴窗⼝，在独⽴窗⼝的view菜单中选择Dock ……菜单项就可以将独⽴的窗⼝重新防⽌的桌⾯上。
MATLAB总复习题(14151)
diag(diag(A))产生的矩阵是[1 0 0;0 5 0;0 0 9]，若有命令 B=A([3 1],:)，则 B=[7
8 9;1 2 3].
24. MATLAB 的基本功能有 数值计算 ， 符号计算 ， 绘图 ，
编程
、 GUI 开发 。
25.
利用冒号“：”产生矢量
100
， 100
B．在模块旁双击鼠标后出现的编辑框中
4
实现
C．在模块的下边或旁边的编辑框中 D．在模块的属性窗口中
19、以下 4 种语句中 MATLAB 独有的语句是( C )。
A．if 语句
B．for 语句
C．try 语句
D．switch 语句
20、以下 MATLAB 变量名，合法的是( A )。
A．flower2
1) array1(2,:)= 0.3000 8.1000 7.6000 1.8000
2) array1(1:3,[3 3 4])= [4.1000,4.1000,-2.7000;-7.6000,-7.6000,1.8000;0.5000,0.5000,0.4000]
3) array1([1 1],[3 4])=[4.1000,-2.7000;4.1000,-2.7000] 27. 设 a,b,c 和 d 的定义如下，计算如下表达式的值： a=3,b=2,c=9,d=1 1)a*b^2>a*c ans = 0 2)d|b>a ans = 1 3)(d|b)>a ans = 0
选择题
1、产生四维单位矩阵的语句为( B )。
A．ones(4)
B．eye(4)
C．zeros(4)
D．rand(4)
2、求解线性方程组 Ax b ，当 detA 0 ，方程的解是( A )。
MATLAB)课后实验答案-精简版
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MA TLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z e=+(2) 21ln(2z x =+，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦(3) 0.30.330.3sin(0.3)ln,3.0, 2.9,,2.9,3.022aaee a z a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.52. 已知：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 和B1234530166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D 。
(3) 查看MA TLAB 工作空间的使用情况。
4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
(2). 建立一个字符串向量 例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R R S A OS +⎡⎤=⎢⎥⎣⎦。
Matlab练习
Matlab 练习习题1． 设a=1.2,b=-4.6,c=8.0,e=-4.0,计算⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+=d bc e a t ππ22arctan2． 设a=5.67,b=7.811,计算)lg(b a e ba ++3． 已知园半径为15，求其直径、周长和面积。
4． 已知三角形三边a=8.5,b=14.6,c=18.5,求三角形面积2/)(,))()(((c b a s c s b s a s s area ++=---=5． 下列命令执行后，L1、L2、L3、L4的值分别是多少？A=1:9；B=10-A ； L1=A==B ； L2=A<=5； L3=A>3&A<7;L4=find(A>3&A<7)； 习题1． 设矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=321212113A⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=111012111B求：（1）2A+B ；（2）4A 2-3B 2；（3）AB ；（4）BA ；（5）AB-BA2.设三阶矩阵A 、B ，满足A -1BA=6A+BA⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=71000410031A 求矩阵B2． 设（2E-C -1B ）A T =C -1,其中E 是4阶单位矩阵，A T 是4阶矩阵A 的转置，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=1021000032231021B ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=1021000002101021C求矩阵A3． 有一4阶魔方矩阵a ，找出矩阵中大于7的元素，并将它们重新排列成列向量b 。
4． 给定矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=054000031A ,删去整行和整列的0。
5． 设2阶矩阵A 、B 、X ，满足X-2A=B-X,⎥⎦⎤⎢⎣⎡--=2112A⎥⎦⎤⎢⎣⎡--=0220B 求矩阵X6． 求矩阵的主对角元素、逆矩阵、行列式的值、秩、特征值和特征向量⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=163053064A ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=1124111221B7． 分别用矩阵求逆、矩阵除法求方程组的解⎪⎩⎪⎨⎧=+--=+-=+-1074453932z y x z y x z y x 8． 已知多项式P 1(x)=3x+2,P 2(x)=5x 2-x ＋2，P 3(x)=x 2-0.5,求：（1） P(x)=P 1(x)+P 2(x)+P 3(x) （2） P(x)=P 1(x)*P 2(x)*P 3(x) （3） P(x)=0的全部根计算x i =0.2*i,i=0,1,2各点上的P(x i )。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1、利用MATLAB提供的randn函数生成符合正态分布的10×5随机矩阵A,进行如下操作： (1)A各列元素的均值和标准方差。 (2)A的最大元素和最小元素。 (3)求A每行元素的和以及全部元素之和。 (4)分别对A的每列元素按升序、每行元素按降序排序。 A=randn(10,5);disp('各列元素的均值：');mean(A) disp('各列元素的标准方差：');std(A) disp('A的最大元素：');max(max(A)) disp('A的最小元素：');min(min(A)) disp('A每行元素之和：');sum(A,2) disp('全部元素之和：');sum(sum(A)) disp('每列元素按升序：');Y=sort(A) disp('每行元素按降序：');Y=sort(A,2,'descend')
各列元素的均值： ans = -0.1095 0.1282 -0.2646 0.3030 -0.2464 各列元素的标准方差： ans = 0.9264 1.2631 0.8129 0.8842 1.3151 A的最大元素：ans =2.5855 A的最小元素：ans =-1.9330 A每行元素之和： ans = -2.2970 1.2545 0.0661 5.0489 -0.6988 1.1002 -2.9310 -2.0595 -1.6878 0.3112 全部元素之和：ans =-1.8932 每列元素按升序： Y = -1.2141 -1.4916 -1.4224 -1.1658 -1.9330 -1.1135 -1.0891 -1.4023 -0.8045 -1.7947 -0.8637 -1.0616 -0.7648 -0.2437 -1.1480 -0.7697 -0.7423 -0.6156 0.1978 -0.6669 -0.2256 0.0326 -0.1961 0.2157 -0.4390 -0.0068 0.0859 -0.1924 0.2916 -0.0825 0.0774 0.5525 -0.1774 0.6966 0.1049 0.3714 1.1006 0.4882 0.8351 0.1873 1.1174 1.5442 0.7481 1.4193 0.7223 1.5326 2.3505 0.8886 1.5877 2.5855 每行元素按降序： Y = 1.4193 -0.6156 -0.8637 -1.0891 -1.1480 0.7481 0.2916 0.1049 0.0774 0.0326 0.7223 0.5525 0.1978 -0.1924 -1.2141 2.5855 1.5877 1.1006 0.8886 -1.1135 1.5442 -0.0068 -0.6669 -0.7648 -0.8045 1.5326 0.6966 0.1873 0.0859 -1.4023
习 题 六
0.8351 -0.0825 -0.7697 -1.4224 -1.4916 0.4882 0.3714 -0.2437 -0.7423 -1.9330 0.2157 -0.1774 -0.2256 -0.4390 -1.0616 2.3505 1.1174 -0.1961 -1.1658 -1.7947 2、按要求对指定函数进行插值和拟合。
(1)按表6.4用三次样条方法插值计算0~90°范围内整数点的正弦值和0~75°范围内整数点的正切值，然后用5次多项式拟合方法计算相同的函数值，并将两种计算结果进行比较。 表6.4 特殊角的正弦与正切值表 α（度） 0 15 30 45 60 75 90
sinα 0 0.2588 0.5000 0.7071 0.8660 0.9659 1.0000 tanα 0 0.2679 0.5774 1.0000 1.7320 3.7320
a=[0 15 30 45 60 75 90]; A=0:90; sing=[0 0.2588 0.5000 0.7071 0.8660 0.9659 1.0000]; b=[0 15 30 45 60 75]; B=0:75; tang=[0 0.2679 0.5774 1.0000 1.7320 3.7320]; SC=interp1(a,sing,A,'spline'); TC=interp1(b,tang,B,'spline'); P1=polyfit(a,sing,5); SN=polyval(P1,A); P2=polyfit(b,tang,5); TN=polyval(P2,B); subplot(1,2,1); plot(A,SC,'.',A,SN,'-*'); subplot(1,2,2); plot(B,TC,'.',B,TN,'-*');
05010000.20.40.60.811.21.402040608000.511.522.533.54
(2)按表6.5 用3次多项式方法插值计算1~100内特殊值的平方根表 表6.5 1~100内特殊值的平方根表
N 1 4 9 16 25 36 49 64 81 100
N 1 2 3 04 05 06 07 08 09 10
N=[1 4 9 16 25 36 49 64 81 100]; SQRT=[1 2 3 4 5 6 7 8 9 10]; n=1:100; SQ=interp1(N,SQRT,n,'cubic')
SQ = Columns 1 through 9 1.0000 1.3729 1.7125 2.0000 2.2405 2.4551 2.6494 2.8292 3.0000 Columns 10 through 18 3.1636 3.3186 3.4661 3.6069 3.7422 3.8729 4.0000 4.1237 4.2435 Columns 19 through 27 4.3599 4.4730 4.5832 4.6907 4.7958 4.8988 5.0000 5.0993 5.1966 Columns 28 through 36 5.2921 5.3857 5.4777 5.5681 5.6570 5.7446 5.8309 5.9160 6.0000 Columns 37 through 45 6.0829 6.1647 6.2454 6.3249 6.4035 6.4810 6.5577 6.6334 6.7082 Columns 46 through 54 6.7823 6.8556 6.9281 7.0000 7.0712 7.1416 7.2113 7.2804 7.3487 Columns 55 through 63 7.4164 7.4835 7.5500 7.6159 7.6812 7.7459 7.8102 7.8739 7.9372 Columns 64 through 72 8.0000 8.0623 8.1242 8.1855 8.2464 8.3068 8.3668 8.4263 8.4854 Columns 73 through 81 8.5441 8.6024 8.6603 8.7178 8.7749 8.8317 8.8881 8.9442 9.0000 Columns 82 through 90 9.0555 9.1107 9.1655 9.2201 9.2744 9.3284 9.3821 9.4354 9.4884 Columns 91 through 99 9.5412 9.5935 9.6456 9.6973 9.7486 9.7996 9.8502 9.9005 9.9505 Column 100 10.0000 3、已知一组实验数据如表6.6所示。
表6.6 一组实验数据 i 1 2 3 4 5
ix 165 123 150 123 141
iy 187 126 172 125 148
求它的线性拟合曲线。 X=[165 123 150 123 141];Y=[187 126 172 125 148]; P=polyfit(X,Y,1);Y1=polyval(P,X); plot(X,Y,'r*',X,Y1,'b')
120125130135140145150155160165120130140150160170180190
4、已知多项式5.0)( ,25)( ,23)(23221xxPxxxPxxP
, 求：
(1) )()()()(321xPxPxPxP
P1=[3 2];P2=[5 -1 2];P3=[1 0 -0.5]; P1=[0 P1];P12=conv(P1,P2);P3=[0 P3]; P=conv(P12,P3)