matlab实习作业

MATLAB实习报告
班级：*********
姓名：*******
学号：********
习题一：输入10个数，求其中最大数和最小数。

要求分别用循环结构和调用MATLAB的max函数、min函数来实现。

解：（1）循环结构方式
M文件：
indata = input(' Please input the number 10 ');
maxdata = indata(1);
mindata = indata(1);
for(i=1:10)
if(indata(i) > maxdata)
maxdata = indata(i);
end
if(indata(i) < mindata)
mindata = indata(i);
end
end
maxdata
mindata
运行结果：
（2）调用函数方式
M文件：
indata = input('Please input the number 10 ');
maxdata = max(indata)
mindata = min(indata)
运行结果：
习题二：求Fibonacci数列
(1)大于4000的最小项。

(2)5000之内的项数
解：（1）M文件：
f=[1,1];
i = 1;
while f(i)+f(i+1)< 4000
f(i+2)=f(i)+f(i+1);
i=i+1;
end
mindata = min(f)
运行结果：
（2）M文件：
f=[1,1];
i = 1;
while f(i)+f(i+1)< 5000
f(i+2)=f(i)+f(i+1);
i=i+1;
end
data = f
运行结果：
习题三：解方程组Ax ＝b ，分别用求逆解法与直接解法求其解。

解：（1）直接解法 M 文件： A=[ 2,-1, 0, 0, 0; -1, 2,-1, 0, 0; 0,-1, 2,-1, 0; 0, 0,-1, 2,-1; 0, 0, 0,-1, 2] b = [1;0;0;0;0] X = A\b
运行结果：
⎥⎥⎥⎥
⎥⎥⎦
⎤
⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎦
⎤
⎢⎢⎢
⎢⎢
⎢⎣⎡--------=00001,21
00
12100012100012
10001
2b A
（2）逆解法
M文件：
A=[ 2,-1, 0, 0, 0;
-1, 2,-1, 0, 0;
0,-1, 2,-1, 0;
0, 0,-1, 2,-1;
0, 0, 0,-1, 2]
b = [1;0;0;0;0]
X = inv(A)*b
运行结果：
习题四：编一个m程序，求N阶方阵A的行列式的值。

解：M文件
det(rand(input('Random input matrix order ')))
运行结果：
习题五：今有多项式P1(x)=x4-2x+1,P2(x)=x2+4x-0.5，要求先求得P(x)=P1(x)+P2(x)，然后计算xi=0.2*i各点上的P(xi)（i=0,1,2,…,5）值。

解：M文件：
P1 = [1,0,0,-2,1]; %P1(x)多项式的系数
P2 = [0,0,1,4,-0.5];%P2(x) 多项式的系数
P = P1 + P2 %P1(x)+P2(x) 多项式的系数
i = 0:5;
Xi = i * 0.2 %Xi
y = polyval(P,Xi) %输出
运行结果：
习题六：试编一个m程序，将一维数组x中的N个数按颠倒的次序重新存储。

如N=5，原来x为：
x=[ 1 3 5 7 9 ]
而经过颠倒处理后x中数据的次序应该为：
x=[ 9 7 5 3 1 ]
解：M文件：
num = input('num is ');
indata = input('indata is')
if mod(num,2)
display('The input is odd');
for i = 1:1:(num-1)/2
temp = indata(i);
indata(i) = indata(num+1-i);
indata(num+1-i) = temp;
end
else
display('The input is even');
for i = 1:1:(num/2)
temp = indata(i);
indata(i) = indata(num+1-i);
indata(num+1-i) = temp;
end
end
display('Reverse the order after the results is');
display(indata);
运行结果：
当输入的数为偶数时：
当输入的数为奇数时：
习题七：在[0 2π]范围内绘制二维曲线图y=sin(x)*cos(5x)。

解： M文件：
x = linspace(0,2*pi,60);
y = sin(x).*cos(5*x)
plot(x,y,'b','linewidth',3);
axis([0 2*pi -1.5 1.5]);
title('sin(x)cos(5x)');
xlabel('x');
ylabel('y');
运行结果：
习题八：把自己姓名中的各个拼音字母（小写）所对应的ASCII 值转换为二进制序列，并用二进制脉冲绘图表示出来
解：M文件：
name = 'wangcunyue'
ASCII = double(a)
c = dec2bin(b,8)%Convert decimal integer to a binary string
c = c'; %??
stem(str2num(c(:)))%Convert string matrix to numeric array
name = 'wangcunyue'
b = double(name)
c = dec2bin(b,8)%Convert decimal integer to a binary string
c = c'; %??
stem(str2num(c(:)))%Convert string matrix to numeric array
运行结果：
心得体会：通过这几天的学习，我对matlab具有的功能震撼了，应用于各个领域。

我很喜欢它的编程风格，它与面向过程的C语言有一定的相似又有很大的不同，它的Simulink仿真平台与我之前接触过的Labview编程风格很相识，我认为也可以把它称之为G语言（图形语言），它的GUI图形界面又与我之前接触的面向对象的VB语言很相识，建立起图形界面很方面，只需要在相应的控件的回调函数来编写的它的相应的操作即可，它的强大功能深深的吸引了我，在资金不足的情况下在网上买了一本精通
matlab的书籍，打算利用空余时间好好学习这个拥有强大魅力的软件。