1.随机生成一组长度为10000的服从正态分布的随机数,其中均值为10,方差是2,然后检验随机数的性质。

1.随机生成一组长度为10000的服从正态分布的随机数,其中均值为10,方差是2,然后检验随机数的性质。（1）计算其均值u、方差a和中值；（2）分别计算这组随机数在[u-a,u+a]，[u-2a,u+2a]，[u-3a,u+3a]范围内占总数的比例。
2.根据所给数据data.mat，用四种插值法求出[5.5,6.3,7.6,8.7,9.8,10.1]六个位置的函数值，并画出函数曲线和插值点。
question1

%1.随机生成一组长度为10000的服从正态分布的随机数,其中均值为10,方差是2,然后检验随机数的性质。
% （1）计算其均值[图片]、方差[图片]和中值；
% （2）分别计算这组随机数在[图片]，[图片]，[图片]范围内占总数的比例。
clc;
clear;
A=normrnd(10,sqrt(2),1,10000);%生成10行5列均值为0方差为6的正态分布随机数
u=mean(A);%求均值
b=var(A);%求方差
c=median(A);%求中位数
k1=find(A>u-b&A<u+b);
s1=size(k1,2)/size(A,2);
k2=find(A>u-2*b&A<u+2*b);
s2=size(k2,2)/size(A,2);
k3=find(A>u-3*b&A<u+3*b);
s3=size(k3,2)/size(A,2);
fprintf('均值u=%2f 方差b=%f 中位数c=%2f 第一个占比s1=%2f 第二个占比s2=%2f 第三个占比s3=%2f',u,b,c,s1,s2,s3)

question2

syms y1 y2 y3 y4
clc;
clear;
load('data.mat');
A=a;
x=1:201;
y=A(:,4);
c=5.5;
%interp1对sin函数进行分段线性插值，调用interp1的时候，默认的是分段线性插值
z1=interp1(x,y,c,'nearest');
y1=interp1(x,y,x,'nearest');%method是最邻近插值
figure
plot(x,y,'o',x,y1,'-')
title('临近插值')

%临近插值
z2=interp1(x,y,c,'linear');
y2=interp1(x,y,x,'linear');
figure
plot(x,y,'o',x,y2,'r');
title('线性插值')

%球面线性插值
z3=interp1(x,y,c,'spline');
y3=interp1(x,y,x,'spline');
figure
plot(x,y,'o',x,y3,'-')
title('球面插值')

%三次多项式插值法
z4=interp1(x,y,c,'cubic');
y4=interp1(x,y,x,'cubic');
figure
plot(x,y,'o',x,y4,'-');
title('三次多项式插值')
fprintf('(第一种)5.5位置函数值=%2f (第二种)5.5位置函数值=%f (第三种)5.5位置函数值=%2f (第四种)5.5位置函数值=%2f',z1,z2,z3,z4)