习题 1
1. 执行下列指令,观察其运算结果, 理解其意义:
(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时所花费的时间为
(单位:年)
用MATLAB表达式写出该公式并用下列数据计算:r=2, p=0.5, n=12.
4.已知函数f(x)=x42x 在(-2, 2)内有两个根。取步长h=0.05, 通过计算函数值求得函数的最小值点和两个根的近似解。(提示:求近似根等价于求函数绝对值的最小值点)
5. (1) 用z=magic(10)得到10阶魔方矩阵;
(2) 求z的各列元素之和;
(3) 求z的对角线元素之和(提示:先用diag(z)提取z的对角线);
(4) 将z的第二列除以;
(5) 将z的第3行元素加到第8行。
6. 先不用MATLAB判断下面语句将显示什么结果?size(B)又得出什么结果?
B1={1:9;' David Beckham '};
B2={180:-10:100; [100,80,75,;77,60,92;67 28 90;100 89 78]};
B=[B1, B2];
B{1,2}(8)
D=cell2struct(B,{'f1','f2'},2);
[a,b]=D.f1
然后用MATLAB验证你的判断。进一步,察看变量类型和字节数,并用Workspace工具栏显示B和D的具体内容。
习题 2
1. 设x为一个长度为n的数组,编程求下列均值和标准差
, n>1
2. 求满足>100的最小m值。
3. 用循环语句形成Fibonacci数列 F1 = F2 =1, Fk = Fk-1 + Fk-2 , k=3,4,…。并验证极限
. (提示:计算至两边误差小于精度 10-8)
4. 分别用for和while循环结构编写程序,求出。并考虑一种避免循环语句的程序设计,比较不同算法的运行时间。
5.假定某天的气温变化记录如下表,试作图描述这一天的气温变化规律。
时刻t(h) | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 |
温度oC(t) | 15o | 14o | 14o | 14o | 14o | 15o | 16o | 18o | 20o | 22o | 23o | 25o | 28o |
时刻t(h) | 13 | 14 | 15 | 16 | 17 | 18 | 如何用matlab将已知点连线 19 | 20 | 21 | 22 | 23 | 24 | |
温度oC(t) | 31o | 32o | 31o | 29o | 27o | 25o | 24o | 22o | 20o | 18o | 17o | 16o | |
6. 作出下列函数图象
(i) 曲线y = x2 sin (x2 - x - 2), -2 x 2 (要求分别使用plot或fplot完成)
(ii) 椭圆x2/4 + y2/9 = 1
(iii) 抛物面z = x2 + y2 , x<3, y<3
(iv) 曲面 z=x4+3x2+y2-2x-2y-2x2y+6, |x|<3, -3<y<13
(v) 空间曲线x=sint, y=cost, z=cos(2t), 0<t<2
(vi) 半球面 x=2sincos, y=2sinsin, z=2cos, 03600, 0900
(vii) 三条曲线合成图y1=sinx, y2=sinxsin(10x), y3= sinx , 0<x<
7.作下列分段函数图
8. 查询trapz的功能和用法:查trapz.m文件所在目录,查看trapz.m的程序结构,查看trapz.m文件所在目录还有哪些文件?
9. 用MATLAB函数表示下列函数,并作图。
10. 已知连续时间Lyapunov方程为
AX+XA’= C
其中A=, C=. 试通过lookfor和help的帮助用MATLAB求解。
习题 3
1. 设a=(1,2,3),b=(2,4,3), 分别计算a./b, a.\b, a/b, a\b, 分析结果的意义。
2. 用矩阵除法解下列线性方程组,并判断解的意义
(1)
(2)
(3)
(4)
3. 求第2题第(4)小题的通解。
4. (人口流动趋势)对城乡人口流动作年度调查,发现有一个稳定的朝向城镇流动的趋势,
每年农村居民的5%移居城镇而城镇居民的1%迁出,现在总人口的20%位于城镇。假如城乡总人口保持不变,并且人口流动的这种趋势继续下去,那么
(1)一年以后住在城镇人口所占比例是多少?两年以后呢?十年以后呢?
(2)很多年以后呢?
(3)如果现在总人口70%位于城镇,很多年以后城镇人口所占比例是多少?
(4)计算转移矩阵的最大特征值及对应的特征向量,与问题(2)(3)有何关系?
5. (经济预测)在某经济年度内,各经济部门的投入产出表如下表3.5(单位:亿元)
消耗部门 | 最后需求 | 总产值 | ||||
工 业 | 农 业 | 第三产业 | ||||
生 产 部 门 | 工 业 | 6 | 2 | 1 | 16 | 25 |
农 业 | 2.25 | 1 | 0.2 | 1.55 | 5 | |
第三产业 | 3 | 0.2 | 1.8 | 15 | 20 | |
假设某经济年度工业,农业及第三产业的最后需求均为17亿元,预测该经济年度工业,农业及第三产业的产出(提示:对于一个特定的经济系统而言,直接消耗矩阵和Leontief矩阵可视作不变)。
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论