设直径为1的圆为⊙O, 分别做⊙O的外切正n边型, 及内接正n边型, 分别记它们周长 为An和Bn, 则圆的周长S满足 An < S < Bn, 其An,Bn 满足如下迭代关系式:A2n = 2Bn*An/(Bn + An) , B2n =sqrt( Bn*A2n )帮我用matlab将上述的迭代点(Bn, An), n = 3, 6, . . . , 3 × 2^k, ...分别画出两个折线图, 并用不同的点型, 线性, 颜色表示.
时间: 2023-03-27 09:00:37 浏览: 68
好的,我可以回答这个问题。以下是用MATLAB绘制迭代点的代码:
n = 3;
k = 1;
Bn = zeros(1, 2^k);
An = zeros(1, 2^k);
for i = 1:2^k
Bn(i) = sqrt(2)/2;
An(i) = n*Bn(i);
n = n*2;
end
figure;
plot(Bn, '-o', 'LineWidth', 2, 'MarkerSize', 10, 'MarkerFaceColor', 'r');
hold on;
plot(An, '-s', 'LineWidth', 2, 'MarkerSize', 10, 'MarkerFaceColor', 'b');
xlabel('n');
ylabel('迭代点');
legend('Bn', 'An');
title('外切正n边型和内接正n边型的迭代点');
其中,'-o'表示用圆圈标记点,'-s'表示用正方形标记点,'LineWidth'表示线条宽度,'MarkerSize'表示标记大小,'MarkerFaceColor'表示标记填充颜色。您可以根据需要自行调整这些参数。
相关问题
设直径为1的圆为⊙O, 分别做⊙O的外切正n边型, 及内接正n边型, 分别记它们周长 为Pn和pn, 则圆的周长S满足 pn < S < Pn, 其中, pn, Pn满足如下迭代关系式, P2n = 2pn*Pn/ (pn + Pn) , p2n = sqrt( pn*P2n ) 帮我用matlab将上述的迭代点(pn, Pn), n = 3, 6, . . . , 3 × 2^k, ...分别画出两个折线图, 并用不同的点 型, 线性, 颜色表示.
我可以回答这个问题。以下是用 Matlab 画出迭代点的代码:
n = 3;
k = ;
while n <= 3 * 2^k
Pn = 1;
pn = cos(pi/n);
for i = 1:10
P2n = 2 * pn * Pn / (pn + Pn);
pn = sqrt(pn * P2n);
Pn = P2n;
end
plot(n, pn, 'ro');
hold on;
plot(n, Pn, 'bx');
n = n * 2;
k = k + 1;
end
xlabel('n');
ylabel('Pn, pn');
legend('pn', 'Pn');
grid on;
其中,cos(pi/n) 是正 n 边型的半径,10 是迭代次数,可以根据需要调整。折线图中,红色圆圈表示 pn,蓝色叉号表示 Pn。可以根据需要修改点型、线性和颜色。
设直径为1的圆为⊙O, 分别做⊙O的外切正n边型, 及内接正n边型, 分别记它们周长 为Pn和pn, 则圆的周长S满足pn < S < Pn,其中, pn, Pn满足如下迭代关系式, P2n = 2pn*Pn /(pn + Pn ), p2n = sqrt(pn*P2n) 帮我用matlab将上述的迭代点(pn, Pn), n = 3, 6, . . . , 3 × 2 ^k 分别画出两个折线图, 并用不同的点型, 线性, 颜色表示.
我可以回答这个问题。以下是用MATLAB绘制迭代点的代码:
n = 3:3:3*2^k; % 定义n的取值范围
Pn = zeros(1,length(n)); % 初始化Pn和pn
pn = zeros(1,length(n));
Pn(1) = 2; % 初始值
pn(1) = sqrt(3);
for i = 2:length(n)
Pn(i) = 2*pn(i-1)*Pn(i-1)/(pn(i-1)+Pn(i-1)); % 计算Pn和pn
pn(i) = sqrt(pn(i-1)*Pn(i));
end
S = pi; % 圆的周长
for i = 1:length(n)
S = S + Pn(i); % 计算S
end
figure; % 绘制折线图
plot(n,Pn,'-o','LineWidth',2,'MarkerSize',8,'MarkerFaceColor','r'); % 绘制Pn的折线图
hold on;
plot(n,pn,'-s','LineWidth',2,'MarkerSize',8,'MarkerFaceColor','b'); % 绘制pn的折线图
legend('P_n','p_n');
xlabel('n');
ylabel('Value');
title('Iteration Points');
grid on;
阅读全文