实验8 Matlab综合实验

实验目的：
1、掌握matlab基本命令输入；
2、掌握函数文件的使用方法；
3、掌握命令文件的编写方法；
4、掌握基本数值计算方法与matlab常用工具箱函数的使用。
实验内容：
1.已知，
，利用最小二乘法进行曲线拟合，说明拟合过程中多项式阶数的选择对拟合结果的影响，并给出较为合理的拟合多项式。

x=[0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1];
[~,k]=size(x);
y=[-0.447,1.978,3.28,6.16,7.08,7.34,7.66,9.56,9.48,9.30,11.2];
for n=1:9
    ANSS=polyfit(x,y,n);  
    for i=1:n+1 
       answer(i,n)=ANSS(i);
   end
    x0=0:0.01:1;
    y0=ANSS(1)*x0.^n;
    for num=2:1:n+1     
        y0=y0+ANSS(num)*x0.^(n+1-num);
    end
    subplot(3,3,n)
    plot(x,y,'*')
    hold on
    plot(x0,y0)
end
suptitle('不同次数方程曲线拟合结果，从1到9阶')

2.在同一窗口画出以下信号的曲线，每条曲线使用不同的颜色和点形，并在每个图中加上图例和说明。
（1）
；（2）；
（3）；（4）
注：x轴取值范围自选。

t1=-5:0.01:5;
t2=-0.05:0.001:0.05;
x1=-6:0.001:6;
x2=-0.06:0.001:0.06;
y1=sin(2*t1).*sin(9*t1);
y2=exp(3*sqrt(t2));
y3=(0.8+(3*cos(x1))/(1+x1.^3))*sin(x1);
y4=0.6*exp(-0.5*x2).*sin(10*x2);
subplot(2,2,1),plot(t1,y1,'m-d','linewidth',0.001,'MarkerSize',10);title('第一题')
legend('sin(2t1)sin(9t1)')
xlabel('t1轴')
ylabel('y1轴')
text(0,0,'原式=0')
subplot(2,2,2),plot(t2,y2,'r:*','linewidth',0.001,'MarkerSize',3);title('第二题')
legend('exp(3*sqrt(t2))')
xlabel('t2轴')
ylabel('y2轴')
text(0,1,'原式=1')
subplot(2,2,3),plot(x1,y3,'b:+','linewidth',0.001,'MarkerSize',3);title('第三题')
legend('(0.8+(3*cos(x1))/(1+x1.^3))*sin(x1)')
xlabel('x1轴')
ylabel('y3轴')
text(0,0,'原式=0')
subplot(2,2,4),plot(x2,y4,'k-h','linewidth',0.001,'MarkerSize',10);title('第四题')
legend('0.6*exp(-0.5*x2).*sin(10*x2)')
xlabel('x2轴')
ylabel('y4轴')

3.读入一副图像，说明图像的格式，并将此图像变换为黑白格式，加入随机噪声。

I=imread('D:\a_study\matlab\sy\1.png');
inf=imfinfo('D:\a_study\matlab\sy\1.png');
A=rgb2gray(I);
H=imnoise(A,'speckle',0.1);
imshow(H)