MATLAB数值运算（6）

实验目的

1）熟悉数据处理与统计方法；
2）掌握常用的数值运算方法；
3）掌握MATLAB中的常用数值运算函数；
4）熟悉运用MATLAB编程实现数值计算。

实验要求

1）实验前，要认真预习实验任务，了解实验目的和实验内容；
2）实验时，要利用MATLAB语言独立编写程序代码，并调试程序使其能正确运行；
3）实验后，按要求编写实验报告，源程序要有适当的注释，以提高程序的可读性。

![在这里插入图片描述]

在这里插入代码片A=[90 89 91 82 75;88 82 85 89 95;79 89 93 82 76;
   91 89 87 79 85;90 79 78 91 82;86 85 84 89 88;
   89 91 79 90 94;75 85 95 90 87;85 78 75 90 89;98 89 87 86 80];
max(A)  %每门课的最高分
min(A)  %每门课的最低分
B=mean(A) %每门课的平均分
D=sum(A,2)%每个学生的各课总分
max(D) %总分最高分
min(D) %总分最低分
mean(D)%学生总分的平均分
% max(A,[],2) %每位学生的最高分
% min(A,[],2) %每位学生的最低分
% mean(A,2) %每位学生的平均分
x=1:10;
plot(x,0*x+B(:,1),'r')
hold on
plot(x, A(:,1), 'r');  %数学成绩 红色 r   
plot(x, A(:,2), 'g');  %语文成绩 绿色 g
plot(x, A(:,3), 'b');  %外语成绩 蓝色 b  
plot(x, A(:,4), 'c');  %物理成绩 青色 c
plot(x, A(:,5), 'k');  %化学成绩 黑色 k   
plot(x,0*x+B(:,2),'b')
plot(x,0*x+B(:,3),'g')
plot(x,0*x+B(:,4),'c')
plot(x,0*x+B(:,5),'r')
hold off

P1=[1 0 1];
P2=[1 3];
P3=[1 1];
s1=conv(P1,conv(P2,P3)) %conv用于多项式P2和P3的乘积
s2=[1 0 2 1];
[q,r]=deconv(s1,s2)  
% deconv求s1与s2商与余项，q是商，r余项
%r与s1的维度相同

a=[1 2 4 0 5];
b=[0 1 2];
c=[1 2 3];
e=a+conv(b,c)  %poly2sym(e)的系数
poly2sym(e)  %x^4 + 2*x^3 + 5*x^2 + 3*x + 10
x=roots(e) %求根
A=[-1 1.2 -1.4;0.75 2 3.5;0 5 2.5];
polyval(e,A)%求A矩阵每个元素对应求e多项式的值

h1=6:2:18; %原测量时间
h2=6.5:2:17.5;  %现测量时间
t1=[18,20,22,25,30,28,24];  %原室内温度
t2=[15,19,24,28,34,32,30];  %原室外温度
T1=interp1(h1,t1,h2,'spline') %T1=interp1(x,Y,xx,'way')，x和Y是原矢量，xx是要求各点位置矢量，way方法，T1是生成的新值
T2=interp1(h1,t2,h2,'spline')
plot(h2,T1,'b')
hold on
plot(h2,T2,'r')
hold off
xlabel('时间(h)')
ylabel('温度（℃）')
legend('室内温度','室外温度','Location','NorthEast')
%legend（）使用格式：legend(string1,string2,string3, ...，n)string1——对应第一条曲线的字符串

syms x
y=sin(x)^3+cos(x)^3;
y1=diff(y); %对y求1阶导
s1=subs(y1,pi/6)  %R = subs(S, new) 利用new的值代替符号表达式S中的默认符号。
s2=subs(y1,pi/4)
s3=subs(y1,pi/3)
s4=subs(y1,pi/2)

 x=fminsearch('(x.^3+cos(x)+x.*log(x))./exp(x)',0,1)  %0到1的最小值

% syms x;
% y=exp(x).*(1+exp(x)).^2;
% I1=int(y,0,log(2));  %进行定积分 算不出来
format long
x=0:0.01:log(2);%选择不同的布长，结果的精度不一
y=exp(x).*(1+exp(x)).^2;
I1=trapz(x,y)  %进行积分

实验结论

第1题最主要是学会画一条直线方法。第2题和第2题主要是代入公式就可以了，特别是学会polyval。第4题学会3次样条插值，学会interp1就可以了。第567题也是带入公式即可。这一章主要了解各个函数的使用就行可以了，相对简单，但有两个画图的，就是对上一个章的巩固，最主要是学会了画直线，这个很重要。这一章对程序要求很高，有利于提高程序编写能力。