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)