1樓:匿名使用者
使用圖中給出的遞推公式計算:
n = 8;
s = zeros(n+1,1);
s(1) = log(6)-log(5);
for i = 2 : n+1,
s(i) = 1/i - 5 * s(i-1);
ends(1) = ;
num2str(s,'%.3f')
得到的結果為
ans =
-0.412
2.391
-11.707
58.734
-293.505
1467.667
-7338.211
36691.166
但如果直接採用數值積分的方法計算:
s = zeros(n,1);
for i = 1 : n
s(i) = quadl(@(x)x.^i./(x+5),0,1);
ends
得到結果為
s =0.0884
0.0580
0.0431
0.0343
0.0285
0.0243
0.0212
0.0188
第二種方法得到的結果是正確的,第一種方法雖然遞推公式是精確成立的,但該公式卻是數值不穩定的,不能用於實際的計算。
2樓:塗智華
編寫一函式calcu如下
function s=calcu(n)
s=0.182;
for i=1:n
s=1/i-5*s;
ends
用matlab求π的近似值 π/4=1-1/3+1/5-1/7……直到某一項的絕對值小於10^-6
3樓:匿名使用者
syms n
x=(-1)^(n-1)/(2*n-1);
n=1;
xx=eval(x);%%%可換為xx=subs(x)y=0;
while abs(xx)>10^(-6)y=y+xx;
n=n+1;
xx=eval(x);
enddisp(['π的近似值為' num2str(y*4)])上面的方法由於涉及符號運算,所以特別特別慢,我試了一下15分鐘都沒算出來,但程式肯定沒問題,因為當精度降低後還是有結果的,下面這種方法很快,因為這是matlab所擅長的數值計算。
y=0;
for n=1:1000000%%%這裡要盡量的大xx=(-1)^(n-1)/(2*n-1);
if abs(xx)>10^(-6)
y=y+xx;
else
break
endendfprintf('π的近似值為%.6f ,共需迭代%d步\n',y*4,n)
matlab求π的近似值π/2=(1+1/(1*3))*(1+1/(3*5))*...
4樓:
half_pi = 1;
a = 1; b = 3;
for k = 1:1e4
half_pi = half_pi*(1+1/(a*b));
a = b;
b = b+2;
endhalf_pi
求教各位大神matlab如何使用牛頓迭代法求下面方程的近似值?
5樓:匿名使用者
牛頓迭代法的原理是,
實現過程如下:
function y = fun_main(x) %定義函式y =x^3+4*x^2-10;
endfunction y = dfun_main(x) %定義導函式
y =x^3+4*x^2-10;
y = 3*x^2 + 8*x;
endsyms x
ezplot('y=x^3+4*x^2-10',[-5 4 -11 5]),grid on,hold on
title('y=x^3+4x^2-10','fontsize',12)
ylabel('y','fontsize',12)xlabel('x','fontsize',12)x0=1.5;tol=1e-6; %x0初值,tol精度x=newton(x0,tol) %牛頓迭代函式y=fun_main(x)
plot(x,y,'rp')
執行結果及圖形
matlab作業急求大神,matlab作業急求大神
r 0.2 syms t phi 15 pi t 位移 x r cos phi r y r sin phi 速度來自 dx diff x dy diff y v sqrt dx dx dy dy v v 加速bai度 a diff v t1 0 4 1500 2 15 x1 subs x,t,t1 ...
求白酒業大神指導
瓶裝酒必須有生產許可證,現在很多地方限制白酒發展,不給頒發新的生產許可證了。我是公司業務員,做白酒的,老闆要求寫乙份 與客戶溝通的技巧 求大神幫忙 1 做事要沉穩,語氣要適中。表現出產品安全 負責任的。2 公司實力雄厚 技術力量投入大,產品受到專家認可。3 產品質量與 是配套的,可以與同等 的產品比...
求大神解答指導,速度速度
2a b b 2a 3 a b 3a 4a 3b 3a b b 2a 4a b 3 4a 4ab b 12a 9ab 6a 5ab b 4a b 12a 12ab 3b 12a 9ab 6a 5ab b 10a 8ab 9b 10 1 8 1 2 9 2 10 1 8 2 9 4 10 16 36 ...