Matlab求解包含阶跃函数的数值积分 求MATLAB数值积分
\u8bf7\u6559\u5404\u8def\u5927\u795e\uff0cmatlab\u7ea7\u6570\u6c42\u548c\u540e\uff0c\u6c42\u89e3\u6570\u503c\u79ef\u5206\u3002\u65b9\u6cd5/\u6b65\u9aa4
\u6700\u5e38\u89c1\u7684\u7ea7\u6570\u5f62\u5f0f\uff0c\u5982\u4e0b\u6240\u793a\u3002
S \uff1a\u7ea7\u6570\u7684\u548c
i\uff1a \u81ea\u53d8\u91cf\uff0c\u503c\u57df\u4e3a[a,b]
f(i)\uff1a\u4e3a\u5173\u4e8e\u81ea\u53d8\u91cfi\u7684\u51fd\u6570
\u67e5\u9605matlab\u5e2e\u52a9\u53ef\u4ee5\u77e5\u9053\uff0c\u53ef\u4ee5\u7528symsum\u51fd\u6570\u6c42\u89e3\u7ea7\u6570\uff0c\u51fd\u6570\u6709\u4ee5\u4e0b4\u79cd\u4f7f\u7528\u65b9\u5f0f\uff1a
\u4f5c\u4e3a\u793a\u4f8b\uff0c\u5c0f\u7f16\u91c7\u7528\u7b2c\u56db\u79cd\u65b9\u5f0f\u6c42\u89e3\u4ee5\u4e0b\u7ea7\u6570\uff1a
clc;clear;
% symsum(expr,v,a,b)
syms v
syms f
a = 1;
b = 100;
f = v^2;
S = symsum(f,v,a,b)
\u70b9\u51fb\u56fe\u4e2d\u7ea2\u8272\u7bad\u5934\u6240\u6307\u7684\u8fd0\u884c\u6309\u94ae\uff0c\u8ba1\u7b97\u7ed3\u679c\u4e3a\uff1a
S = 338350
\u901a\u8fc7int()\u51fd\u6570\u548cvpa()\u51fd\u6570\u6c42\u5f97K\u503c\uff08\u6570\u503c\u89e3\uff09\uff0c\u7136\u540e\u7528polt()\u51fd\u6570\u753b\u51fab\u3001K\u7684\u5173\u7cfb\u56fe\u3002
\u4ee3\u7801\u89c1\u9644\u4ef6\u3002
syms x
F=heaviside(x);
int(F,x,-10,10)
结果:
ans =
10
绛旓細L(x'')+0.2L(x')-0.4L(x)=0.2L(u(t))(p^2)*X(p)-p*x(0)-x'(0)+0.2(p*X(p)-x(0))-0.4*X(p)=0.2/p,(p^2+0.2*p-0.4)*X(p)=0.2/p,X(p)=0.2/(p^3+0.2*(p^2)-0.4p)鐒跺悗锛屼綘鏌aplace Transform 琛ㄥ嵆鍙傛垨鑰鐢∕Atlab 姹傝В,浠g爜濡備笅锛欶=0....
绛旓細姝ょ▼搴忕敤浜庣敓鎴愰樁璺冧俊鍙 t 鎬绘椂闀 s 姝ラ暱 t0 闃惰穬寮濮嬫椂鍒 a0 闃惰穬鍓嶅箙鍊 a1 闃惰穬鍚庣殑骞呭 j=0;for i=0:s:t j=j+1;if i < t0 y(j)=a0;else y(j)=a1;end end 鍐插嚮淇″彿浣犲彲浠ュ娉曠偖鍒躲傛柟娉曚笁锛氫綘涔熷彲浠ヨ繖鏍峰仛锛Matlab涓鍙疯绠椾腑鎻愪緵 鍗曚綅闃惰穬鍑芥暟heaviside(t-a)鍗曚綅鑴夊啿鍑芥暟dirac...
绛旓細impulse鍑芥暟鏄啿鍑诲嚱鏁帮紱step鍑芥暟鏄闃惰穬鍑芥暟
绛旓細xx=fft(x);yy=fft(y);hh=yy./xx;h=ifft(hh);濡傛灉瀵规偍鏈夊府鍔╋紝璇疯寰楅噰绾充负婊℃剰绛旀锛岃阿璋紒绁濇偍鐢熸椿鎰夊揩锛乿aela
绛旓細n1=0:10;Xn1=0.8*exp(n1).*(heaviside(n1)-heaviside(n1-5));n2=0:4;Xn2=ones(1,5);nys=n1(1)+n2(1);nyf=n1(end)+n2(end);y=conv(Xn1,Xn2);ny=[nys:nyf];figure(1)stem(n1,Xn1);figure(2);stem(n2,Xn2);figure(3)stem(y,ny);...
绛旓細浠庢媺姘忓彉鎹㈢殑瀹氫箟 鍙煡閬擄紝鍙樻崲鐨勭Н鍒嗚寖鍥存槸浠0鍒版鏃犵┓ 鎴栬呰鏄粯璁や簡x(t)褰搕<0鏃朵负0锛屽張鎴栬呮媺姘忓彉鎹㈠氨瑕佹眰鍘鍑芥暟鐨瀹氫箟鍩熸槸0鍒版鏃犵┓ u(t-a)褰揳澶т簬闆剁殑鏃跺欙紝鏄痷(t)鍚戝彸绉诲姩锛屾弧瓒硉<0鏃秛(t)=0 鑰屽綋a灏忎簬闆剁殑鏃跺欙紝鏄痷(t)鍚戝乏绉诲姩锛屾湁涓閮ㄥ垎t<0鏃秛(t)涓嶄负0 濡傛灉榛樿t<0...
绛旓細matlab涓湁涓撻棬鐨闃惰穬鍑芥暟銆俿tep鍑芥暟銆
绛旓細t=0:0.1:10;%璁惧畾t鐨勮寖鍥 f=3-exp(t);plot(t,f)
绛旓細num=[0,0,0,20000];den=[1,205,1000,20000];t=0:0.05:2.5;sys=tf(num,den)Transfer function:20000 --- s^3 + 205 s^2 + 1000 s + 20000 鎴戠殑娌℃湁闂锛屼綘杈撳叆help tf 鐪嬬湅鏈夋病鏈夎繖涓鍑芥暟鐨鐩稿叧淇℃伅鍑虹幇銆傛病鏈夌殑璇濓紝璇存槑浣犵殑matlab娌℃湁杩欎釜宸ュ叿绠憋紝瑕佸畨瑁呰繖涓伐鍏风 ...
绛旓細闃惰穬鍑芥暟 [y0,t0]=step(G,0:5:500)浜х敓涓涓嚎鎬фā鍨婫鐨勯樁璺冨搷搴旓紝鏃堕棿浠0鍒500锛屾闀挎槸5锛岃繑鍥炶緭鍑哄搷搴攜0,鏃堕棿鐭㈤噺t0