您的位置：360文档中心› 习题六(Matlab数值计算)课后习题

习题六(Matlab数值计算)课后习题

合集下载

数值计算与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课后习题与答案

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课后习题

完整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的过程中，许多人会遇到各种习题，下面我们将给出一些常见习题的答案，希望能帮助大家更好地理解和掌握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数学实验第二版课后练习题含答案课后练习题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数值计算）课后习题

习题六（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教材习题的参考答案，以帮助读者更好地学习和应用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。

matlab课后习题答案1到6章

matlab课后习题答案1到6章

欢迎共阅习题二1.如何理解“矩阵是MATLAB 最基本的数据对象”？ 答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。

(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’); 或E 。

答：B=A(1:3,:); C=A(:,1:2); D=A(2:4,3:4);E=B*C;(2) 分别求E<D 、E&D 、E|D 、~E|~D 和find(A>=10&A<25)。

答：E<D=010001⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，E&D=110111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，11⎡⎤⎢⎥答：student(1).id='0001';student(1).name='Tom'; student(1).major='computer';student(1).grade=[89,78,67,90,86,85]; 8.建立单元矩阵B 并回答有关问题。

B{1,1}=1;B{1,2}='Brenden';B{2,1}=reshape(1:9,3,3); B{2,2}={12,34,2;54,21,3;4 ,23,67};(1)size(B)和ndims(B)的值分别是多少？答：size(B) 的值为2, 2。

ndims(B) 的值为2。

(2)B(2)和B(4)的值分别是(2)建立5×6随机矩阵A，其元素为[100,200]范围内的随机整数。

matlab课后习题答案刘

matlab课后习题答案刘

matlab课后习题答案刘MATLAB课后习题答案刘在本学期的MATLAB课程中，我们学习了MATLAB的基础知识和一些高级功能。

为了帮助同学们更好地复习和巩固所学知识，以下是针对课后习题的一些参考答案。

1. 向量和矩阵的基本操作- 创建一个3x3的单位矩阵。

```matlabA = eye(3);```- 计算矩阵A和B的乘积，其中B是A的转置。

```matlabB = A';C = A * B;```2. 条件语句和循环- 编写一个循环，打印出1到10的所有奇数。

```matlabfor i = 1:2:10disp(i);end```- 编写一个条件语句，检查一个数是否为素数。

```matlabfunction isPrime(n)if n <= 1isPrime = false;elsefor i = 2:sqrt(n)if mod(n, i) == 0isPrime = false;return;endendisPrime = true;endend```3. 函数和脚本- 编写一个函数，计算给定半径的圆的面积。

```matlabfunction area = circleArea(radius)area = pi * radius^2;end```- 编写一个脚本来调用上述函数，并打印出半径为5的圆的面积。

```matlabradius = 5;disp(['Area of circle with radius ', num2str(radius), ' is ', num2str(circleArea(radius))]);```4. 数据结构和文件I/O- 读取一个文本文件中的数据到一个cell数组。

```matlabfilename = 'data.txt';fileID = fopen(filename, 'r');data = textscan(fileID, '%s');fclose(fileID);```- 将一个矩阵写入到一个CSV文件中。

Matlab 课后习题集

Matlab 课后习题集

Matlab程序设计课程实验目录Matlab操作实验（第1章）实验指导思考练习矩阵运算实验（第2章）实验指导思考练习程序设计实验（第3章）实验指导思考练习高级绘图实验（第4章）实验指导思考练习数据分析与多项式计算实验（第5章）实验指导思考练习解方程与最优化问题求解实验（第6章）实验指导思考练习数值运算实验（第7章）实验指导思考练习符号实验（第8章）实验指导思考练习图形句柄实验（第9章）实验指导思考练习图形用户界面设计实验（第10章）实验指导思考练习Simulink动态仿真实验（第11章）实验指导思考练习外部接口技术实验（第12章）实验指导思考练习1、Matlab操作实验一、实验内容1. 第一章实验内容1：path(path,’系统盘名:\matlab文件名\要添加的文件名’)如：path(path,’c:\matlab701\my work’)实验内容-2例1-1：x=[0:0.5:360]*pi/180;plot(x,sin(x),x,cos(x));例1-2p=[3,7,9,0,-23];x=roots(p)例1-3quad(‘x.*log(1+x)’,0,1)例1-4：实验内容3Help 函数名如：help inv help plot等实验内容4（1）x=0:pi/10:2*pi;y=sin(x);二、思考练习思考题1启动：有三种常见方法，（1）、在Windows桌面，单机任务栏上的‚开始‛按钮，选择‚所有程序‛菜单项，然后选择MATLAB程序组中的MATLAB （版本号）程序选项，即可启动。

（2）、在MATLAB的安装路径中找到MATLAB系统启动程序matlab.exe，然后运行它。

（3）、利用建立快捷方式功能，将MATLAB系统启动程序以快捷方式的形式放在Windows桌面上，以后只要在桌面上双击该图标即可启动MATLAB。

关闭：常见的也有3种方法，（1）、在MATLAB主窗口file菜单中选择exit MATLAB命令。

matlab课后习题答案(1-9章)要点

matlab课后习题答案(1-9章)要点

1 数字1.5e2，1.5e3 中的哪个与1500相同吗？1.5e32 请指出如下5个变量名中，哪些是合法的？abcd-2xyz_3 3chan a 变量ABCDefgh 2、5是合法的。

3 在MATLAB 环境中，比1大的最小数是多少？ 1+eps4 设 a = -8 , 运行以下三条指令，问运行结果相同吗？为什么？w1=a^(2/3) w2=(a^2)^(1/3) w3=(a^(1/3))^2w1 = -2.0000 + 3.4641i ;w2 = 4.0000 ;w3 =-2.0000 + 3.4641i 5 指令clear, clf, clc 各有什么用处？clear 清除工作空间中所有的变量。

clf 清除当前图形。

clc 清除命令窗口中所有显示。

第二章1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”符号对象？3/7+0.1双; sym(3/7+0.1)符; sym('3/7+0.1') 符;; vpa(sym(3/7+0.1)) 符;2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') symvar(sym('sin(w*t)'),1) w a z3 （1）试写出求三阶方程05.443=-x 正实根的程序。

注意：只要正实根，不要出现其他根。

（2）试求二阶方程022=+-a ax x 在0>a 时的根。

（1）reset(symengine)syms x positive solve(x^3-44.5) ans =(2^(2/3)*89^(1/3))/2（2）求五阶方程022=+-a ax x 的实根 syms a positive %注意：关于x 的假设没有去除 solve(x^2-a*x+a^2)Warning: Explicit solution could not be found. > In solve at 83 ans =[ empty sym ]syms x clear syms a positivesolve(x^2-a*x+a^2) ans =a/2 + (3^(1/2)*a*i)/2 a/2 - (3^(1/2)*a*i)/24 观察一个数（在此用@记述）在以下四条不同指令作用下的异同。

MATLAB习题及答案

MATLAB习题及答案

填空题1. MATLAB于1984年由美国Mathworks公司推出，其后每年更新（两次。

2. MATLAB是一种以（矩阵）运算为基础的交互式程序设计语言。

3. MATLAB具有卓越的数值计算能力和符号计算、文字处理、可视化建模仿真和实时控制等众多功能，其每个变量代表一个（矩阵），每个元素都看作（复数）。

4.通过命令（help）、（lookfor），可以查找所有命令或函数的使用方法。

5.执行语句a=1:2:10，得到的一维数组是（1 3 5 7 9）.6.执行语句b=linspace（1,10,10）后，一维数组b包含（10）个元素，最大值是10）7.函数rem（）的功能是取（余）数。

8.若p=[1 0 0;1 1 0],则p|〜p=（[1 1 1;1 1 1]）.（注：填空时请用本题的p的方式表示结果）9.若p=[1 0 0;1 1 0],则all（p）=（[1 0 0]）.10.矩阵的加减运算，要求相加减的矩阵阶数相同。

若A=[1 2 3 4；2 3 1 8]，则执行语句：[n,m]=size（A），则n=（2 ）,m=（4 ）.11.对于一维矩阵，求其长度的函数是（length（））.12.数组和数组之间的运算，尤其是对于乘除运算和乘方运算，如果采用点方式进行计算，表明是数组的（元素）之间的运算关系。

13.求矩阵运算A*B时，要求在维度上，A的（列）数与B的（行）数相等。

二、判断题1.MATLAB只有一种数据类型，一种标准的输入输出语句，不需编译，可直接运行。

（对2.MATLAB的特殊常量是一些预选定义好的数值变量。

（对3.MATLAB变量名不区分大小写。

（错4.i是特殊常量。

（对5.NAN是非数。

（对6.MATLAB中所有的变量都表示一个矩阵或一个向量。

（对7.MATLAB中变量不需要先定义后使用，会自动根据实际赋值的类型对变量类型进行定义。

（对8.clc命令可以从内存中删除一个、多个和所有变量。

matlab课后习题参考答案

matlab课后习题参考答案

Matlab课后习题部分参考答案习题二2．student=struct('name',{'zhang','wang','li'},'age',{18,21,[] },'email',{['zh ang@','zhang@'],'',''})5、a=[1 0 0;1 1 0;0 0 1];b=[2 3 4;5 6 7;8 9 10];c=[a b]d=[a;b]c =1 0 023 41 1 0 5 6 70 0 1 8 9 10d =1 0 01 1 00 0 12 3 45 6 78 9 1010、reshape(c,2,9)1 0 1 0 1 5 3 9 71 0 0 02 8 6 4 10 reshape(d,2,9)ans =1 0 5 0 0 6 0 1 71 2 8 1 3 9 0 4 10 (注意：重新排列矩阵，是将原来的矩阵按列排序)11、a.’12、a+b a.*ba\b (即inv(a)*b，考试时，要求能自己手算出a的逆)ans =2 3 43 3 38 9 1013、c=4*ones(3);a>=cb>=c14．strcat('The picture is ','very good')a=['The picture is ','very good']18．double('very good')1．norm（A），det（A），rank（A）2．b=[1 1 1 1 1];c=b';a=[17,24,1,8,50;23,5,7,14,49;4,6,13,20,43;10,12,19,21,62;11,18,25,2,56]; x=a\c习题四1、a=14*rand(1,10)-5;for i=1:9max=i;for j=i+1:10if a(j)>a(max)max=j;endendtemp=a(i);a(i)=a(max);a(max)=temp;enda将上述脚本保存为script41.m2、function y=function42(n)3*randn(1,n)+3;for i=1:n-1max=i;for j=i+1:nif a(j)>a(max)max=j;endendtemp=a(i);a(i)=a(max);a(max)=temp;endy=a;将上述函数保存为function42.m 3、r=input(‘please input 1/2: ’);if r= =1script41;elsen=input(‘please input the length: ’);y=funtion42(n);end4．function f=function1(x,y)if y==1f=sin(x);elseif y==2f=cos(x);elsef=sin(x).*cos(x);end在命令窗口调用函数f=function1(1,3) 5．function result=function3(x,n)result=0;for i=1:nresult=result+sin(x.*i)+(-1).^i.*cos(i.*x); endfunction result=function4(n,t)result=0;A=[1 2 3;0 1 2;0 0 1];for i=1:nresult=result+A*i.*exp(A*i.*t);end习题五3．x1=-2:0.01:2;x2=-2:0.01:2;y1=x1.*sin(x2);y2=x2.*cos(x1);plot3(x1,x2,y 1,'-d',x1,x2,y2,'-d')14、x=8*rand(1,100)-2;y=reshape(x,10,10);save mydata.mat yclearclcload mydatak=y>ones(10);totel=0;for i=1:100if k(i)==1totel=totel+1;endend15、x=rand(5); %生成5*5均匀分布的随机矩阵（该处只作举例，未按题目要求）fid=fopen(‘text.txt’,’w’);count=fwrite(fid,x,’int32’); closestatus=fclose(fid);清除内存，关闭所有窗口fid=fopen(‘text.txt’,’r’);x=fread(fid,[5,5],’int32’); closestatus=fclose(fid); inv(x)16．随便生成一个矩阵，x1=-1:0.2:0.8;x2=-1:0.2:0.8;y1=x1.*sin(x2)；y2=x2.*cos(x1)；a=reshape([y1;y2],10,10)；fid=fopen('table.txt','w'); fprintf (fid,‘%f’，a); fclose(fid) ;清除内存，关闭所有窗口a=fscanf('table.txt','%f') ; exp(a) ;fclose(fid) ;习题六1.c=ploy2str(A)；B=[2,0,0,1,3,5]；2.x=1:10;B=ployval(A，x)；。

matlab习题6答案

matlab习题6答案

matlab习题6答案Matlab习题6答案Matlab是一种广泛应用于科学计算和工程领域的高级编程语言和环境。

它以其强大的数值计算和数据可视化功能而闻名于世。

在Matlab的学习过程中，习题是提高编程技能和理解语言特性的重要方式。

在本文中，我将为大家提供Matlab习题6的答案，帮助大家更好地理解和掌握Matlab编程。

习题1：编写一个Matlab函数，计算给定向量的均值和标准差。

函数的输入参数为一个向量，输出为该向量的均值和标准差。

解答：```matlabfunction [mean_value, std_value] = calculate_stats(vector)mean_value = mean(vector);std_value = std(vector);end```习题2：编写一个Matlab函数，计算给定矩阵的特征值和特征向量。

函数的输入参数为一个矩阵，输出为该矩阵的特征值和特征向量。

解答：```matlabfunction [eigenvalues, eigenvectors] = calculate_eigen(matrix)[eigenvectors, eigenvalues] = eig(matrix);end```习题3：编写一个Matlab函数，计算给定矩阵的逆矩阵。

函数的输入参数为一个矩阵，输出为该矩阵的逆矩阵。

解答：```matlabfunction inverse_matrix = calculate_inverse(matrix)inverse_matrix = inv(matrix);end```习题4：编写一个Matlab函数，实现矩阵的转置操作。

函数的输入参数为一个矩阵，输出为该矩阵的转置矩阵。

解答：```matlabfunction transpose_matrix = calculate_transpose(matrix)transpose_matrix = matrix';end```习题5：编写一个Matlab函数，计算给定向量的累积和。

matlab课后习题答案

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 。

数值计算课后答案6

数值计算课后答案6

习 题 六 解 答1、在区间［0，1］上用欧拉法求解下列的初值问题，取步长h=0.1。

(1)210(1)(0)2y y y '⎧=--⎨=⎩(2)sin (0)0x y x e y -'⎧=+⎨=⎩解：取h=0.1，本初值问题的欧拉公式具体形式为21(1)(0,1,2,)n n n y y y n +=--=由初值y 0=y(0)=2出发计算，所得数值结果如下： x 0=0，y 0=2；x 1=0.1，2100(1)211y y y =--=-= x 2=0.2，2211(1)101y y y =--=-=指出：可以看出，实际上求出的所有数值解都是1。

2、用欧拉法和改进的欧拉法(预测－校正法)求解初值问题，取步长h=0.1。

22(00.5)(0)1y x y x y '⎧=-≤≤⎨=⎩ 解：由预测校正公式11(,)[(,)(,)]2n n n n n n n n y hf x y hy f x y f x y ++⎧=+⎪⎨=++⎪⎩n+1n+1y y ， 取h=0.1，本初值问题的预测－校正公式的具体形式为122210.1(2)0.05[(2)(2)]nn n n n n n n y x y y x y x y ++⎧=+⨯-⎪⎨=+-+-⎪⎩n+1n+1y y 由初值y 0=y(0)=1出发计算，所得数值结果如下： x 0=0，y 0=1； x 1=0.1，2000220001120.1(2)0.8,0.05[(2)(2)]10.05[(02)(0.120.8]0.82y x y y x y x y =+-==+-+-=+-+-⨯=11y y3、试导出解一阶常微分方程初值问题000(,)()()y f x y x a x b y x y '==≤≤⎧⎨=⎩的隐式欧拉格式111(,)(0,1,2,)n n n n y y hf x y n +++=+=并估计其局部截断误差。

matlab课后练习习题及答案详解

matlab课后练习习题及答案详解

matlab课后习题及答案详解第1章MATLAB概论与其余计算机语言对比较，MATLAB语言突出的特色是什么？MATLAB拥有功能强盛、使用方便、输入简捷、库函数丰富、开放性强等特色。

MATLAB系统由那些部分构成？MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。

安装MATLAB时，在选择组件窗口中哪些部分一定勾选，没有勾选的部分此后怎样补安装？在安装MATLAB时，安装内容由选择组件窗口中个复选框能否被勾选来决定，能够依据自己的需要选择安装内容，但基本平台（即MATLAB选项）一定安装。

第一次安装没有选择的内容在补安装时只要依据安装的过程进行，不过在选择组件时只勾选要补装的组件或工具箱即可。

MATLAB操作桌面有几个窗口？怎样使某个窗口离开桌面成为独立窗口？又怎样将离开出去的窗口从头搁置到桌面上？在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是封闭窗口的Close按钮，一个是能够使窗口成为独立窗口的Undock 按钮，点击Undock按钮就能够使该窗口离开桌面成为独立窗口，在独立窗口的view菜单中选择Dock,,菜单项就能够将独立的窗口从头防备的桌面上。

怎样启动M文件编写/调试器？在操作桌面上选择“成立新文件”或“翻开文件”操作时，M文件编写/调试器将被启动。

在命令窗口中键入edit命令时也能够启动M文件编写/调试器。

储存在工作空间中的数组能编写吗？怎样操作？储存在工作空间的数组能够经过数组编写器进行编写：在工作空间阅读器中双击要编写的数组名翻开数组编写器，再选中要改正的数据单元，输入改正内容即可。

命令历史窗口除了能够察看前方键入的命令外，还有什么用途？命令历史窗口除了用于查问从前键入的命令外，还能够直接履行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。

怎样设置目前目录和搜寻路径，在目前目录上的文件和在搜寻路径上的文件有什么差别？目前目录能够在目前目录阅读器窗口左上方的输入栏中设置，搜寻路径能够经过选择操作桌面的file菜单中的SetPath菜单项来达成。

MATLAB数值计算-习题

MATLAB数值计算-习题

1. 用函数 roots 求方程的根roots([1 -1 -1])2. ，在n 个节点 (n 不要太大，如取5~11)上用分段线性、三次方、样条插值方法，计算m 个插值点(m 可取50~100)的函数值。

(注， n 取 10，m 取 100)x=linspace(0,2*pi,10);y=sin(x);xi=linspace(0,2*pi,100);y1=interp1(x,y,xi);y2=interp1(x,y,xi,'spline');y3=interp1(x,y,xi,'cublic');3.测得某地大气压强随高度变化的一组数据如表3-11 所示，试用插值法和拟合法估量高度为0，100，200，300，......，2000 米时的大气压强值。

表3-11 某地大气压强随高度变化数据高度 /m 0 300 600 1000 1500 2000压强 /Pax=[0 300 600 1000 1500 2000];y=[];xi=0:100:2000;y1=interp1(x,y,xi,'spline');p=polyfit(x,y,3);y2=polyval(p,xi);4. 利用梯形法和辛普森法求定积分的值：梯形法：x=linspace(-3,3,200);y=exp(-x.^2/2)/(2*pi);I1=trapz(x,y)辛普森法：I2=quad('exp(-x.^2/2)/(2*pi)',-3,3)或许：t='exp(-x.^2/2)/(2*pi)';I2=quad(t,-3,3)5.分别用矩形法、梯形法、辛普森法和牛顿 -科茨 4 种方法来近似计算定积分。

矩形法：x=linspace(0,1,100);y=x./(x.^2+4);t=cumsum(y*(1/99));t1=t(100)梯形法：x=linspace(0,1,100);y=x./(x.^2+4);t2=trapz(x,y)辛普森法t3=quad('x./(x.^2+4)',0,1)牛顿 -科茨t4=quadl('x./(x.^2+4)',0,1)6.解以下方程组：a=[5 1 2 1;2 5 1 1;1 2 10 2;1 2 2 10];b=[9; 9; 15; 15]x=a\b9.利用二三阶龙格－库塔方法来求解以下初值问题：先定义函数function f=exe9(x,y)f=2*x .*y在命令窗口输入[x1,y1]=ode23('exe9',[0::],1)。

matlab课后习题答案刘

matlab课后习题答案刘

matlab课后习题答案刘
《MATLAB课后习题答案刘》
MATLAB是一种强大的数学软件，广泛应用于工程、科学和数学领域。

在学习MATLAB的过程中，课后习题是巩固知识、提高能力的重要途径。

为了帮助学生更好地掌握MATLAB知识，我们特为大家整理了一份MATLAB课后习题答案刘，希望对大家有所帮助。

1. 问题一：编写一个MATLAB程序，计算1到100之间所有偶数的和。

解答：可以使用for循环结合if语句来实现这个问题。

具体代码如下：
```matlab
sum = 0;
for i = 1:100
if mod(i, 2) == 0
sum = sum + i;
end
end
disp(sum);
```
2. 问题二：编写一个MATLAB程序，计算给定数组的平均值。

解答：可以使用mean函数来计算数组的平均值。

具体代码如下：
```matlab
A = [1, 2, 3, 4, 5];
avg = mean(A);
disp(avg);
```
3. 问题三：编写一个MATLAB程序，绘制正弦函数的图像。

解答：可以使用plot函数来绘制正弦函数的图像。

具体代码如下：
```matlab
x = 0:0.01:2*pi;
y = sin(x);
plot(x, y);
```
通过以上几个例子，我们可以看到MATLAB的强大之处。

希望大家在学习MATLAB的过程中能够多加练习，不断提高自己的编程能力。

如果还有其他问题，欢迎随时向我们咨询。

祝大家学习进步！。

  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 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、分别用矩阵求逆、矩阵除法以及矩阵分解求线性方程组的解。

相关文档
最新文档