大物实验-冲击电流计测量磁滞回线matlab数据处理（matlab学习实践记录）

通过实验测得45组Ii，di，最后计算要得到HB图。matlab这个软件我也刚刚开始探索，对这个页面还是不太熟悉，但是这个数据确实有点多，还是勉强自己速成一下。

公式

实验测得K=0.01，N1=1100，L1=Π*14.5（钢圆环轴线周长）

matlab代码

  clc %清除命令窗口的内容
  clear %清除工作空间的所有变量
  clear all %清除工作空间的所有变量，函数，和MEX文件
  clf %清除当前的Figure
  close %关闭当前的Figure窗口
  close all %关闭所有的Figure窗口

  data1 = Untitled(1:45,1); %插入Ii数据
  data2 = Untitled(1:45,2); %插入di数据
  Ii = data1;
  di = data2;
  di_sum = 0;
  i = 0;

  %套公式计算
  Hi = Ii * 1100 / (3.14 * 0.145);
  for i = 1:1:45
    di_sum = di_sum + di(i);
    Bi(i) = 0.01 * di_sum;
  end

   %画图
   scatter(Hi,Bi,5)
   hold on
   plot(Hi,Bi)
   title('磁滞回线图')
   xlabel('H')
   ylabel('B')
   %打印数据
   fprintf('Bi=%f\r\n',Bi);
   fprintf('Hi=%f\r\n',Hi);

得到磁滞回线图

Bi，Hi数据（按顺序）

只截了一部分。