24小时热门版块排行榜    

查看: 2210  |  回复: 8

daijiantao

新虫 (正式写手)

[求助] 用Matlab化简一个复杂含参积分式的函数表达式已有1人参与

急需化简一个复杂含参积分的表达式,表达式如上传图片所示:
用Matlab化简一个复杂含参积分式的函数表达式

由于积分式比较复杂,之前用matlab做了数值解,并验证正确了
format long
clc,clear
M=0.230873094;
N=0.025914445;
n=0.86;
exp=114.8;
wf0=46;
tH=300;
F = @(t,p1,p2,p3,p4,p5)(((p1./(1+t)-p2.*log10(1+t)./(1+t)).^p3./p4)-(p1./(1+t)-p2.*log10(1+t)./(1+t))./p5);
dc=quadl(@(x)F(x,M,N,n,exp,wf0),0,tH)

由于数值解分别计算一个一个案例非常麻烦,所以想要寻求一下图片所示公式的化简公式,并可以用于Excel等软件中批量处理。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

“数值解分别计算一个一个案例”,对于不同的案例,积分式什么参数会改变?
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2016-04-07 13:51:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daijiantao

新虫 (正式写手)

引用回帖:
2楼: Originally posted by 月只蓝 at 2016-04-07 13:51:21
“数值解分别计算一个一个案例”,对于不同的案例,积分式什么参数会改变?

M N exp n wf0 tH会变。所以想找一个含参的近似方程,怎样用MATLAB实现呢~
3楼2016-04-07 23:45:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyyacgh

铁杆木虫 (职业作家)

彪悍的人生,不需要签名!
4楼2016-04-08 07:21:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
3楼: Originally posted by daijiantao at 2016-04-07 23:45:21
M N exp n wf0 tH会变。所以想找一个含参的近似方程,怎样用MATLAB实现呢~...

M N exp n wf0 tH会变,写个循环语句就解决了。
原积分式即便近似为一个代数式,M N exp n wf0 tH会变的情况下,一样麻烦。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
5楼2016-04-08 09:30:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daijiantao

新虫 (正式写手)

引用回帖:
5楼: Originally posted by 月只蓝 at 2016-04-08 09:30:09
M N exp n wf0 tH会变,写个循环语句就解决了。
原积分式即便近似为一个代数式,M N exp n wf0 tH会变的情况下,一样麻烦。...

如果我只有积分上限变化(tH),其余都不变,能基于上面的程序帮忙给一个程序吗,最好可以出图?
6楼2016-04-10 23:15:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^

syms p1 p2 p3 p4 p5 tH t
F = (((p1./(1+t)-p2.*log10(1+t)./(1+t)).^p3./p4)-(p1./(1+t)-p2.*log10(1+t)./(1+t))./p5);
ff=int(F,t,0,tH);
latex(ff)

结果(不咋的):
ans =

\int _{0}^{{\it tH}}\! \left( {\frac {{\it p1}}{1+t}}-{\frac {{\it p2}\,\ln  \left( 1+t \right) }{\ln  \left( 10 \right)  \left( 1+t \right) }} \right) ^{{\it p3}}{{\it p4}}^{-1}- \left( {\frac {{\it p1}}{1+t}}-{\frac {{\it p2}\,\ln  \left( 1+t \right) }{\ln  \left( 10 \right)  \left( 1+t \right) }} \right) {{\it p5}}^{-1}{dt}
abcd
7楼2016-04-10 23:56:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
daijiantao: 金币+20, ★★★★★最佳答案, 非常感谢! 2016-04-11 21:34:58
引用回帖:
6楼: Originally posted by daijiantao at 2016-04-10 23:15:42
如果我只有积分上限变化(tH),其余都不变,能基于上面的程序帮忙给一个程序吗,最好可以出图?...

以下代码是积分上限tH从100开始以间隔1变化到600的例子,可自行根据实际修改。
需注意以下代码为m文件格式,请新建一个m文件,把以下代码完整复制进去,再运行。
CODE:
function qiu_jifen
clear all;clc
format long
tH_0=100;    %积分上限tH的最小值
tH_n=500;    %积分上限tH的最大值
num=tH_n-tH_0+1;     %积分上限tH总的变化数目

tH_data=linspace(tH_0,tH_n,num);
for i=1:length(tH_data)
    tH=tH_data(i);
    dc(i)=quadgk(@obj,0,tH);
end
[tH_data' dc']
figure(1)
plot(tH_data,dc,'b-')

function f=obj(t)
M=0.230873094;
N=0.025914445;
n=0.86;
exp=114.8;
wf0=46;
p1=M;
p2=N;
p3=n;
p4=exp;
p5=wf0;
f=(((p1./(1+t)-p2.*log10(1+t)./(1+t)).^p3./p4)-(p1./(1+t)-p2.*log10(1+t)./(1+t))./p5);

计算结果:
积分上限tH          数值积分结果
ans =

  1.0e+002 *

   1.000000000000000  -0.000062764008049
   1.010000000000000  -0.000062774137557
   1.020000000000000  -0.000062783623116
   1.030000000000000  -0.000062792482273
   1.040000000000000  -0.000062800731978
   1.050000000000000  -0.000062808388605
   1.060000000000000  -0.000062815467978
   1.070000000000000  -0.000062821985394
   1.080000000000000  -0.000062827955645
   1.090000000000000  -0.000062833393040
   1.100000000000000  -0.000062838311419
   1.110000000000000  -0.000062842724180
   1.120000000000000  -0.000062846644291
   1.130000000000000  -0.000062850084306
   1.140000000000000  -0.000062853056388
   1.150000000000000  -0.000062855572318
   1.160000000000000  -0.000062857643510
   1.170000000000000  -0.000062859281030
   1.180000000000000  -0.000062860495604
   1.190000000000000  -0.000062861297633
   1.200000000000000  -0.000062861697206
   1.210000000000000  -0.000062861704109
   1.220000000000000  -0.000062861327838
   1.230000000000000  -0.000062860577611
   1.240000000000000  -0.000062859462376
   1.250000000000000  -0.000062857990818
   1.260000000000000  -0.000062856171376
   1.270000000000000  -0.000062854012244
   1.280000000000000  -0.000062851521385
   1.290000000000000  -0.000062848706536
   1.300000000000000  -0.000062845575217
   1.310000000000000  -0.000062842134738
   1.320000000000000  -0.000062838392207
   1.330000000000000  -0.000062834354535
   1.340000000000000  -0.000062830028448
   1.350000000000000  -0.000062825420484
   1.360000000000000  -0.000062820537008
   1.370000000000000  -0.000062815384213
   1.380000000000000  -0.000062809968127
   1.390000000000000  -0.000062804294619
   1.400000000000000  -0.000062798369402
   1.410000000000000  -0.000062792198041
   1.420000000000000  -0.000062785785954
   1.430000000000000  -0.000062779138422
   1.440000000000000  -0.000062772260588
   1.450000000000000  -0.000062765157464
   1.460000000000000  -0.000062757833934
   1.470000000000000  -0.000062750294760
   1.480000000000000  -0.000062742544584
   1.490000000000000  -0.000062734587930
   1.500000000000000  -0.000062726429213
   1.510000000000000  -0.000062718072737
   1.520000000000000  -0.000062709522699
   1.530000000000000  -0.000062700783196
   1.540000000000000  -0.000062691858225
   1.550000000000000  -0.000062682751684
   1.560000000000000  -0.000062673467380
   1.570000000000000  -0.000062664009028
   1.580000000000000  -0.000062654380255
   1.590000000000000  -0.000062644584600
   1.600000000000000  -0.000062634625523
   1.610000000000000  -0.000062624506398
   1.620000000000000  -0.000062614230524
   1.630000000000000  -0.000062603801121
   1.640000000000000  -0.000062593221338
   1.650000000000000  -0.000062582494248
   1.660000000000000  -0.000062571622857
   1.670000000000000  -0.000062560610099
   1.680000000000000  -0.000062549458846
   1.690000000000000  -0.000062538171902
   1.700000000000000  -0.000062526752009
   1.710000000000000  -0.000062515201849
   1.720000000000000  -0.000062503524043
   1.730000000000000  -0.000062491721154
   1.740000000000000  -0.000062479795691
   1.750000000000000  -0.000062467750106
   1.760000000000000  -0.000062455586798
   1.770000000000000  -0.000062443308113
   1.780000000000000  -0.000062430916348
   1.790000000000000  -0.000062418413750
   1.800000000000000  -0.000062405802518
   1.810000000000000  -0.000062393084803
   1.820000000000000  -0.000062380262713
   1.830000000000000  -0.000062367338308
   1.840000000000000  -0.000062354313608
   1.850000000000000  -0.000062341190589
   1.860000000000000  -0.000062327971185
   1.870000000000000  -0.000062314657291
   1.880000000000000  -0.000062301250763
   1.890000000000000  -0.000062287753418
   1.900000000000000  -0.000062274167037
   1.910000000000000  -0.000062260493362
   1.920000000000000  -0.000062246734103
   1.930000000000000  -0.000062232890932
   1.940000000000000  -0.000062218965490
   1.950000000000000  -0.000062204959383
   1.960000000000000  -0.000062190874186
   1.970000000000000  -0.000062176711442
   1.980000000000000  -0.000062162472664
   1.990000000000000  -0.000062148159335
   2.000000000000000  -0.000062133772907
   2.010000000000000  -0.000062119314805
   2.020000000000000  -0.000062104786427
   2.030000000000000  -0.000062090189141
   2.040000000000000  -0.000062075524292
   2.050000000000000  -0.000062060793195
   2.060000000000000  -0.000062045997143
   2.070000000000000  -0.000062031137402
   2.080000000000000  -0.000062016215215
   2.090000000000000  -0.000062001231800
   2.100000000000000  -0.000061986188353
   2.110000000000000  -0.000061971086047
   2.120000000000000  -0.000061955926033
   2.130000000000000  -0.000061940709439
   2.140000000000000  -0.000061925437374
   2.150000000000000  -0.000061910110925
   2.160000000000000  -0.000061894731159
   2.170000000000000  -0.000061879299122
   2.180000000000000  -0.000061863815843
   2.190000000000000  -0.000061848282330
   2.200000000000000  -0.000061832699574
   2.210000000000000  -0.000061817068546
   2.220000000000000  -0.000061801390200
   2.230000000000000  -0.000061785665474
   2.240000000000000  -0.000061769895287
   2.250000000000000  -0.000061754080542
   2.260000000000000  -0.000061738222127
   2.270000000000000  -0.000061722320913
   2.280000000000000  -0.000061706377753
   2.290000000000000  -0.000061690393489
   2.300000000000000  -0.000061674368946
   2.310000000000000  -0.000061658304934
   2.320000000000000  -0.000061642202248
   2.330000000000000  -0.000061626061671
   2.340000000000000  -0.000061609883971
   2.350000000000000  -0.000061593669902
   2.360000000000000  -0.000061577420206
   2.370000000000000  -0.000061561135611
   2.380000000000000  -0.000061544816832
   2.390000000000000  -0.000061528464574
   2.400000000000000  -0.000061512079528
   2.410000000000000  -0.000061495662371
   2.420000000000000  -0.000061479213772
   2.430000000000000  -0.000061462734387
   2.440000000000000  -0.000061446224861
   2.450000000000000  -0.000061429685826
   2.460000000000000  -0.000061413117907
   2.470000000000000  -0.000061396521714
   2.480000000000000  -0.000061379897851
   2.490000000000000  -0.000061363246909
   2.500000000000000  -0.000061346569469
   2.510000000000000  -0.000061329866103
   2.520000000000000  -0.000061313137374
   2.530000000000000  -0.000061296383834
   2.540000000000000  -0.000061279606028
   2.550000000000000  -0.000061262804489
   2.560000000000000  -0.000061245979743
   2.570000000000000  -0.000061229132307
   2.580000000000000  -0.000061212262690
   2.590000000000000  -0.000061195371392
   2.600000000000000  -0.000061178458903
   2.610000000000000  -0.000061161525708
   2.620000000000000  -0.000061144572283
   2.630000000000000  -0.000061127599095
   2.640000000000000  -0.000061110606605
   2.650000000000000  -0.000061093595265
   2.660000000000000  -0.000061076565521
   2.670000000000000  -0.000061059517811
   2.680000000000000  -0.000061042452565
   2.690000000000000  -0.000061025370208
   2.700000000000000  -0.000061008271158
   2.710000000000000  -0.000060991155824
   2.720000000000000  -0.000060974024611
   2.730000000000000  -0.000060956877916
   2.740000000000000  -0.000060939716130
   2.750000000000000  -0.000060922539637
   2.760000000000000  -0.000060905348818
   2.770000000000000  -0.000060888144043
   2.780000000000000  -0.000060870925681
   2.790000000000000  -0.000060853694091
   2.800000000000000  -0.000060836449630
   2.810000000000000  -0.000060819192646
   2.820000000000000  -0.000060801923485
   2.830000000000000  -0.000060784642484
   2.840000000000000  -0.000060767349976
   2.850000000000000  -0.000060750046291
   2.860000000000000  -0.000060732731750
   2.870000000000000  -0.000060715406673
   2.880000000000000  -0.000060698071372
   2.890000000000000  -0.000060680726154
   2.900000000000000  -0.000060663371324
   2.910000000000000  -0.000060646007180
   2.920000000000000  -0.000060628634016
   2.930000000000000  -0.000060611252121
   2.940000000000000  -0.000060593861780
   2.950000000000000  -0.000060576463274
   2.960000000000000  -0.000060559056880
   2.970000000000000  -0.000060541642868
   2.980000000000000  -0.000060524221507
   2.990000000000000  -0.000060506793061
   3.000000000000000  -0.000060489357789
   3.010000000000000  -0.000060471915947
   3.020000000000000  -0.000060454467787
   3.030000000000000  -0.000060437013556
   3.040000000000000  -0.000060419553500
   3.050000000000000  -0.000060402087857
   3.060000000000000  -0.000060384616865
   3.070000000000000  -0.000060367140758
   3.080000000000000  -0.000060349659765
   3.090000000000000  -0.000060332174111
   3.100000000000000  -0.000060314684021
   3.110000000000000  -0.000060297189713
   3.120000000000000  -0.000060279691403
   3.130000000000000  -0.000060262189304
   3.140000000000000  -0.000060244683626
   3.150000000000000  -0.000060227174576
   3.160000000000000  -0.000060209662356
   3.170000000000000  -0.000060192147167
   3.180000000000000  -0.000060174629206
   3.190000000000000  -0.000060157108668
   3.200000000000000  -0.000060139585745
   3.210000000000000  -0.000060122060624
   3.220000000000000  -0.000060104533491
   3.230000000000000  -0.000060087004530
   3.240000000000000  -0.000060069473922
   3.250000000000000  -0.000060051941842
   3.260000000000000  -0.000060034408467
   3.270000000000000  -0.000060016873969
   3.280000000000000  -0.000059999338517
   3.290000000000000  -0.000059981802279
   3.300000000000000  -0.000059964265420
   3.310000000000000  -0.000059946728102
   3.320000000000000  -0.000059929190484
   3.330000000000000  -0.000059911652725
   3.340000000000000  -0.000059894114980
   3.350000000000000  -0.000059876577402
   3.360000000000000  -0.000059859040141
   3.370000000000000  -0.000059841503346
   3.380000000000000  -0.000059823967175
   3.390000000000000  -0.000059806431749
   3.400000000000000  -0.000059788897221
   3.410000000000000  -0.000059771363731
   3.420000000000000  -0.000059753831418
   3.430000000000000  -0.000059736300416
   3.440000000000000  -0.000059718770860
   3.450000000000000  -0.000059701242881
   3.460000000000000  -0.000059683716610
   3.470000000000000  -0.000059666192174
   3.480000000000000  -0.000059648669685
   3.490000000000000  -0.000059631149296
   3.500000000000000  -0.000059613631116
   3.510000000000000  -0.000059596115264
   3.520000000000000  -0.000059578601860
   3.530000000000000  -0.000059561091021
   3.540000000000000  -0.000059543582862
   3.550000000000000  -0.000059526077497
   3.560000000000000  -0.000059508575038
   3.570000000000000  -0.000059491075595
   3.580000000000000  -0.000059473579278
   3.590000000000000  -0.000059456086193
   3.600000000000000  -0.000059438596446
   3.610000000000000  -0.000059421110141
   3.620000000000000  -0.000059403627380
   3.630000000000000  -0.000059386148265
   3.640000000000000  -0.000059368672895
   3.650000000000000  -0.000059351201369
   3.660000000000000  -0.000059333733782
   3.670000000000000  -0.000059316270230
   3.680000000000000  -0.000059298810807
   3.690000000000000  -0.000059281355605
   3.700000000000000  -0.000059263904715
   3.710000000000000  -0.000059246458228
   3.720000000000000  -0.000059229016231
   3.730000000000000  -0.000059211578812
   3.740000000000000  -0.000059194146056
   3.750000000000000  -0.000059176718049
   3.760000000000000  -0.000059159294872
   3.770000000000000  -0.000059141876609
   3.780000000000000  -0.000059124463341
   3.790000000000000  -0.000059107055146
   3.800000000000000  -0.000059089652104
   3.810000000000000  -0.000059072254292
   3.820000000000000  -0.000059054861786
   3.830000000000000  -0.000059037474662
   3.840000000000000  -0.000059020092992
   3.850000000000000  -0.000059002716850
   3.860000000000000  -0.000058985346308
   3.870000000000000  -0.000058967981437
   3.880000000000000  -0.000058950622305
   3.890000000000000  -0.000058933268983
   3.900000000000000  -0.000058915921537
   3.910000000000000  -0.000058898580034
   3.920000000000000  -0.000058881244539
   3.930000000000000  -0.000058863915118
   3.940000000000000  -0.000058846591834
   3.950000000000000  -0.000058829274749
   3.960000000000000  -0.000058811963925
   3.970000000000000  -0.000058794659424
   3.980000000000000  -0.000058777361305
   3.990000000000000  -0.000058760069627
   4.000000000000000  -0.000058742784448
   4.010000000000000  -0.000058725505825
   4.020000000000000  -0.000058708233816
   4.030000000000000  -0.000058690968475
   4.040000000000000  -0.000058673709857
   4.050000000000000  -0.000058656458015
   4.060000000000000  -0.000058639213004
   4.070000000000000  -0.000058621974875
   4.080000000000000  -0.000058604743680
   4.090000000000000  -0.000058587519469
   4.100000000000000  -0.000058570302292
   4.110000000000000  -0.000058553092198
   4.120000000000000  -0.000058535889236
   4.130000000000000  -0.000058518693453
   4.140000000000000  -0.000058501504897
   4.150000000000000  -0.000058484323612
   4.160000000000000  -0.000058467149645
   4.170000000000000  -0.000058449983040
   4.180000000000000  -0.000058432823842
   4.190000000000000  -0.000058415672093
   4.200000000000000  -0.000058398527837
   4.210000000000000  -0.000058381391115
   4.220000000000000  -0.000058364261969
   4.230000000000000  -0.000058347140439
   4.240000000000000  -0.000058330026565
   4.250000000000000  -0.000058312920387
   4.260000000000000  -0.000058295821944
   4.270000000000000  -0.000058278731274
   4.280000000000000  -0.000058261648413
   4.290000000000000  -0.000058244573400
   4.300000000000000  -0.000058227506271
   4.310000000000000  -0.000058210447061
   4.320000000000000  -0.000058193395805
   4.330000000000000  -0.000058176352538
   4.340000000000000  -0.000058159317295
   4.350000000000000  -0.000058142290108
   4.360000000000000  -0.000058125271011
   4.370000000000000  -0.000058108260036
   4.380000000000000  -0.000058091257215
   4.390000000000000  -0.000058074262579
   4.400000000000000  -0.000058057276159
   4.410000000000000  -0.000058040297986
   4.420000000000000  -0.000058023328089
   4.430000000000000  -0.000058006366497
   4.440000000000000  -0.000057989413240
   4.450000000000000  -0.000057972468345
   4.460000000000000  -0.000057955531841
   4.470000000000000  -0.000057938603755
   4.480000000000000  -0.000057921684114
   4.490000000000000  -0.000057904772944
   4.500000000000000  -0.000057887870271
   4.510000000000000  -0.000057870976121
   4.520000000000000  -0.000057854090519
   4.530000000000000  -0.000057837213489
   4.540000000000000  -0.000057820345056
   4.550000000000000  -0.000057803485244
   4.560000000000000  -0.000057786634075
   4.570000000000000  -0.000057769791574
   4.580000000000000  -0.000057752957761
   4.590000000000000  -0.000057736132661
   4.600000000000000  -0.000057719316293
   4.610000000000000  -0.000057702508680
   4.620000000000000  -0.000057685709843
   4.630000000000000  -0.000057668919801
   4.640000000000000  -0.000057652138576
   4.650000000000000  -0.000057635366187
   4.660000000000000  -0.000057618602654
   4.670000000000000  -0.000057601847995
   4.680000000000000  -0.000057585102230
   4.690000000000000  -0.000057568365376
   4.700000000000000  -0.000057551637451
   4.710000000000000  -0.000057534918474
   4.720000000000000  -0.000057518208462
   4.730000000000000  -0.000057501507431
   4.740000000000000  -0.000057484815398
   4.750000000000000  -0.000057468132379
   4.760000000000000  -0.000057451458391
   4.770000000000000  -0.000057434793449
   4.780000000000000  -0.000057418137569
   4.790000000000000  -0.000057401490764
   4.800000000000000  -0.000057384853051
   4.810000000000000  -0.000057368224444
   4.820000000000000  -0.000057351604955
   4.830000000000000  -0.000057334994601
   4.840000000000000  -0.000057318393393
   4.850000000000000  -0.000057301801345
   4.860000000000000  -0.000057285218470
   4.870000000000000  -0.000057268644781
   4.880000000000000  -0.000057252080290
   4.890000000000000  -0.000057235525009
   4.900000000000000  -0.000057218978950
   4.910000000000000  -0.000057202442124
   4.920000000000000  -0.000057185914543
   4.930000000000000  -0.000057169396217
   4.940000000000000  -0.000057152887158
   4.950000000000000  -0.000057136387376
   4.960000000000000  -0.000057119896881
   4.970000000000000  -0.000057103415682
   4.980000000000000  -0.000057086943791
   4.990000000000000  -0.000057070481216
   5.000000000000000  -0.000057054027967
用Matlab化简一个复杂含参积分式的函数表达式-1
附图1.png

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
8楼2016-04-11 10:00:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daijiantao

新虫 (正式写手)

引用回帖:
8楼: Originally posted by 月只蓝 at 2016-04-11 10:00:43
以下代码是积分上限tH从100开始以间隔1变化到600的例子,可自行根据实际修改。
需注意以下代码为m文件格式,请新建一个m文件,把以下代码完整复制进去,再运行。

function qiu_jifen
clear all;clc
format l ...

非常感谢
9楼2016-04-11 21:59:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 daijiantao 的主题更新
信息提示
请填处理意见