matlab離散資料積分怎麼做?高手請進

2021-03-10 22:39:35 字數 1702 閱讀 9678

1樓:匿名使用者

t=0:0.01:2*pi;

x=sin(t);

s_x=cumtrapz(t,x);%s_x就是copyx的數值積分

plot(t,x,t,s_x)

___________

樓主請您明白一點,sin(t)的不定積分是-cos(t)+c,c是任意常數,而用數值方法作出來的是定積分,是sin(t)從0積到t的積分,因此在t=0時,就是0。根據牛頓-萊布尼茨定理,其實積出來的是-cos(t)+cos(0)即1-cos(t),這才是定積分的函式曲線。從你的貼圖上看,matlabguy的結果和我的結果也是一樣的,都是1-cos(t),而不是您所說的cos(t)+1。

你一定要得到-cos(t)的話,就把s_x再減去1即可。還有,-cos(t)和sin(t)的曲線不是重合曲線,是兩條不同的曲線,甚至相互正交。要處理離散資料積分,在這幾個回答之中也只有我這個方法可行了,因為你沒有函式關係,無法呼叫simpson函式。

2樓:匿名使用者

將下面內容儲存

bai為simpson.m

function y= simpson(f,a,b,m)% f被積函式;a積分下du

限;zhib積分上限;m子區間dao

個數回(將x分為多少個區間)

if m==0

y=0;

else

h=(b-a)/(2*m);

s1=0;

s2=0;

for i=1:m

x=a+(2*i-1)*h;

s1=s1+feval(f,x);

endfor j=1:(m-1)

x=a+2*j*h;

s2=s2+feval(f,x);

endy=h/3*(feval(f,a)+2*s2+4*s1+feval(f,b));

end在命令

答行輸入

f=@(t)sin(t);

i=;m=floor(2*pi/0.01);

for i=0:m

y=simpson(f,0,2*pi/m*i,i);

i=[i,y];

endt=0:0.01:2*pi;

plot(t,i);

我服了你了,sint和-cost曲線重合?sint積分得到的是-cost+c,常數c由初始條件決定,一般我們都認為從0開始積分,此時c=1。

資料積分的話用greatdju的方法就行了,他那個結果對,你自己琢磨吧。

3樓:匿名使用者

數值bai積分有很多種方法

du, 如下圖所示(來zhi自wikipeida)前面幾位的回答,分dao別對應了裡

回面一種演算法,另外答,數值積分得到的是一個值,如果你需要畫圖的話,那是代數積分。

儲存為mytest.m

function mytest

a=0;

b=2*pi;

n=100;

f=@(x)sin(x);

y=intsimpsons(f,a,b,n);

t=linspace(a,b,n);

y=intsimpsons(f,a,b,n)function i = intsimpsons(f,a,b,n)h=(b-a)/n; xi=a:h:b;

i= h/3*(f(xi(1))+2*sum(f(xi(3:2:end-2)))+4*sum(f(xi(2:2:end)))+f(xi(end)));

橋樑施工竣工怎麼做,橋樑施工竣工資料怎麼做?

按 市政基礎設施施工檔案管理規定 做 俗稱221檔案 上面有註名如何做竣工資料 根據 建築工程檔案管理基礎 規定,一般是這麼分的 a冊 綜合檔案 開工 復工 竣工報告 竣工驗收記錄 強制性條文,驗收表 設計變更通知 圖紙會審 施工組織設計等 b冊 質量保證資料 出廠合格證 檢驗報告 複試報告等 c冊...

市政園林怎麼做呢,市政園林資料怎麼做呢?

驗收資料主要包括 單位 分部 分項工程質量驗收評定表,竣工圖,工程量匯版總表 交工驗收權證書等等,不過各地區和不同系統的要求有所不同,比如不同省份的規範要求不 一 市政綠化和公路綠化 地產園林又有不同,具體的可以找監理單位要相對應的資料。如何做市政園林綠化工程資料 你好,市政園林綠化的資料並不難,比...

cosx 的n次方求積分怎麼做

答案如圖所示 一次bai積分完的話,那個du 函式zhi不是初等dao的 拓展資料 內 我只是來儲存答案的哈 cosx和sinx的n次方求積分的公式是什麼?0,2 cos x ndx 0,2 sin x ndx n 1 n n 3 n 2 4 5 2 3,n為奇數 n 1 n n 3 n 2 3 4...