xianxingnihe.m

来自「数字信号处理课程的有关实验」· M 代码 · 共 43 行

clear all;
clc;
data=[0   498
      0.1 490
      0.2 486
      0.3 480
      0.4 473
      0.5 456
      0.6 452
      0.7 445
      0.8 433
      0.9 417
%       1.0 398
%       1.1 385
%       1.2 370
%       1.3 353
%       1.4 341
%       1.5 330
%       1.6 324
%       1.7 309
%       1.8 299
%       1.9 285
%       2.0 277
  ];
x=data(:,1);
y=data(:,2);
A(:,1)=ones(size(x));A(:,2)=x;
c=A\y;
xx=0:0.01:2;
g=c(1)*ones(size(xx))+c(2)*xx;
c
figure(1),plot(x,y,'o',x,data,'k');axis([0,0.91,410,500]);title('电流变化曲线');
legend('离散电流值','电流变化曲线')
xlabel('电阻烧蚀长度（单位：mm）')
ylabel('电流值（单位：mA)')
figure(2),plot(x,y,'s',xx,g,':');axis([0,0.91,410,500]);
hold on
plot(x,y,'o',x,data,'k');axis([0,0.91,410,500]);
title('函数组合曲线拟合');
ylabel('电流值（单位：mA)')
legend('电流变化曲线','拟合后函数曲线')
xlabel('电阻烧蚀长度（单位：mm）')
ylabel('电流值（单位：mA)')