f(x)=input('請輸入函式f(x)=');
a=input('請輸入區間下限: ');
b=input('請輸入區間上限: ');
n=input('請輸入n（即將區間幾等分）: ');
h=(b-a)/n;
for i=1:n+1 %將區間n等分後，計算各節點的橫座標
xdata(i)=a+(i-1)*h;
endydata=subs(f,'x',xdata);
for k=1:n+1 %計算cotes係數
l=1;
for j=1:n+1
if(j~=k)
l=l*(t-j+1)/(k-j);
endend
s(k)=int(l,0,n)/n;
end
yy1=sum(s.*ydata);
yy2=int(f,a,b); %直接用matlab求積
disp('精確值');
vpa(yy2,8)
if=(b-a)*yy1; %用newton_cptes求積公式求得
disp('用newton_cptes求積公式求得：');
vpa(if,8)
wucha=yy2-if; %計算誤差
if (wucha<0)
wucha=-wucha;
end
disp('誤差為： ')
將區間2等分，即為梯形公式公式，n取1
顯然誤差相對較大
將區間3等分，即simpson公式，n取2
matlab實現kmeans演算法
kmeans是一種聚類演算法 無監督學習 演算法分為兩步 1.隨機選取k個聚類中心。2.計算每個樣本點離哪個聚類中心最近 距離計算 就將該樣本分為這個類。3.重新計算這k個類的聚類中心。一種簡單的計算方法為 計算每個類的平均值即為新的聚類中心。重複執行步驟2，直到聚類中心的變化小於給定閾值，或者達到...
Apriori演算法Matlab實現
clc clear 最小支援度設定 min sup 2 最小置信度 min conf 0.7 讀取檔案，當前的檔案型別是txt檔案，事務資料用數字來表示的，測試資料可以用 資料探勘概念與技術 第三版中的資料為樣本 fid fopen d matlabfile apriori dataapriori....
指派問題 MATLAB實現
適用於任意n階係數矩陣 clear all c 2 1513 4 10414 15 9 1416 13 7811 9 效率矩陣c n size c,1 計算c的行列數n c c 計算目標函式係數，將矩陣c按列排成乙個列向量即可。a b 沒有不等式約束 ae zeros 2 n,n 2 計算等約束的係...