《电机与拖动基础及MATLAB仿真》陈亚爱版课件(1-7章)+代码.zip

clc clear Hdata=[1.38, 1.40, 1.42, 1.44, 1.46, 1.48, 1.50, 1.52, 1.54, 1.56,... 1.58, 1.60, 1.62, 1.64, 1.66, 1.69, 1.71, 1.74, 1.76, 1.78,... 1.81, 1.84, 1.86, 1.89, 1.91, 1.94, 1.97, 2.00, 2.03, 2.06,... 2.10, 2.13, 2.16, 2.20, 2.24, 2.28, 2.32, 2.36, 2.40, 2.45,... 2.50, 2.55, 2.60, 2.65, 2.70, 2.76, 2.81, 2.87, 2.93, 2.99,... 3.06, 3.13, 3.19, 3.26, 3.33, 3.41, 3.49, 3.57, 3.65, 3.74,... 3.83, 3.92, 4.01, 4.11, 4.22, 4.33, 4.44, 4.56, 4.67, 4.80,... 4.93, 5.07, 5.21, 5.36, 5.52, 5.68, 5.84, 6.00, 6.16, 6.33,... 6.52, 6.72, 6.94, 7.16, 7.38, 7.62, 7.86, 8.10, 8.36, 8.62,... 8.90, 9.20, 9.50, 9.80, 10.1, 10.5, 10.9, 11.3, 11.7, 12.1,... 12.6, 13.1, 13.6, 14.2, 14.8, 15.5, 16.3, 17.1, 18.1, 19.1,... 20.1, 21.2, 22.4, 23.7, 25.0, 26.7, 28.5, 30.4, 32.6, 35.1,... 37.8, 40.7, 43.7, 46.8, 50.0, 53.4, 56.8, 60.4, 64.0, 67.8];% 磁场强度数据 Bdata=0.40:0.01:1.69; % 磁感应强度数据 len=length(Hdata); % 计算数组长度 Hmax=Hdata(len); % 提取数组最大值 a=polyfit(Hdata,Bdata,13); % 计算多项式拟合系数a for n=1:151 % 提取拟合数据 Hfit(n)=Hmax*(n-1)/150; Bfit(n)=a(1)*Hfit(n)^13+ a(2)*Hfit(n)^12+ a(3)*Hfit(n)^11+... a(4)*Hfit(n)^10+ a(5)*Hfit(n)^9+ a(6)*Hfit(n)^8+ a(7)*Hfit(n)^7+... a(8)*Hfit(n)^6+ a(9)*Hfit(n)^5+ a(10)*Hfit(n)^4+ a(11)*Hfit(n)^3+... a(12)*Hfit(n)^2+ a(13)*Hfit(n)+ a(14); end plot(Hdata,Bdata, '*') % 对原数据描点绘图 hold on % 保持当前坐标轴和图形 plot(Hfit,Bfit) % 绘制多项式拟合曲线 hold on % 保持当前坐标轴和图形 title('磁化曲线') % 标题'磁化曲线' xlabel('{\itH}/(A/cm)') % x坐标标签“H(A/cm)” ylabel('{\itB}/T') % y坐标标签“B/T”