习题六(Matlab数值计算)课后习题
数值计算与MATLAB方法课后答案
第一章习题1. 序列满足递推关系，取及试分别计算，从而说明递推公式对于计算是不稳定的。
n1 1 0.01 0.00012 0.01 0.0001 0.0000013 0.0001 0.000001 0.000000014 0.000001 0.0000000110-105 0.00000001 10-10n1 1.000001 0.01 0.0000992 0.01 0.000099 -0.000099013 0.000099 -0.00009901-0.010000994 -0.00009901 -0.01000099-1.00015 -0.01000099-1.0001初始相差不大，而却相差那么远，计算是不稳定的。
2. 取y0=28，按递推公式，去计算y100，若取（五位有效数字），试问计算y100将有多大误差？y100中尚留有几位有效数字？解：每递推一次有误差因此，尚留有二位有效数字。
3．函数，求f(30)的值。
若开方用六位函数表，问求对数时误差有多大？若改用另一等价公式计算，求对数时误差有多大？设z=ln(30-y)，，y*， |E(y)| 10-4z*=ln(30-y*)=ln(0.0167)=-4.09235若改用等价公式设z=-ln(30+y)，，y*， |E(y)|⨯10-4z*=-ln(30+y*)=-ln(59.9833)=-4.094074．下列各数都按有效数字给出，试估计f的绝对误差限和相对误差限。
1)f=sin[(3.14)(2.685)]设f=sin xyx*=3.14, E(x)⨯10-2, y*=2.685, E(y)⨯10-3,sin(x*y*)=0.838147484, cos(x*y*)=-0.545443667⨯(-0.5454) ⨯⨯10-2+3.14(-0.5454) ⨯⨯10-3|⨯10-2⨯10-2|E r(f)| ⨯10-2⨯10-2<10-22)f=(1.56)设f = x y ,x*=1.56, E(x)⨯10-2, y*=3.414, E(y)⨯10-3,⨯⨯⨯10-2⨯⨯⨯10-3|⨯⨯⨯10-2⨯⨯⨯10-3|=0.051|E r(f)| =0.01125．计算，利用下列等式计算，哪一个得到的结果最好，为什么？6．下列各式怎样计算才能减少误差？7. 求方程x2-56x+1=0的二个根，问要使它们具有四位有效数字，至少要取几位有效数字？如果利用伟达定理， 又该取几位有效数字呢？解一:若要取到四位有效数字,如果利用伟达定理,解二:由定理二,欲使x1,x2有四位有效数字,必须使由定理一知，∆至少要取7位有效数字。
matlab课后习题与答案
习题二1.如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
因此，矩阵是MATLAB最基本、最重要的数据对象。
2.设A和B是两个同维同大小的矩阵，问：(1)A*B和A.*B的值是否相等？答：不相等。
(2)A./B和B.\A的值是否相等？答：相等。
(3)A/B和B\A的值是否相等？答：不相等。
(4)A/B和B\A所代表的数学含义是什么？答：A/B等效于B的逆右乘A矩阵，即A*inv(B)，而B\A等效于B矩阵的逆左乘A矩阵，即inv(B)*A。
3.写出完成下列操作的命令。
(1)将矩阵A第2~5行中第1, 3, 5列元素赋给矩阵B。
答：B=A(2:5,1:2:5); 或B=A(2:5,[1 3 5])(2)删除矩阵A的第7号元素。
答：A(7)=[](3)将矩阵A的每个元素值加30。
答：A=A+30;(4)求矩阵A的大小和维数。
答：size(A);ndims(A);(5)将向量t的0元素用机器零来代替。
答：t(find(t==0))=eps;(6)将含有12个元素的向量x转换成34矩阵。
答：reshape(x,3,4);(7)求一个字符串的ASCII码。
答：abs(‘123’); 或double(‘123’);(8)求一个ASCII码所对应的字符。
答：char(49);4. 下列命令执行后，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);答：L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0]L2的值为[1, 1, 1, 1, 1, 0, 0, 0, 0]L3的值为[0, 0, 0, 1, 1, 1, 0, 0, 0]L4的值为[4, 5, 6]5. 已知23100.7780414565532503269.5454 3.14A -⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥-⎣⎦完成下列操作：(1) 取出A 的前3行构成矩阵B ，前两列构成矩阵C ，右下角32⨯子矩阵构成矩阵D ，B 与C 的乘积构成矩阵E 。
完整word版,Matlab课后习题
习题 11. 执行下列指令，观察其运算结果, 理解其意义： (1) [1 2;3 4]+10-2i(2) [1 2; 3 4].*[0.1 0.2; 0.3 0.4] (3) [1 2; 3 4].\[20 10;9 2] (4) [1 2; 3 4].^2 (5) exp([1 2; 3 4]) (6)log([1 10 100]) (7)prod([1 2;3 4])(8)[a,b]=min([10 20;30 40]) (9)abs([1 2;3 4]-pi)(10) [1 2;3 4]>=[4,3;2 1](11)find([10 20;30 40]>=[40,30;20 10])(12) [a,b]=find([10 20;30 40]>=[40,30;20 10]) (提示：a 为行号，b 为列号) (13) all([1 2;3 4]>1) (14) any([1 2;3 4]>1) (15) linspace(3,4,5) (16) A=[1 2;3 4];A(:,2)2. 执行下列指令，观察其运算结果、变量类型和字节数，理解其意义： (1) clear; a=1,b=num2str(a),c=a>0, a= =b, a= =c, b= =c (2) clear; fun='abs(x)',x=-2,eval(fun),double(fun)3. 本金K 以每年n 次，每次p %的增值率(n 与p 的乘积为每年增值额的百分比)增加，当增加到rK 时所花费的时间为)01.01ln(ln p n rT +=(单位：年)用MA TLAB 表达式写出该公式并用下列数据计算：r =2, p =0.5, n =12.4．已知函数f (x )=x 4-2x 在(-2, 2)内有两个根。
取步长h =0.05, 通过计算函数值求得函数的最小值点和两个根的近似解。
matlab教程习题答案
matlab教程习题答案
MATLAB教程习题答案
MATLAB是一种用于数学计算、数据分析和可视化的强大工具。
它被广泛应用于工程、科学和金融领域，因此掌握它的使用对于提高工作效率和解决问题至关重要。
在学习MATLAB的过程中，许多人会遇到各种习题，下面我们将给出一些常见习题的答案，希望能帮助大家更好地理解和掌握MATLAB的使用。
1. 编写一个MATLAB程序，计算并输出1到100之间所有奇数的和。
答案：
```matlab
sum = 0;
for i = 1:100
if mod(i, 2) == 1
sum = sum + i;
end
end
disp(sum);
```
2. 编写一个MATLAB程序，生成一个5x5的随机矩阵，并计算矩阵每一行的平均值。
答案：
```matlab
A = rand(5, 5);
row_means = mean(A, 2);
disp(row_means);
```
3. 编写一个MATLAB程序，求解以下线性方程组的解：
```
2x + 3y = 8
4x - 2y = 2
```
答案：
```matlab
A = [2, 3; 4, -2];
B = [8; 2];
X = A \ B;
disp(X);
```
通过以上习题的答案，我们可以看到MATLAB的强大功能和灵活性。
掌握这些基本的使用技巧，可以帮助我们更好地应用MATLAB解决实际问题，提高工作效率。
希望大家在学习和使用MATLAB的过程中能够不断积累经验，不断提升自己的技能水平。
MATLAB数学实验第二版课后练习题含答案
MATLAB数学实验第二版课后练习题含答案课后练习题MATLAB数学实验第二版的课后练习题如下：第一章课后练习题1.编写MATLAB程序，计算并输出下列公式的结果：y = \\frac{1}{\\sqrt{2\\pi\\sigma^2}} e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}其中，x, $\\mu$, $\\sigma$ 分别由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');mu=input('请输入 mu 的值:');sigma=input('请输入 sigma 的值:');y=1/sqrt(2*pi*sigma^2) *exp(-(x-mu)^2/ (2*sigma^2));fprintf('y = %.2f\', y);2.编写MATLAB程序，求解下列方程的解：4x + y = 11\\\\x + 2y = 7答案如下：A= [4,1;1,2];B= [11;7];X=inv(A) *B;fprintf('x = %.2f, y = %.2f\', X(1), X(2));第二章课后练习题1.编写MATLAB程序，计算下列多项式的值：P(x) = x^4 - 2x^3 + 3x^2 - x + 1其中，x 由用户输入。
要求输出结果精确至小数点后两位。
答案如下：x=input('请输入 x 的值:');y=x^4-2*x^3+3*x^2-x+1;fprintf('P(%.2f) = %.2f\', x, y);2.编写MATLAB程序，绘制下列函数的图像：f(x) = \\begin{cases} x + 1, & x < 0 \\\\ x^2, & 0 \\leq x < 1 \\\\ 2x - 1, & x \\geq 1 \\end{cases}答案如下：x=-2:0.01:2;y1=x+1;y2=x.^2.* ((x>=0) & (x<1));y3=2*x-1;plot(x,y1,x,y2,x,y3);legend('y1 = x + 1','y2 = x^2','y3 = 2x - 1');总结本文提供了《MATLAB数学实验第二版》的部分课后练习题及其答案。
习题六（Matlab数值计算）课后习题
习题六（Matlab数值计算）课后习题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.3151A 的最⼤元素：ans =2.5855A 的最⼩元素：ans =-1.9330A 每⾏元素之和：ans =-2.29701.25450.06615.0489-0.69881.1002-2.9310-2.0595-1.68780.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.08250.0774 0.5525 -0.1774 0.6966 0.10490.3714 1.1006 0.4882 0.8351 0.18731.1174 1.5442 0.7481 1.4193 0.72231.53262.3505 0.8886 1.5877 2.5855每⾏元素按降序：Y =1.4193 -0.6156 -0.8637 -1.0891 -1.14800.7481 0.2916 0.1049 0.0774 0.03260.7223 0.5525 0.1978 -0.1924 -1.21412.5855 1.5877 1.1006 0.8886 -1.11351.5442 -0.0068 -0.6669 -0.7648 -0.80451.5326 0.6966 0.1873 0.0859 -1.40230.8351 -0.0825 -0.7697 -1.4224 -1.49160.4882 0.3714 -0.2437 -0.7423 -1.93300.2157 -0.1774 -0.2256 -0.4390 -1.06162.3505 1.1174 -0.1961 -1.1658 -1.79472、按要求对指定函数进⾏插值和拟合。
matlab教材习题答案
matlab教材习题答案Matlab是一种广泛应用于科学与工程领域的计算机编程语言和环境。
它具备强大的数值计算和数据可视化功能，被广泛用于数据分析、信号处理、图像处理、机器学习等领域。
对于初学者而言，掌握Matlab的基本语法和常用函数是非常重要的，而教材习题则是帮助学生巩固所学知识的重要资源。
本文将为大家提供一些Matlab教材习题的参考答案，以帮助读者更好地学习和应用Matlab。
1. 基本语法练习题1.1 计算并输出1到10的平方for i = 1:10fprintf('%d的平方是：%d\n', i, i^2);end1.2 计算并输出1到10的阶乘for i = 1:10fact = 1;for j = 1:ifact = fact * j;endfprintf('%d的阶乘是：%d\n', i, fact);end2. 数值计算练习题2.1 求解一元二次方程的根a = 1;b = -3;c = 2;delta = b^2 - 4*a*c;x1 = (-b + sqrt(delta))/(2*a);x2 = (-b - sqrt(delta))/(2*a);fprintf('一元二次方程的根为：%f, %f\n', x1, x2);2.2 求解线性方程组的解A = [1 2; 3 4];B = [5; 6];X = inv(A) * B;fprintf('线性方程组的解为：%f, %f\n', X(1), X(2));3. 数据处理练习题3.1 统计一个数组中的最大值、最小值和平均值data = [1, 2, 3, 4, 5];max_value = max(data);min_value = min(data);average_value = mean(data);fprintf('最大值：%f\n最小值：%f\n平均值：%f\n', max_value, min_value, average_value);3.2 对一个矩阵进行排序matrix = [4 2 3; 1 5 6; 9 8 7];sorted_matrix = sort(matrix);fprintf('排序后的矩阵为：\n');disp(sorted_matrix);4. 图像处理练习题4.1 读取并显示一张图片image = imread('image.jpg');imshow(image);4.2 对一张图片进行灰度化处理gray_image = rgb2gray(image);imshow(gray_image);5. 信号处理练习题5.1 生成并绘制正弦信号t = 0:0.01:2*pi;x = sin(t);plot(t, x);5.2 对一段音频信号进行傅里叶变换[y, fs] = audioread('audio.wav');Y = fft(y);plot(abs(Y));通过以上几个例子，我们可以看到Matlab的强大功能和灵活性。
数值计算课后习题答案(全)
习 题 一 解 答1．取3.14，3.15，227，355113作为π的近似值，求各自的绝对误差，相对误差和有效数字的位数。
分析：求绝对误差的方法是按定义直接计算。
求相对误差的一般方法是先求出绝对误差再按定义式计算。
注意，不应先求相对误差再求绝对误差。
有效数字位数可以根据定义来求，即先由绝对误差确定近似数的绝对误差不超过那一位的半个单位，再确定有效数的末位是哪一位，进一步确定有效数字和有效数位。
有了定理2后，可以根据定理2更规范地解答。
根据定理2，首先要将数值转化为科学记数形式，然后解答。
解：（1）绝对误差:e(x)=π－3.14＝3.14159265…－3.14＝0.00159…≈0.0016。
相对误差：3()0.0016()0.51103.14r e x e x x -==≈⨯有效数字：因为π＝3.14159265…=0.314159265…×10，3.14＝0.314×10，m=1。
而π－3.14＝3.14159265…－3.14＝0.00159…所以│π－3.14│＝0.00159…≤0.005=0.5×10－2＝21311101022--⨯=⨯所以，3.14作为π的近似值有3个有效数字。
（2）绝对误差:e(x)=π－3.15＝3.14159265…－3.14＝－0.008407…≈－0.0085。
相对误差：2()0.0085()0.27103.15r e x e x x --==≈-⨯有效数字：因为π＝3.14159265…=0.314159265…×10，3.15＝0.315×10，m=1。
而π－3.15＝3.14159265…－3.15＝－0.008407…所以│π－3.15│＝0.008407……≤0.05=0.5×10－1＝11211101022--⨯=⨯所以，3.15作为π的近似值有2个有效数字。
（3）绝对误差:22() 3.14159265 3.1428571430.0012644930.00137e x π=-=-=-≈-相对误差：3()0.0013()0.4110227r e x e x x--==≈-⨯有效数字：因为π＝3.14159265…=0.314159265…×10， 223.1428571430.3142857143107==⨯，m=1。
- 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.3151A 的最大元素：ans =2.5855A 的最小元素：ans =-1.9330A 每行元素之和：ans =-2.29701.25450.06615.0489-0.69881.1002-2.9310-2.0595-1.68780.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.08250.0774 0.5525 -0.1774 0.6966 0.10490.3714 1.1006 0.4882 0.8351 0.18731.1174 1.5442 0.7481 1.4193 0.72231.53262.3505 0.8886 1.5877 2.5855每行元素按降序：Y =1.4193 -0.6156 -0.8637 -1.0891 -1.14800.7481 0.2916 0.1049 0.0774 0.03260.7223 0.5525 0.1978 -0.1924 -1.21412.5855 1.5877 1.1006 0.8886 -1.11351.5442 -0.0068 -0.6669 -0.7648 -0.80451.5326 0.6966 0.1873 0.0859 -1.40230.8351 -0.0825 -0.7697 -1.4224 -1.49160.4882 0.3714 -0.2437 -0.7423 -1.93300.2157 -0.1774 -0.2256 -0.4390 -1.06162.3505 1.1174 -0.1961 -1.1658 -1.79472、按要求对指定函数进行插值和拟合。
(1)按表6.4用三次样条方法插值计算0~90°范围内整数点的正弦值和0~75°范围内整数点的正切值，然后用5次多项式拟合方法计算相同的函数值，并将两种计算结果进行比较。
表6.4 特殊角的正弦与正切值表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,'-*');内特殊值的平方根表SQRT=[1 2 3 4 5 6 7 8 9 10];n=1:100;SQ=interp1(N,SQRT,n,'cubic')SQ =Columns 1 through 91.0000 1.3729 1.71252.0000 2.2405 2.4551 2.6494 2.82923.0000Columns 10 through 183.1636 3.3186 3.4661 3.6069 3.7422 3.87294.0000 4.12374.2435Columns 19 through 274.3599 4.4730 4.5832 4.6907 4.7958 4.89885.0000 5.09935.1966Columns 28 through 365.2921 5.3857 5.4777 5.5681 5.6570 5.7446 5.8309 5.91606.0000Columns 37 through 456.0829 6.1647 6.2454 6.3249 6.4035 6.4810 6.5577 6.63346.7082Columns 46 through 546.7823 6.8556 6.92817.0000 7.0712 7.1416 7.2113 7.28047.3487Columns 55 through 637.4164 7.4835 7.5500 7.6159 7.6812 7.7459 7.8102 7.87397.9372Columns 64 through 728.0000 8.0623 8.1242 8.1855 8.2464 8.3068 8.3668 8.42638.4854Columns 73 through 818.5441 8.6024 8.6603 8.7178 8.7749 8.8317 8.8881 8.9442 9.0000Columns 82 through 909.0555 9.1107 9.1655 9.2201 9.2744 9.3284 9.3821 9.43549.4884Columns 91 through 999.5412 9.5935 9.6456 9.6973 9.7486 9.7996 9.8502 9.90059.9505Column 10010.00003、已知一组实验数据如表6.6所示。
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')4、已知多项式5.0)( ,25)( ,23)(23221-=+-=+=x x P x x x P x x P , 求：(1) )()()()(321x P x P x P x P =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)P =0 0 15.0000 7.0000 -3.5000 0.5000 -2.0000 -2.0000(2) 的全部根。
0)(=x Px=roots(P)x =0.70710.1000 + 0.6245i0.1000 - 0.6245i-0.7071-0.6667(3) 计算)()10 , ,2 ,1 ,0( 2.0i i x P i i x 各点上的 ==。
xi=0.2*(0:10);PZ=polyval(P,xi)PZ =Columns 1 through 9-2.0000 -2.3920 -2.6112 -1.7024 2.7104 15.0000 42.1120 94.1408 184.9056Columns 10 through 11332.5264 560.00005、求函数在指定区间的极值。
(1) 附近的最小值。
在]1 ,1[23),(22y xy x y x f ++=f=@(x)3*x(1).^2+2*x(1).*x(2)+x(2).^2;[U,fmin]=fminsearch(f,[1,1])U =1.0e-004 *-0.0675 0.1715fmin =1.9920e-010 (2) 内的最大值。
在] ,0[cos sin )(2πx x x f +=f=inline('-sin(x)-cos(x.^2)');fminbnd(f,0,pi)ans =0.73106、求函数在指定点的数值导数。
(1) 2 ,3 ,4 ,6 ,cos sin )(22ππππ=+=x x x x ff=inline('(sin(x)).^2+(cos(x)).^2');dx=diff(f([0 pi/6 pi/4 pi/3 pi/2]))dx =0 0 0 0 (2) 3 ,2 ,1 ,1)(2=+=x x x ff=inline('sqrt(x.^2+1)');dx=diff(f([0 1 2 3]))dx =0.4142 0.8219 0.92627、求数值积分(1) xdx x 5sin sin 05⎰πf=inline('(sin(x)).^5.*sin(5*x)');I=quadl(f,0,pi)I =0.0982(2) dx x x ⎰-++114211 f=inline('(1+x.^2)./(1+x.^4)');I=quadl(f,-1,1)I =2.2214 (3) dx xx x ⎰+π02cos 1sin f=@(x)[x.*sin(x)./(1+(cos(x)).^2)];I=quadl(f,0,pi)I =2.4674 (4) ππ≤≤≤≤Ω+⎰⎰Ωy x dxdy y x 0 ,0: ,)cos(f=@(x,y)[abs(cos(x+y))];I=dblquad(f,0,pi,0,pi)I =6.28328、图并绘制相应的振幅频率做快速傅里叶变换，采样，用从对取已知FFT s t N t e t h t 5~0,64,0,)(=≥=-t=linspace(0,5,64);h=exp(-t);dt=t(2)-t(1);f=1/dt;f=f*(1:64)/64;H=fft(h);plot(f,abs(H),'-*');xlabel('频率');ylabel('振幅');9、分别用矩阵求逆、矩阵除法以及矩阵分解求线性方程组的解。