求大神指導用matlab求近似值的問題(跪求)

2021-08-11 11:21:24 字數 2191 閱讀 2336

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