《MATLAB程序设计与应用》实验指导书
实验一 matlab 集成环境使用与运算基础
1，先求下列表达式的值，然后显示matlab 工作空间的使用情况并保存全部变量。 （1）0
12
2sin 851z e
=+
程序：
.>> z1=2*sin(85*pi/180)/(1+exp(2)) 结果： z1 =
0.2375
（2）22121(0.45
52
i z In x x +??
=+
=?
?-??
其中
程序：
>> x=[2,1+2*i;-0.45,5]; >> z2=0.5*log(x+sqrt(1+x*x)) 结果： z2 =
0.7114 - 0.0253i 0.8968 + 0.3658i 0.2139 + 0.9343i 1.1541 - 0.0044i
（3）0.3,9.2,8.2,...,8.2,9.2,0.3,2
3.0)3.0sin(2
3.03.03---=+++-=
-a a In
a e e
z a
a
提示：利用冒号表达式生成a 向量，求各点函数值时用点乘运算。 程序：
>> a=-3.0:0.1:30;
>> z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin((a+0.3)*pi/180)+log((0.3+a)/2) 结果： z3 =
1.0e+003 *
Columns 1 through 4
0.0003 + 0.0031i 0.0003 + 0.0031i 0.0003 + 0.0031i 0.0002 + 0.0031i
Columns 5 through 8
0.0002 + 0.0031i 0.0001 + 0.0031i 0.0001 + 0.0031i 0.0000 + 0.0031i Columns 9 through 12
-0.0000 + 0.0031i -0.0001 + 0.0031i -0.0001 + 0.0031i -0.0002 + 0.0031i Columns 13 through 16
-0.0003 + 0.0031i -0.0003 + 0.0031i -0.0004 + 0.0031i -0.0005 + 0.0031i Columns 17 through 20
-0.0006 + 0.0031i -0.0007 + 0.0031i -0.0008 + 0.0031i -0.0009 + 0.0031i Columns 21 through 24
-0.0010 + 0.0031i -0.0012 + 0.0031i -0.0014 + 0.0031i -0.0016 + 0.0031i Columns 25 through 28
-0.0019 + 0.0031i -0.0023 + 0.0031i -0.0030 + 0.0031i -0.0370 Columns 29 through 32
-0.0030 -0.0023 -0.0019 -0.0016 Columns 33 through 36
-0.0014 -0.0012 -0.0010 -0.0009 Columns 37 through 40
-0.0008 -0.0007 -0.0006 -0.0005 Columns 41 through 44
-0.0004 -0.0003 -0.0003 -0.0002 Columns 45 through 48
-0.0001 -0.0001 -0.0000 0.0000 Columns 49 through 52
0.0001 0.0001 0.0002 0.0002 Columns 53 through 56
0.0003 0.0003 0.0003 0.0004 Columns 57 through 60
0.0004 0.0005 0.0005 0.0005 Columns 61 through 64
0.0006 0.0006 0.0006 0.0007 Columns 65 through 68
0.0007 0.0007 0.0008 0.0008 Columns 69 through 72
0.0008 0.0008 0.0009 0.0009 Columns 73 through 76
0.0009 0.0010 0.0010 0.0010 Columns 77 through 80
0.0011 0.0011 0.0011 0.0011 Columns 81 through 84
0.0012 0.0012 0.0012 0.0013 Columns 85 through 88
0.0013 0.0013 0.0013 0.0014 Columns 89 through 92
0.0014 0.0014 0.0015 0.0015 Columns 93 through 96
0.0015 0.0016 0.0016 0.0016 Columns 97 through 100
0.0017 0.0017 0.0017 0.0018 Columns 101 through 104
0.0018 0.0018 0.0019 0.0019 Columns 105 through 108
0.0020 0.0020 0.0020 0.0021 Columns 109 through 112
0.0021 0.0022 0.0022 0.0023 Columns 113 through 116
0.0023 0.0024 0.0024 0.0025 Columns 117 through 120
0.0025 0.0026 0.0026 0.0027 Columns 121 through 124
0.0027 0.0028 0.0029 0.0029 Columns 125 through 128
0.0030 0.0031 0.0031 0.0032 Columns 129 through 132
0.0033 0.0034 0.0034 0.0035 Columns 133 through 136
0.0036 0.0037 0.0038 0.0039 Columns 137 through 140
0.0040 0.0041 0.0042 0.0043 Columns 141 through 144
0.0044 0.0045 0.0046 0.0047 Columns 145 through 148
0.0049 0.0050 0.0051 0.0053 Columns 149 through 152
0.0054 0.0056 0.0057 0.0059 Columns 153 through 156
0.0060 0.0062 0.0064 0.0066 Columns 157 through 160
0.0068 0.0069 0.0071 0.0074 Columns 161 through 164
0.0076 0.0078 0.0080 0.0083 Columns 165 through 168
0.0085 0.0088 0.0090 0.0093 Columns 169 through 172
0.0096 0.0099 0.0102 0.0105 Columns 173 through 176
0.0108 0.0112 0.0115 0.0119 Columns 177 through 180
0.0123 0.0127 0.0131 0.0135 Columns 181 through 184
0.0139 0.0144 0.0148 0.0153 Columns 185 through 188
0.0158 0.0163 0.0168 0.0174 Columns 189 through 192
0.0180 0.0185 0.0191 0.0198 Columns 193 through 196
0.0204 0.0211 0.0218 0.0225 Columns 197 through 200
0.0233 0.0241 0.0249 0.0257 Columns 201 through 204
0.0265 0.0274 0.0284 0.0293 Columns 205 through 208
0.0303 0.0313 0.0324 0.0335 Columns 209 through 212
0.0346 0.0358 0.0370 0.0382 Columns 213 through 216
0.0395 0.0409 0.0423 0.0437 Columns 217 through 220
0.0452 0.0467 0.0483 0.0500 Columns 221 through 224
0.0517 0.0534 0.0552 0.0571 Columns 225 through 228
0.0591 0.0611 0.0632 0.0654 Columns 229 through 232
0.0676 0.0699 0.0723 0.0748 Columns 233 through 236
0.0773 0.0800 0.0827 0.0856 Columns 237 through 240
0.0885 0.0915 0.0947 0.0979 Columns 241 through 244
0.1013 0.1047 0.1083 0.1121 Columns 245 through 248
0.1159 0.1199 0.1240 0.1282 Columns 249 through 252
0.1326 0.1372 0.1419 0.1467 Columns 253 through 256
0.1518 0.1570 0.1624 0.1679 Columns 257 through 260
0.1737 0.1796 0.1858 0.1921 Columns 261 through 264
0.1987 0.2055 0.2125 0.2198 Columns 265 through 268
0.2273 0.2351 0.2431 0.2514 Columns 269 through 272
0.2600 0.2689 0.2781 0.2876 Columns 273 through 276
0.2974 0.3076 0.3180 0.3289 Columns 277 through 280
0.3401 0.3517 0.3637 0.3761 Columns 281 through 284
0.3889 0.4021 0.4158 0.4299 Columns 285 through 288
0.4446 0.4597 0.4753 0.4915 Columns 289 through 292
0.5082 0.5254 0.5433 0.5617 Columns 293 through 296
0.5807 0.6004 0.6208 0.6418 Columns 297 through 300
0.6636 0.6861 0.7093 0.7333 Columns 301 through 304
0.7581 0.7838 0.8103 0.8377 Columns 305 through 308
0.8660 0.8952 0.9254 0.9567 Columns 309 through 312
0.9890 1.0223 1.0568 1.0924 Columns 313 through 316
1.1292 1.1673 1.2066 1.2472
Columns 317 through 320
1.2892 1.3326 1.3774 1.4237
Columns 321 through 324
1.4715 1.5210 1.5721 1.6249
Columns 325 through 328
1.6794 1.7357 1.7940 1.8541
Columns 329 through 331
1.9163 1.9805
2.0468
（4）??
?
??=<≤<≤<≤+--=5.2:5.0:0,322110,121,222
4t t t t t t t t z 其中
提示：用逻辑表达式求分段函数值。 程序：
>> t=0:0.5:2.5;
>> z1=(t>0&t<=1).*t.^2; >> z2=(t>=1&t<2).*(t.^2-1); >> z3=(t>=2&t<3).*(t.^2-2*t+1); >> z=z1+z2+z3 z =
0 0.2500 1.0000 1.2500 1.0000 2.2500
2，已知 123441
3134
787,2033
65
73
2
7A 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] 程序：
（1）
>> A=[12,34,-4;34,7,87;3,65,7];
>> B=[1,3,-1;2,0,3;3,-2,7];
>> I=eye(3);
>> D=A+6*B
结果：
D =
18 52 -10
46 7 105
21 53 49
程序：
>>D =A-B+I
结果：
D=
12 31 -3
32 8 84
0 67 1
（2）
程序：
>> D=A*B
结果：
D =
68 44 62
309 -72 596
154 -5 241
程序：
>> D=A.*B
结果：
D =
12 102 4
68 0 261
9 -130 49
（3）
>> D=A^3
D =
37226 233824 48604
247370 149188 600766
78688 454142 118820 >> D=A.^3
D =
1728 39304 -64
39304 343 658503
27 274625 343
（4）
>> D=A/B
D =
16.4000 -13.6000 7.6000
35.8000 -76.2000 50.2000
67.0000 -134.0000 68.0000
>> D=B\A
D =
109.4000 -131.2000 322.8000
-53.0000 85.0000 -171.0000
-61.6000 89.8000 -186.2000
（4）
>> D=[A,B]
D =
12 34 -4 1 3 -1
34 7 87 2 0 3
3 65 7 3 -2 7 >> D=[A([1,3],:);B^2]
D =
12 34 -4 3 65 7 4 5 1 11 0 19 20 -5 40
3，设有矩阵A 和B
123453
0166789101769A ,11
12131415023416171819209
7021
22
23
24
254
13
11B ????
????-????????==-????????????????
（1） 求它们的乘积C ；
（2） 将矩阵C 的右下角3*2子矩阵赋给D ； （3） 察看matlab 工作空间的使用情况。 程序： （1）
>> A=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,24,25];
>> B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11]; >> C=A*B 结果：
C =
93 150 77 258 335 237 423 520 397 588 705 557 753 890 717 （2）
>> D=C(3:5,2:3) 结果： D =
520 397 705 557 890 717
4，完成下列操作：
（1）求[100,999]之间能被21整除的数的个数。
提示：先利用冒号表达式，再利用find和length函数。（2）建立一个字符串向量，删除其中的大写字母。
提示：利用find函数和空矩阵
程序：
（1）
>> A=100:1:999;
>> B=rem(A,21)==0;
>> C=find(B==1);
>> D=length(C)
结果：
D =
43
（2）
>> x=input('请输入字符串向量:','s');
请输入字符串向量:cToOmPpuIter
>> y=find(x>='A'&x<='Z');
>> x(y)=[];
>> x
x =
computer
实验二 matlab 矩阵分析与处理
1， 设有分块矩阵??
?
???=2*23
*22*33
*3S O R E
A ，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证??
?
?
??+=2
2S O
RS R E
A 。
>> E=eye(3); >> R=rand(3,2); >> O=zeros(2,3); >> S=diag(1:2); >> A=[E,R;O,S] A =
1.0000 0 0 0.9501 0.4860 0 1.0000 0 0.2311 0.8913 0 0 1.0000 0.6068 0.7621 0 0 0 1.0000 0 0 0 0 0
2.0000
>> A^2
ans =
1.0000 0 0 1.9003 1.4579 0 1.0000 0 0.4623
2.6739 0 0 1.0000 1.2137 2.2863 0 0 0 1.0000 0 0 0 0 0 4.0000
2， 建立一个5×5矩阵，求它的行列式的值、迹、秩和范数。
>> A=[-29,6,18;20,5,12;-8,8,5]; >> b=eig(A) b =
-25.3169 -10.5182
16.8351
>> [v,b]=eig(A) v =
0.7130 0.2803 0.2733 -0.6084 -0.7867 0.8725 0.3487 0.5501 0.4050 b =
-25.3169 0 0 0 -10.5182 0 0 0 16.8351
3，下面是一个线性方程组：
????
?
?????=????????????????????52.067.095.06/15
/14
/15/14/13/14/13/12/1321x x x （1）求方程的解；
（2）将方程右边向量元素b 3改为0.53，再求解，并比较b 3的变化和解的相对变化。
（1）
>> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; >> C=[0.95;0.67;0.52];
>> B=A\C B =
1.2000 0.6000 0.6000
>> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];C=[0.95;0.67;0.52]
C =
0.9500 0.6700 0.5200
（2）
>> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];
>> C=[0.95;0.67;0.53];
>> B=A\C
B =
3.0000
-6.6000
6.6000
3，建立A矩阵，试比较sqrtm(A)和sqrt(A)，分析它们的区别。
>> A=magic(3);
>> sqrtm(A)
ans =
2.7065 + 0.0601i 0.0185 + 0.5347i 1.1480 - 0.5948i
0.4703 + 0.0829i 2.0288 + 0.7378i 1.3739 - 0.8207i
0.6962 - 0.1430i 1.8257 - 1.2725i 1.3511 + 1.4155i
>> sqrt(A)
ans =
2.8284 1.0000 2.4495
1.7321
2.2361 2.6458
2.0000
3.0000 1.4142
区别：sqrt是求矩阵里每个元素的平方根，并组成一个新的矩阵，而sqrtm是求矩阵的方根并组成一个新的矩阵。
实验三 选择结构程序设计
1，求分段函数的值。
??
???≠≠<≤≠<--+--+=其它及且－且325x 03x 0,1,65,622
2x x x x x x x x x y
用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。 程序：
>> m=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0]; for x=m
if x<0&x~=-3 y1=x^2+x-6;
disp(['y= ',num2str(y1)])
elseif x>=0&x<5&x~=2&x~=3 y2=x^2-5*x+6; disp(['y= ',num2str(y2)]) else
y3=x^2-x-1;
disp(['y= ',num2str(y3)]) end end y= 14 y= 11 y= 2 y= 1 y= -0.25 y= 5 y= 19
2，输入一个百分制成绩，要求输出成绩等级A 、B 、C 、D 、E 。其中90分～100分为A ，80分～89分为B ，70～79分为C ，60～69分为D ，60分以下为E 。
要求：
（1） 分别用if 和swtich 语句实现；
（2） 输入百分制成绩后要判断成绩的合理性，对不合理的成绩应输
出出错信息。
（1）if 语句 程序：
>> x=input('请输入一个分数:');
请输入一个分数:90
>> if x>=90
y='A';
elseif x>=80
y='B';
elseif x>=70
y='C';
elseif x>=60
y='D';
else y='E';
end
>> y
y =
A
程序：
>> x=input('请输入一个分数:');
请输入一个分数:120
>> if x>100|x<0
x=input('输入错误，请重新输入成绩:');
end
输入错误，请重新输入成绩:80
>> if x>=90
y='A';
elseif x>=80
y='B';
elseif x>=70
y='C';
elseif x>=60
y='D';
else y='E';
end
>> y
y =
B
（2）
>> x=input('请输入成绩：');
请输入成绩：102
>> while (x>100|x<0)
x=input('输入错误，输入错误请重新输入成绩');
end
输入错误，输入错误请重新输入成绩68
>> switch fix(x/10)
case {9,10}
y='A';
case {8}
y='B';
case {7}
y='C';
case {6}
y='D';
case {0,1,2,3,4,5}
y='E';
end
>> y
y =
D
4，建立5×6矩阵，要求输出矩阵第n行元素，当n值超过矩阵的行数时，自动转为输出矩阵最后一行元素，并给出出错信息。
程序：
>> a=rand(5,6)
a =
0.0153 0.4186 0.8381 0.5028 0.1934 0.6979
0.7468 0.8462 0.0196 0.7095 0.6822 0.3784
0.4451 0.5252 0.6813 0.4289 0.3028 0.8600
0.9318 0.2026 0.3795 0.3046 0.5417 0.8537
0.4660 0.6721 0.8318 0.1897 0.1509 0.5936
>> n=input('请输入要输出的第n：');
if n>5|n<0
b=a(5:5,:);
disp(['输入错误，程序输出最后一行的元素：',num2str(b)])
else
c=a(n:n,:);
disp(['输出结果是：',num2str(c)])
end
请输入要输出的第n：6
输入错误，程序输出最后一行的元素：0.46599 0.67214 0.8318 0.18965
电子科技大学电子工程学院标准实验报告（实验）课程名称MATLAB与数值分析 学生姓名：李培睿 学号：2013020904026 指导教师：程建
一、实验名称 《MATLAB与数值分析》第一次上机实验 二、实验目的 1. 熟练掌握矩阵的生成、加、减、乘、除、转置、行列式、逆、范数等运算 操作。（用.m文件和Matlab函数编写一个对给定矩阵进行运算操作的程序） 2. 熟练掌握算术符号操作和基本运算操作，包括矩阵合并、向量合并、符号 转换、展开符号表达式、符号因式分解、符号表达式的化简、代数方程的符号解析解、特征多项式、函数的反函数、函数计算器、微积分、常微分方程的符号解、符号函数的画图等。（用.m文件编写进行符号因式分解和函数求反的程序） 3. 掌握Matlab函数的编写规范。 4、掌握Matlab常用的绘图处理操作，包括：基本平面图、图形注释命令、 三维曲线和面的填充、三维等高线等。（用.m文件编写在一个图形窗口上绘制正弦和余弦函数的图形，并给出充分的图形注释） 5. 熟练操作MATLAB软件平台，能利用M文件完成MATLAB的程序设计。 三、实验内容 1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。并以x， y为坐标显示图像 x(n+1) = a*x(n)-b*(y(n)-x(n)^2); y(n+1) = b*x(n)+a*(y(n)-x(n)^2) 2. 编程实现奥运5环图，允许用户输入环的直径。 3. 实现对输入任意长度向量元素的冒泡排序的升序排列。不允许使用sort 函数。 四、实验数据及结果分析 题目一： ①在Editor窗口编写函数代码如下：
1.以下两种说法对吗？ （1）MATLAB进行数值的表达精度与其指令窗中的数据显示精度相同。 （2）MATLAB指令窗中显示的数据有效位数不超过七位。 2.历史指令窗所记录的内容与diary指令所产生的“日志”内容有什么不同？ DIARY filename causes a copy of all subsequent command window input and most of the resulting command window output to be appended to the named file. If no file is specified, the file 'diary' is used. DIARY OFF suspends it. DIARY ON turns it back on. DIARY, by itself, toggles the diary state. Use the functional form of DIARY, such as DIARY('file'), when the file name is stored in a string. 3.如何把用户自己的“工作目录”永久地设置在MATLAB的搜索路径上？“位于搜索路 径上的目录”与“当前目录”在MATLAB中的功用相同吗？ 4.如何向MATLAB工作空间输入一个含有100个左右元素的一维或二维数值数组？用直 接键入法？用数组编辑器？用M文件编辑器？ 5.运用数组算术运算符去掉下面程序里的for/end循环： x=11:15 for k=1:length(x) z(k)=x(k)^2+2.3*x(k)^0.5; end x=11:15 z1=x.^2+2.3*x.^0.5 6.不使用数组算术运算符，重写下面的程序代码： x=[2 1 4] z=1./(1+x.^2) x=2;k=1; while i<=4, z2(k)=1/(1+i^2); i=i+1;x=x+1; end 7．某公司销售电脑打印机的价格方案如下： （）如果顾客只买一台打印机，则一台的基本价格为$150。 （）如果顾客购买两台以上打印机，则第二台价格为$120。 （）第三台以后，每台$110。
实验指导 1, >> figure('Color',[1,0,0],'WindowButtonDownFcn','text(0.5,0.5,''Left Button Pressed'')');axis off 2, (1)默认属性 >> x=linspace(0,600,100000);y=(log(x+sqrt(1+x.*x)))/2;line(x,y) 句柄操作 >> x=linspace(0,600,100000); y=(log(x+sqrt(1+x.*x)))/2; line(x,y,'linewidth',3,'linestyle','-.','color','r'); text(300,3.5,'文字标注') (2)默认属性 >> t=linspace(0,20,500);x=t.*t;y=5.*t.*t;line(x,y) 句柄操作 >> t=linspace(0,20,500); x=t.*t;y=5.*t.*t;line(x,y); line(x,y,'linewidth',3,'linestyle','-.','color','r'); text(250,1600,'文字标注') 3, (1) >> x=linspace(0,50,1000); [x,y]=meshgrid(x); z=x.*x+y.*y-5.*sin(x.*y); axes('view',[-37.5,30]); surf(x,y,z); light('position',[10,20,4000]); shading interp (2) >> x=linspace(0,50,1000); [x,y]=meshgrid(x); z=y.*y.*y; axes('view',[-37.5,30]); surf(x,y,z); light('position',[30,20,1300000]); shading interp 4, >> x=-2*pi:0.01:2*pi; y1=sin(x); y2=cos(x); axes('position',[0.1,0.6,0.2,0.2]); plot(x,y1); axes('position',[0.6,0.6,0.2,0.2]); plot(x,y2); axes('position',[0.1,0.1,0.2,0.2]); fplot('tan(x)',[-1.5,1.5]); axes('position',[0.6,0.1,0.2,0.2]);
MATLAB程序设计语言 实 验 报 告 专业及班级：电子信息工程 姓名：王伟 学号：1107050322 日期 2013年6月20日
实验一 MATLAB 的基本使用 【一】 实验目的 1.了解MATALB 程序设计语言的基本特点，熟悉MATLAB 软件的运行环境； 2.掌握变量、函数等有关概念，掌握M 文件的创建、保存、打开的方法，初步具备将一般数学问题转化为对应计算机模型处理的能力； 3.掌握二维图形绘制的方法，并能用这些方法实现计算结果的可视化。 【二】 MATLAB 的基础知识 通过本课程的学习，应基本掌握以下的基础知识： 一. MATLAB 简介 二. MATLAB 的启动和退出 三. MATLAB 使用界面简介 四. 帮助信息的获取 五. MATLAB 的数值计算功能 六. 程序流程控制 七. M 文件 八. 函数文件 九. MATLAB 的可视化 【三】上机练习 1. 仔细预习第二部分内容，关于MATLAB 的基础知识。 2. 熟悉MATLAB 环境，将第二部分所有的例子在计算机上练习一遍 3. 已知矩阵???? ??????=??????????=123456789,987654321B A 。求A*B ，A .* B ，比较二者结果是否相同。并利用MATLAB 的内部函数求矩阵A 的大小、元素和、长度以 及最大值。 程序代码： >> A=[1 2 3;4 5 6;7 8 9]; >> B=[9 8 7;6 5 4;3 2 1]; >> A*B ans =
30 24 18 84 69 54 138 114 90 >> A.*B ans = 9 16 21 24 25 24 21 16 9 两者结果不同 >> [m,n]=size(A) m = 3 n = 3 >> b=sum(A) b = 12 15 18 >> a=length(A) a = 3 >>max(A) ans =
<<<<<<精品资料》》》》》 Matlab 上机测试题答案 班级 姓名 学号 要求: （1）请将程序和运行结果粘贴到题目下面，写清楚哪部分是程序，哪部分是运行结果，运行结果包括图形，图形请适当缩小，然后将word 文档传到教育在线，注意保存格式是2003的doc ，文件名”自111_110110_张三” （2）在最开始前请将输入Matlab 命令窗口中输入”diary ”命令，随后在workspace 窗口中可以发现一个名为“diary ”的文件。在老师要求结束练习后，在命令窗口输入”diary off ”命令，然后将文件“diary ”上传到教育在线。 （3）注意题号不要弄乱。 --------------------------------------------------------------------- 1、 矩阵Y= ???? ? ???? ???3472123100451150425 ，给出元素1的全下标和单下标，并用函数练习全下标和单下标的转换，求出元素100的存储位置。取出子矩阵? ? ? ? ??21301，并求该矩阵的维数。 解：命令为： Y=[5,2,4;0,15,1;45,100,23;21,47,3] Y(2,3) Y(10) sub2ind([4 3],2,3) [i,j]=ind2sub([4 3],10) find(Y==100) sub2ind([4 3],3,2) B=Y(2:2:4,3:-2:1) 或 B=Y([2 4],[3 1]) [m n]=size(Y) --------------------------------------------------------------------- 2、已知矩阵A=[1 0 -1 ；2 4 1； -2 0 5]，B=[0 -1 0；2 1 3；1 1 2] 求2A+B 、A 2-3B 、A*B 、B*A 、A .*B ，A/B 、A\B 解：命令为： A=[1 0 -1 ；2 4 1； -2 0 5] B=[0 -1 0；2 1 3；1 1 2] E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B --------------------------------------------------------------------- 3、利用函数产生3*4阶单位矩阵和全部元素都为8的4*4阶矩阵，并计算两者的乘积。 解：命令为： A=eye(3,4) B=8*ones(4) C=A*B --------------------------------------------------------------------- 4、创建矩阵a=????? ???????------7023021.5003.120498601 ,取出其前两列构成的矩阵b ，取出前两行构成矩阵c ，转置矩阵b 构成矩阵d ，计算a*b 、c 练习题 1.求函数在指定点的数值导数 x=sym('x'); >> y=[x x.^2 x.^3;1 2*x 3*x.^2;0 2 6*x]; >> x=1; >> eval(diff(y)) ans = 1 2 3 0 2 6 0 0 6 >> x=2; >> eval(diff(y)) ans = 1 4 12 0 2 12 0 0 6 >> x=3; >> eval(diff(y)) ans = 1 6 27 0 2 18 0 0 6 2.求下列函数导数 （1） x=sym('x'); >> y=x^10+10^x+(log(10))/log(x); >> diff(y) ans = 10*x^9+10^x*log(10)-2592480341699211/1125899906842624/log(x)^2/x （2） x=sym('x'); >> y=log(1+x); >> x=1; >> eval(diff(y,2)) %在x=1的条件下对y表达式求两次导数后导函数的值 ans = -0.2500 3.用数值方法求下列积分 首先先讲一下trapz的用法，如下题 t=0:0.001:1; >> y=t; >> trapz(t,y) ans = 0.5000 （1） >> x=1:0.01:5; >> y=(x.^2).*sqrt(2*x.^2+3); >> trapz(x,y) ans = 232.8066 (2) x=pi/4:0.01:pi/3; >> y=x./(sin(x).^2); >> trapz(x,y) ans = 0.3810 第三题拟合曲线题 x=[0:0.1:1]; >> y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; >> a=polyfit(x,y,2); >> x=[0.05:0.2:1.05]; >> y=a(3)+a(2)*x+a(1)*x.^2 %注意x要在y前先赋值，不然y不会运行为最新的x对呀的y值 y = MATLAB 数学实验报告 指导老师： 班级： 小组成员： 时间：201_/_/_ Matlab 第二次实验报告 小组成员： 1 题目：实验四，MATLAB 选择结构与应用实验 目的：掌握if 选择结构与程序流程控制，重点掌握break，return , pause语句的应用。 问题：问题1：验证“哥德巴赫猜想” ，即：任何一个正偶数（n>=6）均可表示为两个质数的和。要求编制一个函数程序，输入一个正偶数，返回两个质数的和。 问题分析：由用户输入一个大于6 的偶数，由input 语句实现。由if 判断语句判断是否输入的数据符合条件。再引用质数判断函数来找出两个质数，再向屏幕输出两个质数即可。 编程：function [z1,z2]=gede(n); n=input('please input n') if n<6 disp('data error'); return end if mod(n,2)==0 for i=2:n/2 k=0; for j=2:sqrt(i) if mod(i,j)==0 k=k+1; end end for j=2:sqrt(n-i) if mod(n-i,j)==0 k=k+1; end end if k==0 fprintf('two numbers are') fprintf('%.0f,%.0f',i,n- i) break end end end 结果分析 如上图，用户输入了大于6的偶数返回两个质数5和31,通过 不断试验，即可验证哥德巴赫猜想。 纪录：if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰，更快的解决问题。 2题目：实验四，MATLAB选择结构与应用实验 目的：用matlab联系生活实际，解决一些生活中常见的实际问 题。 MATLAB上机习题一 请按以下步骤完成上机实验： 1）在FTP上下载“MATLAB上机习题一.doc”文件，所有习题列在该文件内； 2）在MATLAB中完成所有习题，并将屏幕截图粘贴到相应习题后面； 3）如果习题是问答题，请将答案写在题目后； 4）如果有的习题要求提供脚本文件，请将脚本文件内容拷贝到相应习题后； 5）将文件保存并重命名为“自己的学号-姓名”，例如“20110771-张三.doc”； 6）上传该文件到FTP的相关目录。 1. 运行MATLAB软件，拖放、关闭界面上的子窗口，并恢复到原始试图。 2. 采用鼠标及命令两种方式将桌面添加到MATLAB搜索路径列表的起始及最后位置。 3. 采用鼠标及命令方式将当前工作目录设置为桌面指向的文件夹。 4. 通过使用帮助确定内置变量ispc 的含义。 6. 观察MATLAB中关键字、字符串、注释的字体显示颜色。 5. 创建double类型的变量，并进行计算 1）a=87，b=190，计算a+b、a-b、a*b； 2）创建uint8类型的变量，数值与（1）中相同，进行相同的计算，观察计算结果与预想的是否一致。 6. 计算如下表达式： 1）()sin 60 2）3e 3）3cos 4π?? ??? 4）2 7562323336 +?-?? 7. 设u=2，v=3，计算： 1）4log uv v 2）() 2 2u e v v u +- 3 8. 计算如下表达式： 1）()()3542i i -+ 2）()sin 28i - 9. 判断下面语句的运算结果，并与MATLAB 计算结果做比较： 1）420< 2）420≤ 3）420== 4）4~20= 5）''''b B < 10. 设a=39，b=58，c=3，d=7，判断下面表达式的值与MATLAB 计算结果做比较： 1）a b > 2）a b < 3）&&a b b c >> 4）a d == 5）|a b c > 6）~~d 11. 判断下列变量名哪些是合法的MATLAB 变量： 1）fred 7）fred! 2）book_1 8）book-1 3）2ndplace 9）Second_Place 4）#1 10）No_1 5）vel_5 11）vel.5 6）tan 12）while 提示，可以使用isvarname 对上述变量名进行检验。 12. 编写脚本程序，命名为barbell.m ，完成如下计算： 如图1所示的杠铃， 1 求函数在指定点的导数值 () 23 2 123,1,2,3 026 x x x f x x x x x == >> syms x >> a=[x x^2 x^3;1 2*x 3*x^2;0 2 6*x]; >> f=det(a); >> diff(f,1) ans = 6*x^2 >> diff(f,2) ans = 12*x >> diff(f,3) ans = 12 2 符号法求下列函数的导数或积分 1）y=x10+10x+log x 10,求y’ f=('x^10+10^x+log(10)/log(x)') f = x^10+10^x+log(10)/log(x) >> diff(f) ans = 10*x^9+10^x*log(10)-log(10)/log(x)^2/x 2）y=ln(1+x), 求y’’∣ x=1 f=('log(1+x)/log(e)') f = log(1+x)/log(e) >> diff(f,1,2) ans = -1/(1+x)^2/log(e) 3) y=e x/cosx,求y’ f=('exp(x)/cos(x)') f = exp(x)/cos(x) >> diff(f) ans = exp(x)/cos(x)+exp(x)/cos(x)^2*sin(x) 4） function f=fun0(t) f=t*sin(t) int('fun0','0','pi') ans = 1/2*pi^2 5) 已知函数z=sin(xy), 计算 syms x y >> z=('sin(x*y)') z = sin(x*y) >> diff(diff(z,y,2),x) ans = -cos(x*y)*y*x^2-2*sin(x*y)*x 3 用数值方法求定积分 1） function f=fun(x) f=x.^2.*sqrt(2.*x.^2+3) quad('fun',1,5) ans = 232.8057 2） function f=fun(x) f=x./sin(x).^2 quad('fun',pi/4,pi/3) ans = 0.3835 4 已知数据[x,y]如下表，试求2次拟合多项式f(x)，然后求x=0.05,0.25,0.45,0.65,0.85,1.05各点的函数近似值，并绘出拟合曲线及求得 信号与系统MATLAB第一次实验报告 一、实验目的 1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。 2.学会运用MATLAB表示常用连续时间信号的方法 3.观察并熟悉一些信号的波形和特性。 4.学会运用MATLAB进行连续信号时移、反折和尺度变换。 5.学会运用MATLAB进行连续时间微分、积分运算。 6.学会运用MATLAB进行连续信号相加、相乘运算。 7.学会运用MATLAB进行连续信号的奇偶分解。 二、实验任务 将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。 三、实验内容 1.MATLAB软件基本运算入门。 1). MATLAB软件的数值计算： 算数运算 向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn 为结束值。 矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开； 矩阵的不同行之间必须用分号”;”或者ENTER分开。2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。 举例：计算一个函数并绘制出在对应区间上对应的值。 2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名” 2.MATLAB软件简单二维图形绘制 1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y) 2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表 示第p个区域，表达为subplot(mnp)或者subplot(m,n,p) 3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin]) 4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’) 5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’) 6).输出：grid on 举例1： 举例2： (1)在MATLAB的命令窗口中执行_____命令，将命令窗口的显示内容清空。（） A.clear B.clc C.echo off D.cd (2)在MATLAB的命令窗口中执行_____命令，使数据输出显示为十六进制表示。（） A.format long B.format rat C.format hex D.format short e (3)下列变量名中_____是合法的。（） A.x*y,a,1 B.x\y,a1234 C.end,1 bcx D.char_1,i,j (4)已知x=0:5,则x有_____个元素。（） A.5 B.6 C.7 D.8 (5)一下运算符中哪个的优先级最高_____。（） A./ B.^ C.~= D.& (6)使用检测函数isnumeric(10)的结果是_____。（） A.1 B.0 C.false D.true (7)三维图形中默认视角是_____度。（） A.方位角=0 俯仰角=90 B.方位角=90 俯仰角=0 C.方位角=37.5 仰俯角=30 D.方位角=0 仰俯角=180 (8)将符号表达式化简为因式分解因式分解因式分解因式分解形式，使用_____函数。（） A.collect B.expand C.horner D.factor (9)运行以下命令，则_____描述是正确的。（）>>syms a b c d >>A=[a b;c d] A.A占用的内存小于100B B.创建了5个符号变量 C.A占用的内存是a b c d的总和 D.不存在 （10）已知数组a=[1 2 3;4 5 6;7 8 9],则a(:,end)是指_____元素。 （11）运行命令bitor（8,7）的结果是_____。 （12）运行以下命令： >>x=0:10; >>y1=sin(x); >>y2=5*sin(x); >>y3=[10*sin(x );20*sin(x)]; >>plot(x,y1,x,y2,x,y3) 则在一个图形窗口中，可以看到_____条曲线。 （13）符号表达式“g=sym(sin(a*z)+cos(w*v))”中的自由符号变量是_____。 （14）运行以下命令： >>syms t >>f1=1/t >>limitf1_r=limit(f1,'t','0','right'); 则函数limitf1_r趋向0的右极限为_____。 15.在MATLAB的命令窗口中执行______命令，使数值5.3显示为5.300000000000000e+000 A. format long B. format long e C. format short D. format short e 16.下列变量名中______是合法的。A．char_1,i,j B．1_1, a.1 C．x\y,a1234 D.end,1bcx 17.已知x=0:9,则x有_____个元素。 A.12 B.11 C.10 D.9 18.产生对角线上为全1其余为0的2行3列矩阵的命令是______ A. ones(2,3) B. ones(3,2) C. eye(2,3) D. eye(3,2) 19.已知数组a= [1 2 3 4 5 6 7 8 9] ，则运行a(:,1)=[]命令后______ A. a变成行向量 B. a数组为2行2列 C. a 数组为3行2列 D. a数组中没有元素3 20.按含义选出各个函数名：表示4舍5入到整数的是____，表示向最接近0取整的是____,表示向最接近-∞取整的是____,表示向最接近∞取整的是_____ A. round(x) B. fix(x) C. floor(x) D. ceil(x) 21.已知a=0:5,b=1:6,下面的运算表达式出错的为______ A. a+b B. a./b C. a’*b D. a*b 22.已知s=’显示”hello”’,则s的元素个数是______ A. 12 B. 9 C. 7 D.18 第二章 1·求下列表达式的值。 (1）w=sqrｔ(２)*（1+0.34245*10＾-6） ｗ= 1.4１42 （2)a＝3.5;b=５；c=-9．8; x=（2＊pi＊a＋(ｃ+b）/(pi＋a＊ｂ*c)-exp（２))/(ｔａn（ｂ+c)+a)；x x ＝ ０.98２9 （3)a=３．32;b=-7.9； y＝２*pi*a＾２*[(1-ｐi/４)*b-(0.８33３-ｐi／4)*a］; y ｙ= -128.4２7１ (4）t=[2,１-3i；５，－０.65］； z=１／2*exp（2＊t)*ｌog（t+ｓｑrｔ(１＋t^２)); z z= 1.０e+00４＊ 0.0057 － 0.000７i 0．０04９ - 0.0０2７i 1.９8８4 －０.3６9６i 1.7706 - １.0５39i 2,已知a,ｂ,求下列表达式的值。 a=[-1,5,－４；０,7，８;3,6１,7];b=[８,3,-1;２,5，3;-3,2，0]; (1）a+6*b ans= 47２3－10 12 3726 -15 73 7 a^２-b+eye（3) ans= －１8 -２１7 17 22 53３ 109 ２1867 52６ (２)a＊b anｓ＝ 14 14 1６ -10 5１ 21 125 328 １80 a.*b anｓ = -8 １5 ４ 0 35 24 -9 1２2 0 b*a ans ＝ -１１０-15 7 22８５3 3 -1 28 (3）a/b anｓ = 1.２234 -0．925５ 2．9787 －0.９468 2.3511 -0．９574 ４.6170 3.8723 13．8９３6 ｂ\a ａns = -0．51０6 -8.617０ -1.1277 ０.７340 17.5７45 1.8085 －0.8８30 -21.21２8 0.4043 （４)[a,b] anｓ= －１ 5 -4 ８ 3 －１ 0 7 8 2 5 3 ３６1 7 -3 2 0 ［a([1,３],:);b^2] aｎs ＝ -１５ -4 3 61 7 73 3７ 1 17 37 １3 -20 1 ９ 3.已知ａ,完成下列操作。 a=[２3,１0,－0.778,0;4１,-45，65,５;32,5,０,32;6,-９.５4,5４,3.1４］; （1）输出a在[10，25]范围内的全部元素。 k=ｆｉｎd(ａ>1０＆a<2５） ａ(k） k = 1 ans＝ 23 (２)取出a前３行构成矩阵b，前两列构成矩阵c,右下角3＊2子矩阵构成矩阵ｄ,ｂ与c的乘积构成矩阵e。 ｂ＝a(1:3,:) b = ２3.0000 １０.0000 -0.77８0 0 41.0000 -45.0000 65．0０0０５.00００ ３2．０000 5.0000 0 32.0000 c=a（:,1:2) c = 23.00００ 10.０000 41.0０00 -4５.0000 32.00０0 ５.00００ ６.0000 -9.54０0 d=a(2:4,3：4） d = 65．000０ 5．0０0０ 0 32.00００ ５４.00０0 3．１４00 e=b*c e = 1.0e+003 * ０.9141 -0.2２39 1 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。 for m=100:999 m1=fix(m/100); %求m的百位数字 m2=rem(fix(m/10),10); %求m的十位数字 m3=rem(m,10); %求m的个位数字 if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end 2.从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。 sum=0; n=0; val=input('Enter a number (end in 0):'); while (val~=0) sum=sum+val; n=n+1; val=input('Enter a number (end in 0):'); end if (n > 0) sum mean=sum/n end 3. 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求[1,500]之间的全部完数。 for m=1:500 s=0; for k=1:m/2 if rem(m,k)==0 s=s+k; end end if m==s disp(m); end end 4. 从键盘上输入数字星期，在屏幕上显示对应英文星期的单词。 function week n=input('input the number:'); if isempty(n) errror('please input !!') end if n>7|n<1 error('n between 1 and 7') end switch n case 1 disp('Monday') case 2 disp('Tuesday') case 3 disp('Wednesday') case 4 disp('Thursday') case 5 disp('Friday') case 6 disp('Saturday') case 7 disp('Sunday') end 5. 某公司销售电脑打印机的价格方案如下： （）如果顾客只买一台打印机，则一台的基本价格为$150。 （）如果顾客购买两台以上打印机，则第二台价格为$120。 （）第三台以后，每台$110。 写一段程序分别计算出购买1--10台打印机所需的钱数。打印机台数可以在程序开始处指定，或通过input命令读入。运行程序，计算出购买10台打印机的总价格。 写出程序，生成分别购买1--10台打印机所需价格的图表（使用fprintf命令输出图表，不允许手算）。 x=input('请输入购买的打印机台数：'); for m=1:x if m<=1 y(m)=150*m; elseif m<=2 y(m)=150+120*(m-1); else y(m)=150+120+110*(m-2); y(1,m)=y(m); end end y(x) plot(1:m,y,'r*--') Matlab 与科学计算考试样题（客观题） 1 下面的MATLAB 语句中不正确的有： a) 2a ＝pi; b) record_1=3+4i c) a=2.0, d) c=1+6j 2. 已知水的黏度随温度的变化公式如下，其中a=0.03368，b=0.000221，计算温度t 为20，30，40度时的粘度分别是： 02 1at bt μμ=++0μ为0℃水的黏度，值为3 1.78510-?；a 、b 为常数，分别为0.03368、0.000221。 3. 请补充语句以画出如图所示的图形： [x,y]=meshgrid(-2:0.1:2, -2:0.1:2); Z=x.*exp(-x.^2-y.^2); ; a) Plot3(x,y,Z) b) plot3(x,y,Z) c) mesh(x,y,Z) d) plot3(x,y,z) 2 a) 0.4900 1.2501 0.8560 b) 0.8560 1.2501 0.4900 c) -0.6341 3.8189 -3.7749 d) 3.8189 -3.7749 2.8533 解释说明： >> x=0.5:0.5:3.0; >> y=[1.75,2.45,3.81,4.80,8.00,8.60]; >> a=polyfit(x,y,2) a = 0.4900 1.2501 0.8560 >> x1=[0.5:0.25:3.0]; >> y1=a(1)*x1.^2+a(2)*x1+a(3) >> plot(x,y,'*') >> hold on >> plot(x1,y1,'--r') 5. 求方程在x=0.5附近的根. 21x x += a) 0.6180 b) -1.1719e-25 c) －1 d) -1.6180 6. 用Newton-Cotes 方法计算如下积分 1 5x ? （a ）133.6625 (b) 23.8600 (c) 87.9027 (d) -1.6180 7. y=ln(1+x)，求x=1时y" a) -0.25 b) 0.5 c) -0.6137 d) -1.6137 8. 某公司用3台轧机来生产规格相同的铝合金薄板。取样测量薄板的 厚度，精确至‰厘米。得结果如下： 轧机1：0.236 0.238 0.248 0.245 0.243 轧机2：0.257 0.253 0.255 0.254 0.261 轧机3：0.258 0.264 0.259 0.267 0.262 计算方差分析结果，并判定各台轧机所生产的薄板的厚度有无显著的差异？ a) p ＝1.3431e-005，没有显著差异。 b) p ＝0.9688，没有显著差异。 c) p ＝0.4956，有显著差异。 d) p ＝0.9688，有显著差异。 22x y x y e x y e --?-=??-+=?? MATLAB 程序设计实验报告 一、实验目的 1. 通过实验熟悉MATLAB 仿真软件的使用方法； 2. 掌握用MATLAB 对连续信号时域分析、频域分析和s 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征； 3. 掌握用MATLAB 对离散信号时域分析、频域分析和z 域分析的方法，利用绘图命令绘制出典型信号的波形，了解这些信号的基本特征； 4. 通过绘制信号运算结果的波形，了解这些信号运算对信号所起的作用。 二、实验设备 1. 计算机 ： 2. MATLAB R2007a 仿真软件 三、实验原理 对系统的时域分析 信号的时域运算包括信号的相加、相乘，信号的时域变换包括信号的平移、反折、倒相及信号的尺度变换。 （1）信号的相加和相乘：已知信号)(1t f 和)(2t f ，信号相加和相乘记为 )()(1t f t f =)(2t f +；)()(1 t f t f =)(2t f *。 （2）信号的微分和积分：对于连续时间信号，其微分运算是用diff 函数来完成的，其语句格式为：diff(function,’variable’,n)，其中function 表示需要进行求导运算的信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；n 为求导的阶数，默认值为求一阶导数。连续信号的积分运算用int 函数来完成，语句格式为：diff(function,’variable’,a,b)，其中function 表示需要进行被积信号，或者被赋值的符号表达式；variable 为求导运算的独立变量；a,b 为积分上、下限，a 和b 省略时为求不定积分。 （3）信号的平移、翻转和尺度变换 信号的平移包含信号的左移与右移，信号的翻转包含信号的倒相与折叠，平移和翻转信号不会改变信号)(t f 的面积和能量。信号的尺度变换是对信号)(t f 在时间轴上的变化，可使信号压缩或扩展。)(at f 将原波形压缩a 倍，)/(a t f 将原波形扩大a 倍。 ￥ 对系统频率特性的分析 1设置matlab 的工作环境，将工作目录设置为d:\work ，添加搜索目录d:\example 2在matlab 的命令窗口里完成如下计算，其中t 的值分别取-1,0,1，表达式如下： 4/3)2 t y e π-= 3自行产生一个5行5列的数组，得到最中间的三行三行矩阵。 4用magic 产生一个5*5的矩阵，将这个矩阵的第二行与第三行互换位置 5求方程组的根 x 1+4x 2-3x 3=2 2x 1+5x 2-x 3=11 x 1+6x 2+x 3=12 6已知：一个多项式的系数向量是p=[1 -6-72 -27]，求这个多项式的根。 7已经两个多项式的系数分别是：[1 2 3 4]和[1 4 9 16]，请求这两个多项式的乘积，及商和余数。 8给定一个多项式的根是[-5 -3+4i -3-4i],求原来的多项式 9 A=[2 3 4； 1 5 7； 6 2 5] 用什么函数，保证第一列排序的时候，其他列跟着变化。最后的结果是： [1 5 7； 2 3 4； 6 2 5] 10补充题：电路分析 电路分析常常涉及到对方程组的求解，这些方程常常是利用描述电流进入和离开节点的电流方程，或者描述电路中网络回路上的电压的电压方程得出的。下图描述了3个回路电压的方程 方程式如下： 11121222132423432 532()0()()0()0V R I R I I R I I R I R I I R I I R I V -++-=??-++-=??-++=? 假设5个电阻值为已知，2个电压值也为已知，求3个电流值。 11自行产生一个5行5列的数组，用两种方法得到最中间的三行三行矩阵。 12根据a=reshape(-4:4,3,3)做一个矩阵，然后（1）取出所有大于0的元素构成一个向量（可推广到从一个矩阵里筛选出符合条件的元素组成一个向量）（2）将原矩阵中大于0的元素正常显示，而小于等于0的元素全部用0来表示(可推广到将原矩阵中不符合条件的全用0来表示，符合条件的值不变)。(思考：大于0的正常显示，小于等于0的用-10来表示) 13建立如下一个元胞数组，现在要求计算第一个元胞第4行第2列加上第二个元胞+第三个元胞里的第二个元素+最后一个元胞的第二个元素。 a={pascal(4),'hello';17.3500,7:2:100} 14建立一个结构体的数组，包括3个人，字段有姓名，年龄，分数，其中分数由随机函数产生一个3行10列的数据表示了有10门课程，每门课程有三个阶段的分数。问题是： 问题1，如何找到第2个人的分数并显示出来 问题2，如何找到第2个人的每门课程3个阶段的平均分数并显示出来 问题3，全班同学（指这3个学生）的10门课程的每门课程的平均分如何计算出来？要求放到一个数组里。 问题4，找到这个班所有同学的姓名放到一个元胞数组里 15给定一个图像文件，格式是jpg ，通过inportdata 引入这个文件，查看引入后数据保存是一个数组，是100*100*3的一个数组，这是一个三维的数组，表明有100*100个点，每个点有RGB 三个方面决定的。现在，要求将这些点的中间部分30个象素宽度的正方形用黑点表示。 16手动构造一个长宽页是1000*1000*3的图像，每个图像的点有三个0-255之间的随机值构成。然后使用image 函数显示这个图像，再用imwrite 函数保存这个图像，使这个图像用图片浏览器也能够查看。 17给定一个矩阵 a = 1 2 3 4 5 6 编写一个M 函数，要求输入是a ，输出有三个：平均数，标准差，秩。程序运行 实验报告 2. The Branching statements 一、实验目的： 1.To grasp the use of the branching statements; 2.To grasp the top-down program design technique. 二、实验内容及要求： 1．实验内容： 1).编写 MATLAB 语句计算 y(t)的值 (Write the MATLAB program required to calculate y(t) from the equation) ???<+≥+-=0 530 53)(2 2t t t t t y 已知 t 从-5到 5 每隔0.5取一次值。运用循环和选择语句进行计算。 (for values of t between -5 and 5 in steps of 0.5. Use loops and branches to perform this calculation.) 2).用向量算法解决练习 1, 比较这两个方案的耗时。 （tic ，toc 的命令可以帮助你完成的时间计算，请使用'help'函数）。 Rewrite the program 1 using vectorization and compare the consuming time of these two programs. (tic, toc commands can help you to finish the time calculation, please use the …help ? function). 2.实验要求： 在报告中要体现top-down design technique, 对于 3 要写出完整的设计过程。 三、设计思路： 1．用循环和选择语句进行计算： 1).定义自变量t ：t=-5:0.5:5; 2).用循环语句实现对自变量的遍历。 3).用选择语句实现对自变量的判断，选择。 4).将选择语句置入循环语句中，则实现在遍历中对数据的选择，从而实现程序的功能。 2. 用向量法实现： 1).定义自变量t ：t=-5:0.5:5; 2).用 b=t>=0 语句，将t>=0得数据选择出，再通过向量运算y(b)=-3*t(b).^2 + 5; 得出结果。 3).用取反运算，选择出剩下的数据，在进行向量运算，得出结果。 四、实验程序和结果 1.实验程序 实验程序：创建m 文件：y_t.m中南大学材料学院科学计算与MATLAB考试题库
matlab实验报告
MATLAB上机习题一
中南大学matlab题目
MATLAB实验报告(1-4)
MATLAB上机考试题(一)
中南大学matlab课后习题(1)
Matlab 上机题及答案
中南大学Matlab与科学计算样题 (加主观题答案)
MATLAB程序设计实验报告
MATLAB上机题汇总(基本题,抽3题)合肥学院
matlab实验报告
相关主题
文本预览