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...