想用三角函数拟合一组数据结果一直报错感谢大神解答程序就是源程序clear

odayle · 发表于 2020-3-12 13:56

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

想用三角函数拟合一组数据结果一直报错感谢大神解答程序就是源程序clear;
theta=linspace(0,rad2deg(2*pi),181);
y1=[0.512
0.471
0.445
0.408
0.382
0.347
0.33
0.295
0.275
0.2449
0.2222
0.1974
0.1792
0.1573
0.1447
0.1289
0.1155
0.1072
0.0991
0.092
0.0877
0.0841
0.0842
0.0856
0.0869
0.0929
0.1027
0.1064
0.1165
0.1295
0.1447
0.1558
0.1696
0.19
0.21
0.2287
0.2496
0.2719
0.2986
0.3246
0.343
0.372
0.393
0.415
0.444
0.47
0.494
0.513
0.542
0.573
0.6
0.613
0.636
0.652
0.675
0.697
0.71
0.729
0.742
0.76
0.767
0.777
0.785
0.79
0.792
0.799
0.808
0.795
0.794
0.789
0.786
0.785
0.782
0.773
0.766
0.746
0.726
0.709
0.697
0.678
0.658
0.639
0.611
0.597
0.568
0.553
0.522
0.501
0.477
0.45
0.421
0.396
0.373
0.342
0.3196
0.2984
0.2686
0.2559
0.2146
0.206
0.1936
0.169
0.1516
0.1428
0.1234
0.1139
0.1051
0.0962
0.0895
0.0837
0.0802
0.0782
0.081
0.0798
0.0824
0.087
0.0921
0.1017
0.1089
0.1205
0.1345
0.1467
0.16
0.1759
0.1947
0.2144
0.2307
0.2553
0.241
0.3002
0.3268
0.348
0.369
0.399
0.42
0.44
0.466
0.492
0.51
0.535
0.555
0.576
0.592
0.622
0.638
0.648
0.666
0.683
0.703
0.714
0.722
0.749
0.752
0.752
0.759
0.76
0.759
0.754
0.752
0.754
0.751
0.746
0.741
0.728
0.715
0.702
0.684
0.677
0.662
0.635
0.62
0.598
0.586
0.56
0.537
0.513
0.494
0.462
0.443
0.414
0.385]';
% [5000 0 25000 3*pi/55],
a=[0.4 0.1 0.4];
func_sin = @(a,theta) a(1)*(sin((theta)+a(2)))^2 + a(3);
A= lsqcurvefit( func_sin,a,theta,y1);
disp(A)

错误类型：

出错 Untitledexperiment (line 187)
A= lsqcurvefit( func_sin,a,theta,y1);

原因:
Failure in initial user-supplied objective function evaluation. LSQCURVEFIT cannot continue.

House · 发表于 2020-3-12 15:24

fx=@(a,x)a(1)*(sin(a(2)*x+a(3))).^2 + a(4);
SSy = 10.965
a =[0.6987103424 0.01759903346 -0.791168842 0.08655566292]
RSS = 0.0708197085559
MSe = 0.00040238
R^2 = 0.99354
也可用fx=(a(1)+a(2)*x+a(3)*x.^2).*sin(a(4)*x+a(5))+a(6);拟合，效果略好些。

帐号		自动登录	找回密码
密码			注册

想用三角函数拟合一组数据 结果一直报错 感谢大神解答 程序就是源程序clear

EDA365欢迎您登录！

推荐内容 /1

想用三角函数拟合一组数据结果一直报错感谢大神解答程序就是源程序clear