|
[求助]
1stopt高版本代跑 (微分方程与代数方程参数拟合)已有2人参与
7参数拟合问题,含有微分方程与代数方程,请求高版本代跑一下程序。
Parameter fac1=[1e-5,],fac2=[1e-5,],fac3=[1e-5,],fac4=[1e-5,],fac5=[,-1e-5],fac6=[,-1e-5],fac7=[1e-5,];
Variable t,lamb,lambdif,pb;
ODEFunction lambv'=if(lambdif<=0,(lamb/3/fac7)*(fac1*((lamb/lambv)^(fac2-1)-(lambv/lamb)^(0.5*fac2+1))+fac3*((lamb/lambv)^(fac4-1)-(lambv/lamb)^(0.5*fac4+1))+fac5*((lamb/lambv)^(fac6-1)-(lambv/lamb)^(0.5*fac6+1))),(lamb/(0.001*3*fac7))*(fac1*((lamb/lambv)^(fac2-1)-(lambv/lamb)^(0.5*fac2+1))+fac3*((lamb/lambv)^(fac4-1)-(lambv/lamb)^(0.5*fac4+1))+fac5*((lamb/lambv)^(fac6-1)-(lambv/lamb)^(0.5*fac6+1))));
Function pb=fac1*(lamb^(fac2-1)*lambv^(-fac2)-lamb^(-0.5*fac2-1)*lambv^(0.5*fac2))+fac3*(lamb^(fac4-1)*lambv^(-fac4)-lamb^(-0.5*fac4-1)*lambv^(0.5*fac4))+fac5*(lamb^(fac6-1)*lambv^(-fac6)-lamb^(-0.5*fac6-1)*lambv^(0.5*fac6));
Data;
//0.1Hz
0.00 0.960000 0.012566 0.057751
0.20 0.962433 0.012467 0.063943
0.40 0.964893 0.012172 0.063200
0.60 0.967303 0.011684 0.063611
0.80 0.969577 0.011012 0.065567
1.00 0.971710 0.010166 0.065701
1.20 0.973627 0.009160 0.065639
1.40 0.975377 0.008010 0.064222
1.60 0.976853 0.006733 0.062565
1.80 0.978063 0.005351 0.061493
2.00 0.979000 0.003883 0.058928
2.20 0.979637 0.002355 0.058219
2.40 0.979960 0.000789 0.055017
2.60 0.979960 -0.000789 0.052617
2.80 0.979660 -0.002355 0.047199
3.00 0.979037 -0.003883 0.042338
3.20 0.978137 -0.005351 0.032877
3.40 0.976940 -0.006733 0.025563
3.60 0.975463 -0.008010 0.018672
3.80 0.973743 -0.009160 0.010299
4.00 0.971830 -0.010166 0.004338
4.20 0.969690 -0.011012 -0.004050
4.40 0.967423 -0.011684 -0.010569
4.60 0.965040 -0.012172 -0.019444
4.80 0.962587 -0.012467 -0.026517
5.00 0.960073 -0.012566 -0.031025
5.20 0.957553 -0.012467 -0.039910
5.40 0.955093 -0.012172 -0.044413
5.60 0.952710 -0.011684 -0.050038
5.80 0.950440 -0.011012 -0.053758
6.00 0.948293 -0.010166 -0.059845
6.20 0.946357 -0.009160 -0.062465
6.40 0.944640 -0.008010 -0.063932
6.60 0.943153 -0.006733 -0.064208
6.80 0.941933 -0.005351 -0.066548
7.00 0.940993 -0.003883 -0.063021
7.20 0.940363 -0.002355 -0.059452
7.40 0.940043 -0.000789 -0.055989
7.60 0.940033 0.000789 -0.051224
7.80 0.940323 0.002355 -0.039363
8.00 0.940947 0.003883 -0.025688
8.20 0.941860 0.005351 -0.009293
8.40 0.943073 0.006733 0.005272
8.60 0.944540 0.008010 0.016966
8.80 0.946260 0.009160 0.028023
9.00 0.948193 0.010166 0.034484
9.20 0.950300 0.011012 0.042393
9.40 0.952567 0.011684 0.048518
9.60 0.954957 0.012172 0.050146
9.80 0.957423 0.012467 0.055289
10.00 0.960000 0.012566 0.057751
Data;
//0.5Hz
0.00 0.960000 0.062831853 0.057592
0.04 0.962137 0.062336405 0.062721
0.08 0.964590 0.060857875 0.067115
0.12 0.967043 0.05841958 0.069002
0.16 0.969327 0.055059973 0.070735
0.20 0.971487 0.050832037 0.072267
0.24 0.973477 0.04580245 0.067785
0.28 0.975240 0.04005053 0.064896
0.32 0.976757 0.03366699 0.065203
0.36 0.978013 0.026752502 0.061654
0.40 0.978987 0.01941611 0.060950
0.44 0.979673 0.011773515 0.055396
0.48 0.980030 0.003945245 0.052443
0.52 0.980090 -0.003945245 0.048458
0.56 0.979813 -0.011773515 0.042921
0.60 0.979240 -0.01941611 0.029484
0.64 0.978350 -0.026752502 0.019637
0.68 0.977177 -0.03366699 0.005958
0.72 0.975753 -0.04005053 -0.008437
0.76 0.974063 -0.04580245 -0.021190
0.80 0.972140 -0.050832037 -0.033505
0.84 0.970023 -0.055059973 -0.046834
0.88 0.967757 -0.05841958 -0.059828
0.92 0.965407 -0.060857875 -0.073131
0.96 0.962903 -0.062336405 -0.084756
1.00 0.960407 -0.062831853 -0.097015
1.04 0.957857 -0.062336405 -0.108150
1.08 0.955407 -0.060857875 -0.116848
1.12 0.952977 -0.05841958 -0.126935
1.16 0.950650 -0.055059973 -0.132818
1.20 0.948500 -0.050832037 -0.140228
1.24 0.946553 -0.04580245 -0.143198
1.28 0.944793 -0.04005053 -0.147031
1.32 0.943243 -0.03366699 -0.148441
1.36 0.941980 -0.026752502 -0.146859
1.40 0.941007 -0.01941611 -0.143694
1.44 0.940333 -0.011773515 -0.137398
1.48 0.939983 -0.003945245 -0.130185
1.52 0.939920 0.003945245 -0.118369
1.56 0.940163 0.011773515 -0.100830
1.60 0.940763 0.01941611 -0.078713
1.64 0.941643 0.026752502 -0.054448
1.68 0.942843 0.03366699 -0.031340
1.72 0.944277 0.04005053 -0.009725
1.76 0.945943 0.04580245 0.004953
1.80 0.947883 0.050832037 0.017859
1.84 0.949980 0.055059973 0.027375
1.88 0.952207 0.05841958 0.037944
1.92 0.954630 0.060857875 0.044558
1.96 0.957083 0.062336405 0.050863
2.00 0.960000 0.062831853 0.057592
Data;
//1Hz
0.00 0.959110 0.125664 0.101944
0.02 0.961607 0.124673 0.101015
0.04 0.964177 0.121716 0.092652
0.06 0.966597 0.116839 0.090688
0.08 0.968997 0.110120 0.087186
0.10 0.971203 0.101664 0.081292
0.12 0.973220 0.091605 0.074067
0.14 0.975063 0.080101 0.072352
0.16 0.976643 0.067334 0.064963
0.18 0.977983 0.053505 0.062168
0.20 0.979003 0.038832 0.059446
0.22 0.979740 0.023547 0.056947
0.24 0.980190 0.007890 0.053218
0.26 0.980303 -0.007890 0.051380
0.28 0.980107 -0.023547 0.049209
0.30 0.979583 -0.038832 0.050486
0.32 0.978767 -0.053505 0.045529
0.34 0.977637 -0.067334 0.042285
0.36 0.976253 -0.080101 0.040140
0.38 0.974577 -0.091605 0.040742
0.40 0.972660 -0.101664 0.028174
0.42 0.970597 -0.110120 0.015944
0.44 0.968337 -0.116839 -0.000255
0.46 0.965933 -0.121716 -0.015326
0.48 0.963463 -0.124673 -0.032873
0.50 0.960933 -0.125664 -0.051744
0.52 0.958353 -0.124673 -0.068495
0.54 0.955860 -0.121716 -0.084269
0.56 0.953367 -0.116839 -0.097367
0.58 0.951017 -0.110120 -0.107719
0.60 0.948813 -0.101664 -0.117079
0.62 0.946750 -0.091605 -0.120800
0.64 0.944930 -0.080101 -0.128407
0.66 0.943330 -0.067334 -0.127666
0.68 0.942017 -0.053505 -0.125147
0.70 0.940943 -0.038832 -0.120878
0.72 0.940243 -0.023547 -0.112602
0.74 0.939793 -0.007890 -0.102641
0.76 0.939683 0.007890 -0.090411
0.78 0.939893 0.023547 -0.066376
0.80 0.940393 0.038832 -0.039621
0.82 0.941243 0.053505 -0.011216
0.84 0.942383 0.067334 0.013463
0.86 0.943800 0.080101 0.032020
0.88 0.945480 0.091605 0.047094
0.90 0.947350 0.101664 0.060936
0.92 0.949430 0.110120 0.076479
0.94 0.951703 0.116839 0.086050
0.96 0.954087 0.121716 0.093122
0.98 0.956563 0.124673 0.102631
1.00 0.959110 0.125664 0.101944
Data;
//2Hz
0.00 0.957683 0.251327 0.103554
0.02 0.962873 0.243431 0.099815
0.04 0.967843 0.220240 0.089030
0.06 0.972333 0.183210 0.079861
0.08 0.976073 0.134668 0.066694
0.10 0.978863 0.077664 0.058112
0.12 0.980390 0.015781 0.053472
0.14 0.980647 -0.047094 0.048643
0.16 0.979643 -0.107010 0.043628
0.18 0.977403 -0.160202 0.036629
0.20 0.974053 -0.203328 0.016363
0.22 0.969813 -0.233678 -0.030502
0.24 0.964933 -0.249346 -0.082703
0.26 0.959757 -0.249346 -0.129067
0.28 0.954600 -0.233678 -0.165265
0.30 0.949830 -0.203328 -0.192737
0.32 0.945690 -0.160202 -0.204648
0.34 0.942393 -0.107010 -0.208398
0.36 0.940157 -0.047094 -0.197350
0.38 0.939237 0.015781 -0.169471
0.40 0.939640 0.077664 -0.110091
0.42 0.941337 0.134668 -0.036682
0.44 0.944197 0.183210 0.015458
0.46 0.948013 0.220240 0.052942
0.48 0.952623 0.243431 0.084528
0.50 0.957683 0.251327 0.103554
Data;
//5Hz
0.00 0.953477 0.628319 0.135643
0.02 0.966487 0.508320 0.098394
0.04 0.976950 0.194161 0.063313
0.06 0.981010 -0.194161 0.047463
0.08 0.977333 -0.508320 0.039043
0.10 0.966537 -0.628319 -0.039222
0.12 0.953520 -0.508320 -0.220273
0.14 0.943087 -0.194161 -0.271013
0.16 0.938930 0.194161 -0.210556
0.18 0.942740 0.508320 0.020174
0.20 0.953477 0.628319 0.135643 |
|