Approximation bis n = 9

Daten: Magnetisierungskennlinie für Eisenkern Dbl: \(TraditionalForm\`H\&^\_μ\)(\(TraditionalForm\`B\&^\))

daten1 = {{0.05, 10}, {0.125, 20}, {0.2, 24}, {0.25, 27}, {0.3, 
            30}, {0.4, 39}, {0.5, 49}, {0.55, 55}, {0.6, 62.5}, {0.65, 
            73}, {0.7, 87}, {0.75, 100}, {0.8, 120}, {0.85, 140}, {0.9, 
            160}, {0.95, 190}, {1. , 228}, {1.05, 275}, {1.1, 335}, {1.15, 
            410}, {1.2, 500}, {1.25, 615}};

bildd = 
        ListPlot[daten1, 
          PlotStyle -> {PointSize[0.015], Hue[0]}];

Approximation Magnetisierungskennlinie Dbl

Hμ9[B_] = Fit[daten1, {B, (B^3), (B^5), (B^7), (B^9)}, B]

bilda9 = 
        Plot[Hμ9[B], {B, 0, 1.25}, 
          PlotStyle -> Hue[0.7]];

Show[bildd, bilda9];

Zeitfunktion mit Bmax

Hμ9Bmaxzeit[Bmax_] = Hμ9[Bmax Sin[ω t]]

Hμ9Bmaxzeitharmon[Bmax_] = 
      TrigReduce[Hμ9Bmaxzeit[Bmax]] // Simplify

Zeitfunktion mit (Subscript[B, max])= 1 T

Sät10 = Hμ9Bmaxzeitharmon[1]

Zeitfunktion mit (Subscript[B, max])= 1.2 T

Sät12 = Hμ9Bmaxzeitharmon[1.2]

Zeitfunktion von (Subscript[H, ν])

Hν[Hkmax_] = Hkmax Cos[ω t]

H9[ω_, Hkmax_, Bmax_] = 
      Hν[Hkmax] + Hμ9Bmaxzeitharmon[Bmax]

H9werte1 = H9[(α/t), 45, 1]

H9werte2 = H9[(α/t), 60, 1.2]

Plot[{H9[(α/t), 45, 1], 
          H9[(α/t), 60, 1.2]}, {α, -2 π, 2 π}, 
        PlotStyle -> {Hue[0], Hue[0.6]}];

Hμ9zeit[ω_, Bmax_] = Hμ9Bmaxzeitharmon[Bmax]

Hμ9zeit[(α/t), 1]

bildp = 
        ParametricPlot[{1.2 Sin[α], 
            Hμ9zeit[(α/t), 1.2]}, {α, 0, 
            2 π}];

Show[bildd, bildp];

bildp = 
        ParametricPlot[{H9[(α/t), 45, 1], 1  Sin[α]}, {α, 
            0, 2 π}];

Aummag1 = 
      Integrate[
        D[1  Sin[α], α] H9[(α/t), 45, 1], {α, 0,  
          2  π}]

Aummag12 = 
      Integrate[
        D[1.2  Sin[α], α] H9[(α/t), 60, 1.2], {α, 
          0,  2  π}]

bildpp = 
        ParametricPlot[{{H9[(α/t), 45, 1], 
              1  Sin[α]}, {H9[(α/t), 60, 1.2], 
              1.2  Sin[α]}}, {α, 0, 2 π}, 
          GridLines -> Automatic];

daten1ge = 
      Table[{daten1[[i, 2]], daten1[[i, 1]]}, {i, 
          Length[daten1]}]

bilddge = 
        ListPlot[daten1ge, 
          PlotStyle -> {PointSize[0.01], Hue[0]}];

Show[bilddge, bildpp];