matlab求解数值积分,知F(x,y),求Fx(x)=∫F(x,y)dy,用什么样的调用格式,x保留符号,积分限已知 matlab数值积分定义被积函数的句柄函数定义遇到问题:被积...
MATLAB\u4e2d\u6c42Y=F(X)\u6570\u503c\u79ef\u5206quadl(@(x) f(x),x1,x2),\u5982 y=sin(x)\u5728\uff081,5\uff09\u4e0a\u7684\u79ef\u5206\uff0cquadl(@(x) sin(x),1,5)
quad2d\u4e0d\u80fd\u4f5c\u7528\u4e8e\u77e9\u9635\u548c\u5411\u91cf\uff1bquadv\u53ef\u4ee5\uff0c\u4e0d\u8fc7\u90a3\u662f\u5bf91\u91cd\u79ef\u5206\u7684\uff1b\u767e\u5ea6\u77e5\u9053\u9650\u5236100\u5b57\uff0c\u6211\u5206\u5f00\u5199\u4e0b\u4ee3\u7801
clear,clc
syms x y
Q = [x,x+y,sin(x);x*y,2*x*y,y];
\u5269\u4f59\u89c1\u8bc4\u8bba
Maple: int(x^(1.3333^y),y,0,1)
Matlab: int((x^(1.3333^y),y,0,1):计算f关于x在区间[a,b]上的定积分。(出来的是数值,必须先定义x,否则就出错了)
因此要就算你那个积分用matlab直接积分是不行的,当然也可以用符号表达式进行
>>syms x
等等……
当时matlab调用符号表达式的时候其实都是maple的库,因此又回到了原点。
int(f,x,a,b):计算f关于x在区间[a,b]上的定积分。
对于Fx(X)=∫x^(1.333^y)dy,积分限为上1,下0,可令f=x^(1.333^y),a=0,b=1。
对于积分限是正负无穷,可令a,或者b为inf。
我才小学
绛旓細Matlab: int((x^(1.3333^y),y,0,1)锛氳绠f鍏充簬x鍦ㄥ尯闂碵a,b]涓婄殑瀹氱Н鍒嗐(鍑烘潵鐨勬槸鏁板硷紝蹇呴』鍏堝畾涔墄锛屽惁鍒欏氨鍑洪敊浜)鍥犳瑕佸氨绠椾綘閭d釜绉垎鐢╩atlab鐩存帴绉垎鏄笉琛岀殑锛屽綋鐒朵篃鍙互鐢ㄧ鍙疯〃杈惧紡杩涜 >>syms x 绛夌瓑鈥︹﹀綋鏃秏atlab璋冪敤绗﹀彿琛ㄨ揪寮忕殑鏃跺欏叾瀹為兘鏄痬aple鐨勫簱锛屽洜姝ゅ張鍥炲埌浜嗗師鐐广
绛旓細鍏朵腑锛宍function_name`鏄绉垎鐨勫嚱鏁帮紝`a`鍜宍b`鏄Н鍒嗙殑涓婁笅闄愩備緥濡傦紝瑕佽绠楀嚱鏁f=x^2鍦ㄥ尯闂碵0, 1]鐨绉垎锛鍙互浣跨敤濡備笅浠g爜锛matlab result = integral x.^2, 0, 1);杩欓噷鐨刞x.^2`琛ㄧずx鐨勫钩鏂广備娇鐢ㄥ尶鍚嶅嚱鏁癅灏嗚緭鍏鏄犲皠鍒板畠鐨勫钩鏂逛笂銆傚洜姝わ紝integral浼氳绠椾粠x=0鍒皒=1鐨勭Н鍒嗗笺2....
绛旓細1銆佷娇鐢╥nt鍑芥暟锛屽嚱鏁扮敱integrate缂╁啓鑰屾潵锛宨nt 鍑芥暟琛ㄨ揪寮忥紝鍙橀噺锛岀Н鍒嗕笂闄愶紝绉垎涓嬮檺銆2銆佹瘮濡姹涓涓狥x = a*x^2,鍦ㄥ尯闂达紙m锛宯锛夊x杩涜绉垎锛棣栧厛瑕佸皢 m,x,a,b 杩欏洓涓彉閲忓畾涔変负绗﹀彿鍙橀噺 syms m x a b;Fx = a*x^2;int(Fx,x,m,n)3銆侀氳繃涓婇潰杩欎釜鏂规硶锛屽氨鑳藉姹傚緱浠绘剰涓涓嚱鏁板湪缁...
绛旓細X = 0:pi/100:pi;Y = sin(X);Z = pi/100*trapz(Y)2銆佽嚜閫傚簲杈涙櫘妫鏁板肩Н鍒璁$畻 quad锛堬級F = @(x)1./(x.^3-2*x-5);Q = quad(F,0,2);3銆佽嚜閫傚簲Lobatto绉垎璁$畻 quadl锛堬級function y = myfun(x)y = 1./(x.^3-2*x-5);end Q = quadl(@myfun,0,2);4銆佽嚜閫傚簲Gauss...
绛旓細quadl(@(x) f(x),x1,x2),濡 y=sin(x)鍦紙1,5锛変笂鐨绉垎锛quadl(@(x) sin(x),1,5)
绛旓細杈撳叆鎸囦护锛歝d G:\matlab\work锛屾槸璇ヨ矾寰勬垚涓哄綋鍓嶈矾寰勩4銆乵atlab涓绠楃Н鍒嗙殑涓ょ鎸囦护锛1.F=quad('fname',a,b,tol,trace) Simpson鏁板肩Н鍒娉 2.F=quad8('fname',a,b,tol,trace) Newton-Cotes鏁板肩Н鍒嗘硶 鍏朵腑锛 fname鏄绉嚱鏁拌〃杈惧紡鎴栧嚱鏁板悕锛宎,b鍒嗗埆鏄笂涓嬮檺锛宼ol鍙互鎺у埗绉垎绮惧害...
绛旓細x=double(solve(f2))鏍规嵁浣犳彁渚涘浘鐗囩畻寮忥紝鍙戠幇浣犵紪鍐欑殑绠楀紡浠g爜鏈夌偣闂銆傛纭柟寮忎负 B=10^6;Tf=2*10^(-3);theta=10^(-3);beta=theta*Tf*B/log(2);syms y x eq0=(1/(x^(1/(beta+1)))*(1/(y^(beta/(beta+1)))*exp(-y/10);f=int(eq0,y,x,+inf)-10;x=double(solv...
绛旓細1銆佸皢瑕佷娇鐢MATLAB璁$畻涓嬪浘涓殑浜岄噸绉垎锛棣栧厛鍦ㄧ數鑴戜笂鎵撳紑MATLAB杞欢銆2銆佹柊寤鸿剼鏈紙Ctrl+N锛夛紝杈撳叆鍥句腑妗嗕綇鐨勪唬鐮佸唴瀹广傚叾涓璔1=dblquad(f,0,2*pi,-pi,pi,1.0e-3)閲囩敤榛樿鏂规硶quad璁$畻浜岄噸绉垎锛岀粷瀵硅绠楃簿搴﹁涓1.0e-3銆俻lot3(x,y,ff,'r','LineWidth',3)鏄粯鍒惰绉垎鍑芥暟ff=x.*sin(y...
绛旓細浜庢棤娉曟眰寰梕xp(x^2)鐨勫師鍑芥暟锛屾垜浠彧鑳界敤鏁板肩畻娉曟潵姹傝В锛鍙互鐢ㄥ鍖栨褰㈠叕寮忋丷omberg鍏紡銆丟auss鍏紡绛夛紝鏈夊ソ澶氱銆傛垜鐢∕atlab缂栦簡涓涓敤Gauss鍏紡姹傝В绉垎鐨勫嚱鏁般俧unction S=GaussIntegrate()杩愮敤Gauss姹傜Н鍏紡璁$畻鏁板肩Н鍒 f涓鸿绉嚱鏁帮紝Rho涓烘潈鍑芥暟锛屼簩鑰呭潎涓虹鍙峰嚱鏁 x=sym('x');f=exp(x^2)...
绛旓細1銆乹uad鍑芥暟鐢ㄤ簬璁$畻鏁板肩Н鍒嗭紝鍑芥暟琛ㄨ揪寮忎腑涓嶈兘鍖呭惈绗﹀彿閲忥紱2銆佽绉嚱鏁扮殑琛ㄨ揪寮忓簲璇ュ啓鎴愬叧浜庤绉彉閲忕殑鍚戦噺鍖栫殑褰㈠紡锛堜篃灏辨槸搴旇鐢ㄧ偣杩愮畻锛夈傚弬鑰冧唬鐮侊細R=1;syms L;rr = 0 : 0.1 : 1;for ii = 1 : length(rr) r = rr(ii); f = @(l)(acos((1+l*l-r*r)/(2*l))+r*...