1 temp=a;
a=b;
b=temp;
2 a=input('a=');
b=input('b=');
if a>=b
fprintf('较大的数是a,a=%g\n',a)
else
fprintf('较大的数是b,b=%g\n',b)
end
3 function [x1,x2 ] = xunlian24_3( a,b,c)
p=b*b-4*a*c;
x1=(-b+sqrt(p))/2*a;
x2=(-b-sqrt(p))/2*a;
disp(x1)
disp(x2)
end
4 %函数文件
function [s c]=xunlian24_4(r)
r=input('请输入圆的半径r=');
s=pi*r*r;c=2*pi*r;
fprintf('面积是s=%.4f\n',s)
fprintf('周长是c=%.4f\n',c)
end
%脚本文件
r=input('r=');
s=pi*r*r;c=2*pi*r;
fprintf('s=%.4f\n',s)
fprintf('c=%.4f\n',c)
5
function xunlian24_5(t)
s=0;
for ii=1:t
s=s+digui(ii);%把ii放到下面那个累乘的函数，累乘完后相加
end
s %相加后输出s
function [ y ] =digui( n )
if n==1 y=1;
else
y=digui(n-1)*n; %对比它小一的数再用一次函数，达到累乘的效果
end
end
end
6 x=input('输入x:');
if x<1
p=x^2;
elseif x>=1&x<2
p=x^2-1;
else
p=x^2-2*x-1;
end
fprintf('p=%g\n',p);
switch p
case 1
1.常见信号的Matlab表示及运算 作业题
1
1、用MATLAB命令函数分别绘出抽样信号、矩形脉冲信号、周期矩形脉冲信号、三角波脉冲信号、三角波周期信号的波形，理解其中参数的意义。
t=-20:0.01:20;
f10=sinc(t);
f20=rectpuls(t,2);
f30=square(t,2);
f40=tripuls(t,0.5);
f50=sawtooth(t,0.5);
subplot(5,1,1);
plot(t,f10);
axis([-10,10,-2,2]);
subplot(5,1,2);
plot(t,f20);
axis([-10,10,-2,2]);
subplot(5,1,3);
plot(t,f30);
axis([-10,10,-2,2]);
subplot(5,1,4);
plot(t,f40);
axis([-10,10,-2,2]);
subplot(5,1,5);
plot(t,f50);
axis([-10,10,-2,2]);
1.常见信号的Matlab表示及运算 作业题
2
t=-20:0.01:20;
f20=rectpuls(t,2);
f21=rectpuls(t,4);
f22=rectpuls(t,8);
f23=rectpuls(t,16);
f24=rectpuls(t,32);
subplot(5,1,1);
plot(t,f20);
axis([-10,10,-2,2]);
subplot(5,1,2);
plot(t,f21);
axis([-10,10,-2,2]);
subplot(5,1,3);
plot(t,f22);
axis([-10,10,-2,2]);
subplot(5,1,4);
plot(t,f23);
axis([-10,10,-2,2]);
subplot(5,1,5);
plot(t,f24);
0 实验一 Matlab使用方法和程序设计
一、实验目的
1、掌握Matlab软件使用的基本方法；
2、熟悉Matlab的数据表示、基本运算和程序控制语句
3、熟悉Matlab绘图命令及基本绘图控制
4、熟悉Matlab程序设计的基本方法调试方法。
二、实验内容：
1、帮助命令
使用help命令，查找 sqrt（开方）函数的使用方法；
2、矩阵运算
（1）矩阵的乘法
已知A=[1 2;3 4]; B=[5 5;7 8];
求A^2*B
（2）矩阵除法
已知 A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3];
求A\B,A/B
（3）矩阵的转置及共轭转置
已知A=[5+i,2-i,1;6*i,4,9-i];
求A.', A'
（4）使用冒号选出指定元素
已知： A=[1 2 3;4 5 6;7 8 9];
求A中第3列前2个元素；A中所有列第2，3行的元素；
3、多项式
（1）求多项式42)(3xxxp的根
（2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]，求矩阵A的特征多项式；
4、基本绘图命令
（1）在同一坐标系中绘制余弦曲线y=cos(t)，y=cos(t-0.25)和正弦曲线y=sin(t-0.5)， t∈[0，2π]
（2）绘制[0，4π]区间上的x1=10sint曲线，并要求：
1）线形为点划线、颜色为红色、数据点标记为加号；
2）坐标轴控制：显示范围、刻度线、网络线；
3）标注控制：坐标轴名称、标题、相应文本；
5、基本程序设计
（1）编写命令文件：计算1+2+⋯⋯+n<2000时的最大n值；
（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。
三、预习要求
利用所学知识，编写实验内容中2到5的相应程序，并写在预习报告上。 1 四、 实验报告要求
第三章
1．要求在闭区间]2,0[上产生具有10个等距采样点的一维数组。试用两种不同的指令实
现。
方法一t1=linspace(0,2*pi,10)
t1=00.69811.39632.09442.79253.49074.18884.88695.58516.2832
方法二t2=0:2*pi/9:2*pi
结果与t1相同
2．由指令rand('state',0),A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5的元素
的位置，分别求出它们的“全下标”和“单下标”。rand('state',0),A=rand(3,5)
[ri,cj]=find(A>0.5);
id=sub2ind(size(A),ri,cj);
ri=ri';
cj=cj';
disp('')
disp('大于0.5的元素的全下标')
disp(['行号',int2str(ri)])
disp(['列号',int2str(cj)])
disp('')
disp('大于0.5的元素的单下标')
disp(id')
A=
0.95010.48600.45650.44470.9218
0.23110.89130.01850.61540.7382
0.60680.76210.82140.79190.1763
大于0.5的元素的全下标
行号132332312
列号112234455
大于0.5的元素的单下标13569111213143．已知矩阵4321A，运行指令B1=A.^(0.5),B2=A^(0.5),可以观察到不同运算方法
所得结果不同。（1）请分别写出根据B1,B2恢复原矩阵A的程序。（2）用指令检验
所得的两个恢复矩阵是否相等。A=[1,2;3,4];
B1=A.^0.5
B2=A^0.5
A1=B1.*B1;
A2=B2*B2;
norm(A1-A2,'fro')
B1=
1.00001.4142
1.73212.0000
B2=
0.5537+0.4644i0.8070-0.2124i
