练习1 基础练习
一、矩阵及数组操作：
1．利用基本矩阵产生3×3和15×8的单位矩阵、全1矩阵、全0矩阵、均匀分布随机矩阵（[-1，1]之间）、正态分布矩阵（均值为1，方差为4）。
解：
d=[1,2,3;4,5,6;7,8,9]
d =
1 2 3
4 5 6
7 8 9
eye(size(d))
ans =
1 0 0
0 1 0
0 0 1
eye(15,8)
ans =
1 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0
0 0 1 0 0 0 0 0
0 0 0 1 0 0 0 0
0 0 0 0 1 0 0 0
0 0 0 0 0 1 0 0
0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 1
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
ones(3,3)
ans =
1 1 1
1 1 1
1 1 1
zeros(3,3)
ans =
0 0 0
0 0 0
0 0 0
x =
0.9003 -0.0280 -0.0871
-0.5377 0.7826 -0.9630
0.2137 0.5242 0.6428
x=-1+2*rand(15,8)
x =
Columns 1 through 6
-0.1106 -0.7222 0.3443 -0.6983 0.0682 -0.6541
0.2309 -0.5945 0.6762 0.3958 0.4542 0.9595
0.5839 -0.6026 -0.9607 -0.2433 -0.3814 -0.4571
0.8436 0.2076 0.3626 0.7200 0.6770 -0.4953
0.4764 -0.4556 -0.2410 0.7073 0.1361 0.7515
-0.6475 -0.6024 0.6636 0.1871 -0.2592 0.4746
-0.1886 -0.9695 0.0056 -0.0069 0.4055 -0.7270
0.8709 0.4936 0.4189 0.7995 0.0931 -0.9765
0.8338 -0.1098 -0.1422 0.6433 -0.1102 0.7878
-0.1795 0.8636 -0.3908 0.2898 0.3891 -0.6017
0.7873 -0.0680 -0.6207 0.6359 0.2426 -0.4026
-0.8842 -0.1627 -0.6131 0.3205 0.5896 0.3229
-0.2943 0.6924 0.3644 -0.3161 0.9137 -0.4312
0.6263 0.0503 -0.3945 -0.4205 0.0452 -0.0616
-0.9803 -0.5947 0.0833 -0.3176 0.7603 -0.8704
Columns 7 through 8
0.9767 -0.0778
0.1656 0.1357
-0.1530 0.5884
0.0310 -0.8816
-0.3321 0.2057
-0.1342 -0.8995
-0.5481 -0.1693
0.1596 -0.3900
0.5207 0.7487
0.0596 -0.9700
0.2811 0.5359
-0.5819 0.9417
-0.2404 0.9802
0.5667 0.5777
0.3617 -0.1227
y=1+2*randn(3)
y =
1.3493 -0.1766 1.2279
0.6266 5.3664 3.1335
2.4516 0.7272 1.1186
y =
Columns 1 through 6
0.1349 1.2279 0.2002 2.0575 1.0001 0.5320
-2.3312 3.1335 2.3800 1.4386 0.3643 1.2369
1.2507 1.1186
2.6312 -0.8438
3.1900 1.6296
1.5754 0.8087
2.4238 -
3.3413 -2.7480 3.8870
-1.2929 -0.6647 3.5805 0.8816 1.8564 0.2981
3.3818 1.5888 2.3372 -1.0213 2.7913 2.2465
3.3783 -1.6724 3.3817 2.2289 2.4619 2.5981
0.9247 2.4286 -1.4049 2.0155 2.1557 2.8818
1.6546 4.2471 0.9604 4.3849 1.0806 -0.9842
1.3493 -0.3836 0.6866
2.1826 2.3542 1.4241
0.6266 2.7160 -2.2082 -0.2872 2.1378 1.4758
2.4516
3.5080 1.5146 1.7607 0.4887 -1.0155
-0.1766 -2.1875 -1.1129 -1.0182 0.2451 -0.4841
5.3664 -1.8819 3.8303 0.9610 0.4082 3.1646
0.7272 2.1423 -0.6102 0.9036 -1.9503 0.7370
Columns 7 through 8
1.7798 1.4681
1.1760 1.0429
-0.2709 -1.0079
-0.1191 -0.8943
1.8873 0.2511
-0.8998 -1.3718
2.5624 -1.1118
2.1379
3.9450
-0.6434 1.1115
0.4688 -1.4346
-1.3756 0.9175
-3.4046 -1.2567
2.9727 -1.6986
-0.0373 0.4778
1.6547
2.9069
2．利用fix及rand函数生成[0，10]上的均匀分布的10×10的整数随机矩阵a，然后统计a中大于等于5的元素个数。
解：
a=fix(0+11*rand(10))
a =
9 2 1 3 0 4 0 8 10 7
4 3 6 4 7 8 4 10 3 9
10 2 2 3 7 9 4 2 0 2
1 5 6 4 3 1 3 7 8 5
9 4 5 5 2 4 9 5 9 0
8 5 5 7 7 5 8 10 3 9
4 6
5 3 8 7 10 7 2 6
6 0 10 1 10 6 6 2 10 3
10 3 3 4 5 5 0 4 1 4
7 6 4 5 9 4 6 1 10 9
k=find(a>=5);
m=a(k);
length(m)
ans =
55
3．在给定的矩阵中删除含有整行内容全为0的行，删除整列内容全为0的列。解：
a=[1,0,0,0;0,0,0,0;0,1,1,0]
a =
1 0 0 0
0 0 0 0
0 1 1 0
a(2,:)=[]
a =
1 0 0 0
0 1 1 0
a(:,4)=[]
a =
1 0 0
0 1 1
二、绘图：
4．在同一图形窗口画出下列两条曲线图像：
y1=2x+5； y2=x^2-3x+1，
并且用legend标注。
解：
x=0:0.01:1000;
y1=2.*x+5;
y2=x.^2-3.*x+1;
plot(x,y1,x,y2)
legend('y1','y2')
5．画出下列函数的曲面及等高线：
z=x^2+y^2+sin(xy)．
解：
[x,y]=meshgrid(0:0.25:4*pi);
z=x.^2+y.^2+sin(x.*y);
meshc(x,y,z);
三、程序设计：
6．编写程序计算（x在[-3，3]，间隔0.01）
解：
x=input('请输入x的值:');
if x>=-3&x<-1
y=(-x.^2-4*x-3)./2;
elseif x>=-1&x<1
y=-x.^2+1;
elseif x>=1&x<=3
y=(-x.^2+4*x-3)./2;
end
y
7．有一列分数序列：
求前15项的和。
解：
a=1;
b=2;
sum=0;
for k=1:15
c=b/a;
sum=sum+c;
t=b;
b=a+b;
a=t;
end
sum
sum =
24.5701
8．用至少三种方法编写函数实现求任意整数n的阶乘。解：
方法一：
function f=factor(n)
n=input('please input n:');
if n<=1
f=1;
else
f=factor(n-1)*n;
end
方法二：
function result = fa(n)
n=input('please input n:');
result = 1;
for i = 1:n
result = result * i;
end
方法三：
n=input('please input n:');
x=1:n;
prod(x)
9．将任意大于6的偶数m写成两个素数p1、p2的和（试着写出所有的m=p1+p2的可能形式）。
解：
function y=gassumption(n);
n=input('please input n:');
if mod(n,2);
error('This numeber is not an even.Please input again.')
elseif n<=6;
error('n must be larger than 6')
else
for m=1:n;
for k=m:n;
if (isprime(m))&(isprime(k))&(m+k==n);
disp([num2str(n),'=',num2str(m),'+',num2str(k)]);
break;
end;
end;
end;
end;
10．是否任意3的倍数m可以写成两个素数p1、p2、p3的和（试着写出所有的m=p1+p2+p3
的可能形式）？
解：
function y=gassumption(n);
n=input('please input n:');
if mod(n,3);
error('Please input again.')
elseif n<6;
error('n must be larger than 6')
else
for m=1:n;
for k=m:n;
for p=k:n
if (isprime(m))&(isprime(k))&(isprime(p))&(m+k+p==n);
disp([num2str(n),'=',num2str(m),'+',num2str(k),'+',num2str(p)]);
break;
end;
end;
end;
end;
end;
四、数据处理与拟合初步：
比。
解：
t=1:10;
y=[4.842,4.362,3.754,3.368, 3.169,3.038,3.034,3.016,3.012,3.005]; u=exp(-t);
p=polyfit(u,y,1);
tt=1:0.05:10;
uu=exp(-tt);
yy1=polyval(p,uu);
z1=polyval(p,u);
wucha1=sqrt(sum((z1-y).^2))
v=t.*u;
q=polyfit(v,y,1);
vv=tt.*uu;
yy2=polyval(q,vv);
z2=polyval(q,v);
wucha2=sqrt(sum((z2-y).^2))
figure(1);
plot(t,y,'+',tt,yy1,t,z1,'o');
figure(2);
plot(t,y,'+',tt,yy2,t,z2,'o');
Figure1的图
Figure2的图
12．计算下列定积分：
解：
第一个：
建立m文件：
function f=jifen1(x)
f=exp(-2*x);
在命令窗口输入：
[s,n]=quad(@jifen1,0,2)
得到结果：
s =
0.4908
n =
25
第二个：
x=0:0.01:2;
y=exp(2*x);
trapz(x,y)
得到结果：
ans =
26.8000
第三个：
x=-1:0.01:1;
y=x.^2-3*x+0.5;
trapz(x,y)
得到结果：
ans =
1.6667
13．微分方程组
当t=0时，x1(0)=1，x2(0)=-0.5，求微分方程t在[0，25]上的解，并画出相空间轨道图像。
解：
14．设通过测量得到时间t与变量y的数据：
t=[0 0.3 0.8 1.1 1.6 2.3];
y=[0.5 0.82 1.14 1.25 1.35 1.41];
分别采用多项式： y=a0+a1t+a2t2
和指数函数 y=b0+b1e^t+b2te^t
进行拟合，并计算均方误差、画出拟合效果图进行比较。
解：
t=[0 0.3 0.8 1.1 1.6 2.3];
y=[0.5 0.82 1.14 1.25 1.35 1.41];
tt=0:0.01:2.3;
a=polyfit(t,y,2)
yy1=polyval(a,tt);
z1=polyval(a,t);
wucha1=sqrt(sum((z1-y).^2))
B=[ones(size(t')) exp(-t)' ( t.*exp(-t))'];
b=B\y'
yy2=b(1)+b(2)*exp(-tt)+b(3)*tt.*exp(-tt);
z2=b(1)+b(2)*exp(-t)+b(3)*t.*exp(-t); wucha2=sqrt(sum((z2-y).^2))
figure(1);
plot(t,y,'+',tt,yy1,t,z1,'o');
figure(2);
plot(t,y,'+',tt,yy2,t,z2,'o');
得到结果：
a =
-0.2346 0.9134 0.5326 wucha1 =
0.0720
b =
1.4132
-0.9136
0.3808
wucha2 =
0.0026
Figure1的图像
Figure2的图像
15．观察函数： y=e^x-1.5cos(2*pi*x)
在区间[-1,1]上的函数图像，完成下列两题：
（1）用函数fzero求解上述函数在[-1,1]的所有根，验证你的结果；
（2）用函数fminbnd求解上述函数在[-1,1]上的极小、极大、最小和最大值，在函数图像
上标出你求得的最小值点作出验证。
注：可以用help fzero命令查看fzero的调用格式，fzero典型的调用方法是：fzero(@myfun,x0) %返回函数myfun在x0附近的根；
fminbnd典型的调用方法是：
fminbnd(@myfun,x1,x2) %返回函数myfun在区间[x1,x2]上的最小值。
练习2 气象观察站调整问题
程序代码：
a=[272.6,324.5,158.6,412.5,292.8,258.4,334.1,303.2,292.9,243.2,159.7,331.2;251.6,287.3,349.5,2 97.4,227.8,453.6,321.5,451.0,446.2,307.5,421.1,455.1;192.7,433.2,289.9,366.3,466.2,239.1,357.4, 219.7,245.7,411.1,357.0,353.2;246.2,232.4,243.7,372.5,460.4,158.9,298.7,314.5,256.6,327.0,296. 5,423.0;291.7,311.0,502.4,254.0,245.6,324.8,401.0,266.5,251.3,289.9,255.4,362.1;466.5,158.9,22 3.5,425.1,251.4,321.0,315.4,317.4,246.2,277.5,304.2,410.7;258.6,327.4,432.1,403.9,256.6,282.9,3 89.7,413.2,466.5,199.3,282.1,387.6;453.4,365.5,357.6,258.1,278.8,467.2,355.2,228.5,453.6,315.6, 456.3,407.2;158.5,271.0,410.2,344.2,250.0,360.7,376.4,179.4,159.2,342.4,331.2,377.7;324.8,406. 5,235.7,288.8,192.6,284.9,290.5,343.7,283.4,281.2,243.7,411.1]；
y=std(a,0,1)
输出结果：
y =
Columns 1 through 6
100.2660 80.9270 108.2444 63.9747 94.1034 94.2002
Columns 7 through 12
38.0479 85.0735 106.4092 57.2472 86.5136 36.8299
数学建模matlab，lingo编程基础题1
请用matlab或者lingo给出如下问题的程序：1）已知非线性整数规划为：
Max z=8*x1^2-x2^2+3*x3^2+4*x4^2+2*x5^2-8*x1 s.t.
0<=xi<=99(i=1,2, (5)
x1+x2+x3+x4+x5<=400
x1+2*x2+2*x3+x4+6*x5<=800
2*x1+x2+6*x3<=800
第1章MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ 答：MATLAB语言有一下特点： 1、起点高。 2、人机界面适合科技人员。 3、强大而简易的作图功能。 4、智能化程度高。 5、功能丰富，可扩展性强。 1.2 MATLAB系统由那些部分组成？ 答：MATLAB系统主要由开发环境、MATLAB数学函数库、MMATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 答：在MATLAB操作桌面上有五个窗口。在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口。在独立窗口的Desktop菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。
1.5 如何启动M文件编辑/调试器？ 答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？ 答：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器。 （2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到
一、填空题 1、求可逆矩阵A的逆矩阵的指令是 inv(A) 2、A=[0 2 3 0],B=[4 5 0 0],则X=A&B的显示结果为 0 1 0 0 3、求矩阵A的转置矩阵的指令是： A’ 4、x=[10 -20 30]; a=sign(x);显示结果为 1 -1 1 5、Notebook是Matlab的文字处理工具，启动Notebook有两种方法，即一 种是在 Word 中启动，另一种是在 Matlab 中启动。 6、Matlab7.3提供了非常强大的帮助系统，包括帮助命令、帮助 窗口、HTML格式帮助、pdf格式帮助和帮助演示等。 7、clc命令和clear命令的区别是 clc清除命令窗口所有显示内容； Clear清除所有显示的变量 . 8、命令窗口中的type命令的功能为显示文件的内容 . 9、matlab7.3提供了三种程序调试的工具，包括直接检测调试、专用调试命 令和使用 M文件编辑/调试器。 10、matlab除了最常用的主函数和子函数外，还提供了嵌套函数、私有函数、 重载函数和匿名函数。 11、在matlab7.3中，创建句柄图形对象的命令格式为____________ h_obj=funname(‘PropertyName’,PropertyValue,……)_ 12、计算5*e3*|cos175°|语句为 5*exp(3)*abs(cos(175*pi/180)). 13、创建随机矩阵的指令用rand,创建稀疏矩阵用sparse 14、为图形添加网格使用grid on，使用legend命令添加图例。 1、MATLAB常用操作界面包括命令窗口、工作空间窗口（浏览器）、 命令历史窗口、当前目录窗口、内存数组编辑器、M文件编辑/调试器、 帮助导航/浏览器、图形窗口等。 2、MATLABInf或inf表示无穷大、NaN或nan表示不是一个数、nargout表 示函数输出宗量数目。 3、MATLAB中逗号主要用作要显示计算结果的指令与其后指令的分隔； 用作输入量与输入量之间的分隔符；用作数组元素分隔符号。 4、工作空间浏览器主要用于内存变量的查阅、保存和编辑。 5、MATLAB实现将全下标转换为单下标的指令为Sub2ind、据单下标换算出全 下标的指令为Ind2sub。 6、二维数组的标识有“全下标”标识、“单下标”标识、“逻辑1”标识。 7、在一个元胞数组A中寻访第2行第3列元胞元素用A(2,3)；寻访数组第 2行第3列元胞中的内容用A{2,3}。 8、MATLAB中clf用于清除图形窗、clc用于清除指令窗中显示内容、clear 用于清除MATLAB工作空间中保存的变量。 1.已知A=[0 9 6;1 3 0];B=[1 4 3;1 5 0];写出下列各指令运行的结果。 A & B 的运行结果ans= ； A ./ B 的运行结果ans= 。 2. 产生4阶全0方阵的命令为；产生3阶全1方阵的命令为。 3. A=rand(2,5);b=size(A);c=length(A)；则b和c的值分别为和。
Matlab 基础练习题 常量、变量、表达式 1、 MATLAB 中，下面哪些变量名是合法的？（ ） （A ）_num （B ）num_ （C ）num- （D ）-num 2、 在MA TLAB 中，要给出一个复数z 的模，应该使用（ ）函数。 （A ）mod(z) （B ）abs(z) （C ）double(z) （D ）angle(z) 3、 下面属于MATLAB 的预定义特殊变量的是？（ ） （A ）eps （B ）none （C ）zero （D ）exp 4、 判断：在MA TLAB 的内存工作区中，存放一个英文字符 'a' 需要占用1个字节，存放 一个中文字符‘啊’需要占用2个字节。（ 错，都是2个字节 ） 5、 判断：MA TLAB 中，i 和j （ 对 ） 6、 判断：MA TLAB 中，pi 代表圆周率，它等于3.14。（ 错，后面还有很多位小数 ） 7、 在MA TLAB 中，若想计算的5 1)3.0sin(21+= πy 值，那么应该在MA TLAB 的指令窗中 输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。 8、 在MA TLAB 中，a = 1，b = i ，则a 占_8__个字节，b 占_16_个字节，c 占________字 节。 9、 在MA TLAB 中，inf 的含义是__无穷大__，nan 的含义是__非数（结果不定）___。 数组 1、 在MA TLAB 中，X 是一个一维数值数组，现在要把数组X 中的所有元素按原来次序 的逆序排列输出，应该使用下面的（ ）指令。 （A ）X[end:1] （B ）X[end:-1:1] （C ）X (end:-1:1) （D ）X(end:1) 2、 在MA TLAB 中，A 是一个字二维数组，要获取A 的行数和列数，应该使用的MATLAB 的命令是（ ）。 （A ）class(A) （B ）sizeof(A) （C ）size(A) （D ）isa(A) 3、 在MATLAB 中，用指令x=1:9生成数组x 。现在要把x 数组的第二和第七个元素都 赋值为0，应该在指令窗中输入（ ） （A ）x([2 7])=(0 0) （B ）x([2,7])=[0,0] （C ）x[(2,7)]=[0 0] （D ）x[(2 7)]=(0 0) 4、 在MA TLAB 中，依次执行以下指令：clear;A=ones(3,4); A(:)=[-6:5];这时， 若在指令窗中输入指令b=A(:,2)'，那么，MATLAB 输出的结果应该是（ ） （A ）b = -3 -2 -1 （B ）b = -2 -1 0 1 （C ）b = -5 -1 3 （D ）b = -5 -2 1 4 5、 在MA TLAB 中，A = 1:9，现在执行如下指令L1 = ~(A>5)，则MATLAB 的执行结果应 该是L1 =___ 1 1 1 1 1 0 0 0 0___。
matlab实现中值滤波去除脉冲噪声matlab小程序（图像处理）2010-04-1612:58:44阅读8评论0字号：大中小 实验原理：中值滤波器是将领域内像素灰度的中值代替该像素的值,对处理脉冲噪声(椒盐噪声)非常有效。为了对一幅图像上的某个点进行中值滤波处理，必须先将掩模内欲求的像素及其领域的像素值排序，确定出中值，主要功能是使拥有不同灰度的点看起来更接近于它的邻近值。 程序说明：函数名为mid(pic_name,s)的函数，其中参数pic_name为读入的图像，s为掩模矩阵的边长，由用户自行决定。 实验说明：随着掩模矩阵的变大，我们可以看到脉冲噪声去除得更加理想，但同时图像会变得更模糊，因为各点像素与其邻域更为接近，因此，进行中值滤波时选择一个适合的掩模矩阵十分重要。另外，我们看到图像的边界处出现了黑色的斑点，这是由于我采用了0来直译边界，这种影响可用镜像反射方式对称地沿其边界扩展来减弱。 另附：其实本实验可以完全由matlab中的函数median或medfilt2简单实现，此处写出内部处理过程，主要是为了让大家理解中值滤波的具体处理过程。 程序源代码： function mid(pic_name,s) close all; s=double(s); X=imread(pic_name); Y1=imnoise(X,'salt&pepper',0.2);%对读入的图像加脉冲噪声 figure; imshow(uint8(Y1)); Y1=double(Y1); [m,n]=size(X); s2=round(s/2); s3=round(s*s/2);%中值像素点的位置
matlab练习题和答案 控制系统仿真实验 Matlab部分实验结果 目录 实验一 MATLAB基本操 作 ............................................................................................ 1 实验二 Matlab编程 .................................................................................................... 5 实验三Matlab底层图形控制 (6) 实验四控制系统古典分析.............................................................................................12 实验五控制系统现代分析 . (15) 实验六 PID控制器的设 计 ...........................................................................................19 实验七系统状态空间设计.............................................................................................23 实验九直流双闭环调速系统仿真 . (25) 实验一 MATLAB基本操作 1 用MATLAB可以识别的格式输入下面两个矩阵 1233，，1443678，i，，,,2357,,,,2335542，i,,,, A,1357B,,,2675342， i,,3239,,,,189543，，,,1894，， 再求出它们的乘积矩阵C，并将C矩阵的右下角2×3子矩阵赋给D矩阵。赋值完成后，调 用相应的命令查看MATLAB工作空间的占用情况。 A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4];
%1、编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。 function f=NO_1(X); Y(1)=X;k=1; while (X~=1) k=k+1; if (mod(X,2)==0) X=X/2; else X=3*X+1; end Y(k)=X; end plot(Y,'b.') end % 2、编制程序产生一个数组，满足：a1=1，a2=1，从第三个元素开始，每个元素等于前两个元素的和，直到数组的前后两个元素的比值比小于1e-4，并且以红色点线的形式画出这个数组。 clear; A(1)=1; A(2)=1; i=3; Z=1; while (abs(Z)>=1e-4) A(i)=A(i-1)+A(i-2); Z=A(i-1)/A(i); i=i+1; end plot (A,'r.') % 3、编写一个函数，能够产生分段函数。function y=test_3_1(X) if (X<=2) y=*X; elseif (X>6) y=; else y=调用分段函数，绘制曲线。 clear; i=1; for j=0::2; x(i)=j;y(i)=test_3_1(j)*test_3_1(j+2); i=i+1; end plot(x,y) % 4、在2pi周期内画正弦函数曲线，并加注坐标轴标识和标题，然后在3pi/4，pi，5pi/4处分别加入带箭头的说明性文本，最后加注图例。 clear; t = 0:pi/50:2*pi; n = length(t); y = sin(t); plot(t,y,'-bo','linewidth',1) xlabel('X');ylabel('Y'); title('正弦函数曲线'); text,, ' \leftarrow 3pi/4','FontSize',18); text,, ' \leftarrow 3pi/4','FontSize',18); text,,' \leftarrow 3pi/4','FontSize',18); hleg1 = legend('sin(x)'); % 5、A为任意一个n*m矩阵，写程序来计算A 中有多少个零元素，并输出个数。 A=input('输入一个矩阵 A = ') n=length(find(A==0)) % 6、A为任意一个向量，写程序找出A中的最小元素，并且输出这个最小元素。 A=input('输入一个向量 A = '); x=length(A); i=1; y=A(i); while (i
MATLAB期末考试试卷及其参考答案 一、填空题（每空1分，20分） 1、MATLAB常用操作界面包括、工作空间窗口、 、、内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。 2、MATLAB中Inf或inf表示、NaN或nan表示、 nargout表示。 3、MATLAB中逗号主要用作； 用作输入量与输入量之间的分隔符；。 4、工作空间浏览器主要用于内存变量的、和。 5、MATLAB实现将全下标转换为单下标的指令为、据单下标换 算出全下标的指令为。 6、二维数组的标识有、、“逻辑1”标识。 7、在一个元胞数组A中寻访第2行第3列元胞元素用；寻访数组第2行第3列元胞中所存的内容用。 8、4、MATLAB中clf用于、clc用 于、 clear用于。 二、简答题（每题5分，共20分） 1、简述MA TLAB历史指令窗的主要作用。 2、简述空数组的功用。 3、简述MATLAB函数的基本结构。 4、简述绘制二维图形的一般步骤。 三、阅读程序并回答问题（每题4分，共28分） 1、写出下列指令运行结果。 A=zeros(2,4); A(:)=1:8;
s=[2 3 5]; A(s) Sa=[10 20 30]' A(s)=Sa 2、写出下列指令运行结果。 A=reshape(1:16,2,8) reshape(A,4,4) s=[1 3 6 8 9 11 14 16]; A(s)=0 3、写出下列指令运行结果。 A=[1,2;3,4]; B=[-1,-2;2;1]; S=3; A.*B A*B S.*A S*B 4、下面的函数主要完成什么功能？ function f=factor(n) if n<=1 f=1; else f=factor(n-1)*n; end 5、写出下列指令运行结果。 ch=‘ABc123d4e56Fg9’; subch=ch(1:5) revch=ch(end:-1:1)
填空题 1、标点符号; %—用来表示该行为注释行。 可以使命令行不显示运算结果， 2、x 为0 ～4pi ，步长为 的向量，使用命令 x=0:*pi:4*pi 创建。 3、输入矩阵A=[1323?57569 ] ，使用全下标方式用A(2,2)取出元素“-5 ”，使用单下标方式用A(5)取出元素“-5 ”。 4、符号表达式sin(2*a+t)+m 中独立的符号变量为 t 。 5、M 脚本文件和M 函数文件的主要区别是M 脚本文件没有函数定义和M 函数文件有函数定义_______。 6. 设x 是一维数组，x 的倒数第3个元素表示为x(_end-2_) 设y 为二维数组，要删除y 的第34行和48列，可使用命令y(34,:)=[] ； y(:,48)=[]； 7. 将变量x 以Ascii 文本格式存储到文件，应使用命令 save _x ； 8. 在while 表达式, 语句体, End 循环语句中，表达式的值 非零 时表示循环条件为真，语句体将被执行，否则跳出该循环语句； 9.要从键盘读入一个字符串并赋值给变量x ，且给出提示“Who is she ”，应使用命令x=input(‘Who is she’,’s’) ； 10．设A=(a ij )m?n 和B=(b ij )m?n 和C=(c ij )m?n 均为m*n 矩阵，且存在于WorkSpace 中，要产生矩阵D=(a ij ?c ij b ij c ij ?)m?n ，可用命令D=(A-C)/B.^C ， 计算|(A ′B )?1| 可用命令det(inv(A’*B) 11. 在MATLAB 命令窗口中的“>>”标志为MATLAB 的 命令行 提示符，“│”标志为 输入 提示符。 12.已知A=[1 2 3；4 5 0；7 8 9]；B=[1 0 3；1 5 0；0 1 2]；写出下列各指令运行的结果。 A+B ； A.*B ； A==B ； ans= [2,2,6;5,10,0;7,9,11] ans= [1,0,9;4,25,0;0,8,18] ans= [1,0,1;0,1,1;0,0,0] 13.已知A 是矩阵，求A 的对角矩阵函数是 diag ； 求A 的下三角矩阵函数是 tril 。 的程序文件和Simulink 模型文件的扩展名分别是.m 、.mdl 最基本的绘图函数为plot() 16. A=[1,2,3;4,5,6]; A(:,[1,3])=[];A= [2；5] 17. fix= -1 , round= -2 . 18.已知A=[1 3;4 6];C=[A,eye(sizeA);A ,ones(sizeA)] 则Ｃ＝ [1 3 1 0；4 6 0 1；1 3 1 1；4 6 1 1]
MATLAB 学习指南 第六章．编写与调用函数 在这一章中，我们讨论如何用多源代码文件来构造一个程序。首先，解释代码文件在MATLAB中如何工作。在编译语言中，例如FORTRAN，C ，或C++，代码被存储在一个或多个源文件中，在进行编译的时候，这些源文件组合在一起 形成了一个单独的可执行文件。作为一种解释型语言，MATLAB以一种更广泛的方式来处理多个源文件。MATLAB代码被放入带有扩展名.m的ASCII文件（或称m-文件）中。MATLAB 6 有一个集成字处理与调试应用程序，尽管会用到其它编辑程序如vi或emacs，集成字处理与调试应用程序仍是编译m-文件的首选程序。 有两种不同的m-文件。一种是脚本文件，它是一种最简单的文件，仅仅将MATLAB中的指令收集在一起。当在交互提示符处输入文件名执行脚本文件时，MATLAB在m-文件内读取并执行指令，就好像指令是我们输入的。而且，似乎我们能够削减m-文件的内容并将削减过的内容传到MATLAB指令窗口中。这种m-文件的用法将在6.1节中给予概述。 在6.2节中要讨论的第二种m-文件包含一个单一函数，此函数名与此m-文件名相同。这种m-文件包含一段独立的代码，这段代码具有一个明确规定的输入/输出界面；那就是说，传给这段代码一列空变量arg1，arg2，…，这段独立代码就能够被调用，然后返回输出值out1，out2，…。一个函数m-文件的第一个非注释行包含函数标头，其形式如下： 此m-文件以返回指令结束，将执行程序返回到函数被调用的位置。或者在交互指令提示符处或者在另一个m-文件内，无论何时用下列指令调用函数代码，函数代码都将被执行。 输入映射到空变量：arg1=var1，arg2=var2，等等。在函数主体内，输出值被分配给了变量out1，out2，等等。当遇到返回值时，当前值out1，out2，…在函数被调用处被映射到变量outvar1，outvar2，…。在用可变长度自变量和输出变量列表编写函数时，MATLAB允许更多的自由。例如，也可以使用下列指令来调用函数。 在此情况下，仅返回一个单一输出变量，这个变量在出口处包含函数变量out1的值。输入和输出自变量可能是字符串，数值，向量，矩阵，或者更高级的数据结构。 为什么使用函数呢？因为从每门计算机科学课程中可知，把一个大的程序分割 成多个可以单独执行一个被明确规定的和被注释过的任务的小程序会使大程序 易读，易于修改，不易于出错。在MATLAB中，先为程序编写一个主文件，或者是一个脚本文件或者更好的话，是一个能够返回一个单一整数的函数m-文件（返回1表示程序执行成功，0表示不完全程序执行，负值表示出现运行误差），这个主文件是程序的进入点。通过把m-文件当作函数来调用，此程序文件可以
第1章 MATLAB 基础 MATLAB 操作桌面有几个窗口如何使某个窗口脱离桌面成为独立窗口又如何将脱离出去的窗口重新放置到桌面上 答：在MATLAB 操作桌面上有五个窗口。在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock 按钮，点击Undock 按钮就可以使该窗口脱离桌面成为独立窗口。在独立窗口的Desktop 菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。 如何启动M 文件编辑/调试器 答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M 文件编辑/调试器将被启动。在命令窗口中键入edit 命令时也可以启动M 文件编辑/调试器。 存储在工作空间中的数组能编辑吗如何操作 答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 命令历史窗口除了可以观察前面键入的命令外，还有什么用途 答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M 文件中。 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别 答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path 菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB 运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 在MATLAB 中有几种获得帮助的途径 答：（1）帮助浏览器：选择view 菜单中的Help 菜单项或选择Help 菜单中的MATLAB Help 菜单项可以打开帮助浏览器。 （2）help 命令：在命令窗口键入“help ” 命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息。 （3）lookfor 命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数。 （4）模糊查询：输入命令的前几个字母，然后按Tab 键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor 和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help 命令显示详细信息。 第2章 MATLAB 矩阵运算基础 在MATLAB 中如何建立矩阵??? ???194375，并将其赋予变量a 答：在Command Window 窗口输入操作: >> a=[5 7 3;4 9 1] 有几种建立矩阵的方法各有什么优点 答：（1）直接输入法，如a=[1 2 3 4]，优点是输入方法方便简捷； （2）通过M 文件建立矩阵，该方法适用于建立尺寸较大的矩阵，并且易于修改； （3）由函数建立，如y=sin(x)，可以由MATLAB 的内部函数建立一些特殊矩阵； （4）通过数据文件建立，该方法可以调用由其他软件产生数据。 在进行算术运算时，数组运算和矩阵运算各有什么要求 答：进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数。 数组运算和矩阵运算的运算符有什么区别 答：在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b 为矩阵乘，a.*b 为数组乘。 计算矩阵??????????897473535与??? ?? ?????638976242之和。 答： 求??? ???+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。 答： 计算??????=572396a 与?? ? ???=864142b 的数组乘积。 答：
MATLAB试卷 一、选择题（每空2分，总共20分） 1.下列哪个变量的定义是不合法的（） (A) abcd-3 (B) xyz_3 (C ) abcdef (D) x3yz 2.下列哪条指令是求矩阵的行列式的值（） (A)inv (B) diag (C ) det (D) eig 3.在循环结构中跳出循环，执行循环后面代码的命令为（） (A) return (B) break (C) continue (D) keyboard 4.已知a=2:2:8, b=2:5，下面的运算表达式中，出错的 为（） (A) a'*b (B) a .*b (C) a*b (D) a-b 5.用round函数四舍五入对数组[2.48 6.39 3.93 8.52]取整，结果为（） (A) [2 6 3 8] (B) [2 6 4
8] (C) [2 6 4 9] (D) [3 7 4 9] 6. 下面的程序执行后array的值为( ) for k=1:10 if k>6 break; else array(k) = k; end end (A) array = [1, 2, 3, 4, 5, 6] (B) array = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10] (C) array =6 (D) array =10. 7.下列关于脚本文件和函数文件的描述中不正确的是( ) A. 去掉函数文件第一行的定义行可转变成脚本文件; B. 函数文件可以在命令窗口直接运行; C. 脚本文件可以调用函数文件; D. 函数文件中的第一行必须以function开始;
8.对应MATLAB提供的绘制平面曲线、离散序列、三维曲线、三维网格曲线、等高线图的函数为( ) A. plot, stem, mesh, plot3, surf B. plot, stem, plot3, mesh, contour C. plot, hist, mesh, plot3, contour D. plot, hist, plot3, meshc, contour 9.设X=[0:0.5*pi:pi], Y=cos(X)且Z=sin(X)，那么Y.*Z, Y*Z’和cat(1,Y,Z)的结果分别为( ) A. [0 0 0] 0 [1 0 -1; 0 1 0] B. [0 0 0] 0 [1 0 -1 0 1 0] C. [0 0 0] [0 1 0; 0 0 0; 0 -1 0] [1 0 -1; 0 1 0] D. [0 1 0; 0 0 0; 0 -1 0] [0 0 0] [1 0 -1 0 1 0] 10.使用下列哪一条指令可以将图形窗体分割成二行三列，并且将第二行第二列的绘图区域设置为当前的绘图区域（） A.subplot(3,2,2) B.subplot(2,3,2) C.subplot(3 ,2,5) D.subplot(2,3,5)
matlab小程序：将txt中十六进制数转为十进制输出function htod(filename) clc [n]=textread(filename,'%2c'); [a b]=size(n) m=zeros(a,b); mm=zeros(a,1); for i=1:a for j=1:b switch n(i,j) case{'0'}m(i,j)=0; case{'1'}m(i,j)=1; case{'2'}m(i,j)=2; case{'3'}m(i,j)=3; case{'4'}m(i,j)=4; case{'5'}m(i,j)=5; case{'6'}m(i,j)=6; case{'7'}m(i,j)=7; case{'8'}m(i,j)=8; case{'9'}m(i,j)=9; case{'A'}m(i,j)=10; case{'B'}m(i,j)=11; case{'C'}m(i,j)=12; case{'D'}m(i,j)=13; case{'E'}m(i,j)=14; case{'F'}m(i,j)=15; otherwise m(i,j)=nan; end end end %m for i=1:a for j=1:b mm(i)=mm(i)+m(i,j)*16^(j-1); end end %mm [a b]=size(mm); size_mm=a mmm=mm'; savefile='C:\Documents and Settings\Administrator\桌面\test.txt'; fid=fopen(savefile,'w');
fprintf(fid,'%4d',mmm) fclose(fid); matlab如何读取二进制、十六进制txt文档 发现matlab如何读取十六进制的和二进制的txt文章不多。今天刚想了一种方法，所以在这里小结一下，所以matlab中文论坛共享一下，没有参考其他的文章哦，觉得好用就帮顶，不好用提意见。 原帖地址https://www.doczj.com/doc/b511257639.html,/thread-23226-1-1.html 本方法同样适合读取十六进制和二进制以外的其他进制文件， txt使用一个最简单的命令就可以读取textread这是一个十分有用，简便的函数（对于fopen fscanf而言） 读取二进制txt文件： 假如txt文档中内容为00010010001101001000，保存在pin.txt文档中 使用a=textread('pin.txt','%s')' a= '0001''0010''0011''0100''1000' 可以看到数据保存为了char格式。 使用bin2dec b=bin2dec(a)' b= 12348 可以看到成功地转换成了十进制文件。 十六进制文件： 00010010001101001000A B C AA a=textread('pin.txt','%s')' a= '0001''0010''0011''0100''1000''A''B''C''AA' 可以看到成功读取了文件。 b=hex2dec(a)' b= 11617256409610 1112170 读取完毕。 小结：本方法以简单使用方便的方法读取二进制、十六进制的txt文档，欢迎大家提出意见
“”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> ('(x)-3*x^2',0) = -2*(-1/6*3^(1/2)) -2*(-11/6*3^(1/2)) -2*(1/6*3^(1/2)) 3、求解下列各题： 1）30 sin lim x x x x ->- >> x;
>> (((x))^3) = 1/6 2） (10)cos ,x y e x y =求 >> x; >> ((x)*(x),10) = (-32)*(x)*(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> x; >> ((((x^2),0,1/2)),17) =
0.54498710418362222 4）4 2 254x dx x +? >> x; >> (x^4/(25^2)) = 125*(5) - 25*x + x^3/3 5）求由参数方程arctan x y t ??=? =??dy dx 与二阶导 数22 d y dx 。 >> t; >> ((1^2))(t); >> ()() = 1
6）设函数(x)由方程e所确定，求y′(x)。>> x y; *(y)(1); >> ()() = (x + (y)) 7） sin2 x e xdx +∞- ? >> x; >> ()*(2*x); >> (y,0) = 2/5
8） 08x =展开（最高次幂为） >> x (1); taylor(f,0,9) = - (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 1 9） 1sin (3)(2)x y e y =求 >> x y; >> ((1)); >> ((y,3),2) =
MATLAB 期末上机考试试题带答案版 姓名： 学号： 成绩： 1.请实现下图： 50 100 150 200 250 -1 -0.8-0.6-0.4-0.200.20.40.60.81x y y=sin(x) x=linspace(0,8*pi,250); y=sin(x); plot(x,y) area(y,-1) xlabel('x') ylabel('y') title('y=sin(x)') 2.请实现下图： x=linspace(0,2*pi,100); y1=sin(x); subplot(2,2,1) plot(x,y1,'k--') grid on xlabel('x') ylabel('y') title('sin(x)') legend('y=sin(x)') y2=cos(x);
subplot(2,2,2) plot(x,y2,'r--') grid on xlabel('x') ylabel('y') title('cos(x)') legend('y=cos(x)') y3=tan(x); subplot(2,2,3) plot(x,y3,'k-') grid on xlabel('x') ylabel('y') title('tan(x)') legend('y=tan(x)') y4=cot(x); subplot(2,2,4) plot(x,y4) grid on xlabel('x') ylabel('y') title('cot(x)') legend('y=cot(x)') 3.解方程组： a=[3 2 1;1 -1 3;2 4 -4]； b=[7;6;-2] ； x=a\b 4.请实现下图： 2 4 6 8 10 12 14 -1-0.8-0.6-0.4-0.200.20.40.60.81x y '自己的姓名' sin(x)sin(2*x) x=linspace(0,4*pi,1000); y1=sin(x);
第一套 华东交通大学2015-2016学年第一学期Matlab期末考试 一、填空题（4*6=24分） 1、在MATLAB命令窗口中的“>>”标志为MATLAB的提示符，“│”标志为提示符。 2、符号表达式sin(2*a+t)+m 中独立的符号变量为__ _。 3、在通常情况下，左除x=a\b是的解，右除x=b/a是的解，一般情况下，。 4、为了使两个plot的图形在同一个坐标显示，可以使用_____命令进行图形保持；可以使用_ __命令为图形添加网格。 5、倘若要是对x进行赋值，从5到25，间隔是0.1,为；倘若要是对x进行赋值，从2到20，中间间隔100个点,为 6、A=[1,2;3,4];B=[1,0;0,1]; A*B= _ , A.*B=_ _； 二、选择题（4*6=24分） 7、如果x=1: 2 : 8,则x(1)和x(4)分别是( ) A．1，8 B．1, 7 C．2, 8 D．2, 7 8、运行如下程序后，输入９回车，命令窗口（command windows）显示的结果为( ) c=input('请输入一个字符','s'); if c>='A' & c<='Z' disp(setstr(abs(c)+abs('a')-abs('A'))); elseif c>='a'& c<='z' disp(setstr(abs(c)- abs('a')+abs('A'))); elseif c>='0'& c<='9' disp(abs(c)-abs('0')); else disp(c); end A. 9 B. 8 C. 7 D. 6 9、MATLAB表达式2*2^3^2的结果是() A．128 B．4096 C. 262144 D．256 10、在循环结构中跳出循环，执行循环后面代码的命令为( ) (A) return (B) break (C) continue (D) keyboard 11、在图形指定位置加标注命令是（） A. title(x,y,?y=sin(x)?); B. xlabel(x,y,?y=sin(x)?); C. text(x,y,?y=sin(x)?); D. legend(x,y,?y=sin(x)?); 12、下列哪个变量的定义是不合法的（）。 (A) abcd-3 (B) xyz_3 (C) abcdef (D) x3yz 三、程序题（52分） 13（10分）、在同一窗口绘制6个子图，已知t=0.1:0.1:2*pi，每个子图分别绘制（其中y1=sin(t),y2=cos(2t),y3=sin(t)+cos(t),y4=3t, y5=2ln(t),y6=et）； 14（10分）、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)： price<200 没有折扣 200≤price<500 3%折扣 500≤price<1000 5%折扣 1000≤price<2500 8%折扣 2500≤price<5000 10%折扣
自己整理MATLAB知识 1入门 例1-1 绘制正弦曲线和余弦曲线。 x=[0:0.5:360]*pi/180; plot(x,sin(x),x,cos(x)); 例1-2 求方程3x4+7x3+9x2-23=0的全部根。 p=[3,7,9,0,-23]; %建立多项式系数向量 x=roots(p) %求根 例1-3 求积分 quad('x.*log(1+x)',0,1) %‘里是被积函数’0,1分 别是积分上下限 例1-4 求解线性方程组。 a=[2,-3,1;8,3,2;45,1,-9]; %方程左面系数 b=[4;2;17]; %方程右面系数 x=inv(a)*b %也可是x=a\b的形式 例1-5 水仙花 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 例1-6 已知，当n=100时，求y的值。程序如下： y=0; n=100; for i=1:n y=y+1/(2*i-1); end y 例1-7 求[100，200]之间第一个能被21整除的整数 for n=100:200 if rem(n,21)~=0 continue end break end n 例1-8 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求[1,500]之间的全部完数。for m=1:500 s=0; for k=1:m/2
第1章MATLAB基础 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 答：在MATLAB操作桌面上有五个窗口。在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口。在独立窗口的Desktop菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。 1.5 如何启动M文件编辑/调试器？ 答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选
定的内容拷贝到剪贴板中、将选定内容直接拷贝到M 文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path 菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB 运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB 中有几种获得帮助的途径？ 答：（1）帮助浏览器：选择view 菜单中的Help 菜单项或选择Help 菜单中的MATLAB Help 菜单项可以打开帮助浏览器。 （2）help 命令：在命令窗口键入“help ” 命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息。 （3）lookfor 命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数。 （4）模糊查询：输入命令的前几个字母，然后按Tab 键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor 和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help 命令显示详细信息。 第2章 MATLAB 矩阵运算基础 2.1 在MATLAB 中如何建立矩阵?? ? ? ??194375，并将其赋予变量a ？ 答：在Command Window 窗口输入操作: