P=（1e+5）*ones(n+1) 和 P=（1e+5）*eye(n+1) 是什么意思？谢谢啦

求极限lim n趋向于无穷（1/n）*n次方根下(n+1)(n+2)⋯(n+n)~

4/e。
记原式=P
P=[(n+1)(n+2)(n+3).(n+n)/n^n]^(1/n)
={[(n+1)/n][(n+2)/n][(n+3)/n].[(n+n)/n]}^(1/n)
=[(1+1/n)(1+2/n)(1+3/n).(1+n/n)]^(1/n)
取自然对数
lnP=(1/n)[ln(1+1/n)+ln(1+2/n)+ln(1+3/n)+.+ln(1+n/n)]
设f(x)=ln(1+x)
则P=[f(1/n)+f(2/n)+...+f(n/n)]/n
当n→∞时
应用分部积分法可求得

则当n→∞时,lnP=ln(4/e)，即P=4/e。
扩展资料：
极限的由来：

与一切科学的思想方法一样，极限思想也是社会实践的大脑抽象思维的产物。极限的思想可以追溯到古代，例如，祖国刘徽的割圆术就是建立在直观图形研究的基础上的一种原始的可靠的“不断靠近”的极限思想的应用；古希腊人的穷竭法也蕴含了极限思想，但由于希腊人“对’无限‘的恐惧”，他们避免明显地人为“取极限”，而是借助于间接证法——归谬法来完成了有关的证明。
到了16世纪，荷兰数学家斯泰文在考察三角形重心的过程中，改进了古希腊人的穷竭法，他借助几何直观，大胆地运用极限思想思考问题，放弃了归缪法的证明。如此，他就在无意中“指出了把极限方法发展成为一个实用概念的方向”。
极限的求法有很多种：
1、连续初等函数，在定义域范围内求极限，可以将该点直接代入得极限值，因为连续函数的极限值就等于在该点的函数值。
2、利用恒等变形消去零因子（针对于0/0型）。
3、利用无穷大与无穷小的关系求极限。
4、利用无穷小的性质求极限。
5、利用等价无穷小替换求极限，可以将原式化简计算。
6、利用两个极限存在准则，求极限，有的题目也可以考虑用放大缩小，再用夹逼定理的方法求极限。
7、利用两个重要极限公式求极限。

由于[(n+1)/e]^n<e(n/e)^n，因为(1+1/n)^n<e（当n趋无穷大时相等），可以转而证明不等式
e*[n/e]^n<n!<e*[(n+1)/e]^(n+1)；（顺带说明一下，这个不等式源自Stirling公式的推导）
取e为底的对数ln，转化成nln(n)-n+1<ln(n!)<(n+1)ln(n+1)-(n+1)+1，下面简单证明之：
因为对数函数单调增，于是int_{n-1}^n[ln x]dx <ln n <int_{n}^{n+1}[ln x]dx，其中int_a^b是积分域为(a,b)的积分；
又因为ln(n!)=ln1+ln2+...+ln n；
于是int_1^n[ln x]dx<ln(n!)<int_1^{n+1}[ln x]dx；（注意ln1=0，我扔掉了(0,1)上的lnx积分）
根据int_1^n[ln x]dx=nln x-n+1，欲证的不等式成立。
也许你需要证明一下(1+1/n)^n<e这个常用不等式吧，因为确实常见，只提示一下，详细的就不写了：
首先证明(1+1/n)^n是随n单调增的，你可以用取对数后泰勒展开作数列比较也可以直接多项式展开比较结合归纳法；其次证明当n趋无穷大时(1+1/n)^n=e，证明则可以取对数后泰勒展开，可证nln(1+1/n)=1。方法也许有许多。

P=（1e+5）*ones(n+1)表示10的5次方乘n+1阶的单位阵；P=（1e+5）*eye(n+1) 表示10的5次方乘n+1阶的对角阵，以下为程序
n=2;P=(1e+5)*ones(n+1)

P =

100000 100000 100000
100000 100000 100000
100000 100000 100000

>> n=2;P=(1e+5)*eye(n+1)

P =

100000 0 0
0 100000 0
0 0 100000

clear
%disp('单变量系统的隐式广义预测控制算法的研究')
%disp('广义预测控制算法初始值')
nn=6; %时域长度p
n=6; %预测长度n
m=2; %控制长度m
t0=0.8; %控制加权系数lamda
a=0.3; %柔化系数
%disp ('最小二乘公式初始值')
t1=1; %input('遗忘因子lamda1=')
d1=1; %(n+1)阶方阵P的形式：0对角阵，1方阵
d2=1; %(n+1)阶方阵P的初始值：0自动赋值1e+5，1键盘输入
if(d1==1)
if(d2==1)
P=6;
else
P=(1e+5)*ones(n+1);
end
else
if(d2==1)
P=6;
P=diag(P);
else
P=(1e+5)*eye(n+1);
end
end
%参数初始值
uuu=0;yyy=0;
uu=zeros(n,1);u=zeros(m,1);
yy=zeros(n,1);y1=zeros(n,1);
Q=zeros(n+1,1);Q(1,1)=1;Q(n+1,1)=1;
%产生周期为100，时间T，幅值为1的方波信号的给定值
T=300;[yr0,t]=gensig('step',100,T,1); %pulse sine
d3=0; %('输出曲线是否去掉前100步：0-不，1-去掉：');
nm=length(t); %确定循环次数
for ij=2:nm
yr=yr0(ij)+1; %产生周期为100，时间为T，幅值在1和2之间变化的方波信号的给定值
%根据系统模型，计算k时刻的输出值y（k）
y=1.49*yy(n,1)-0.49*yy(n-1,1)+0.5*uu(n-1,1);
%(nn=6,n=6,m=2;t0=0.8;a=0.3;t1=1) %example9-1
%y=2.001676*yy(n,1)-1.24339*yy(n-1,1)+0.24171*yy(n-2,1)+0.23589*uu(n,1); %example9-2
%产生均匀分布的白噪声
a9=0;
for i=1:1
a9=a9+rand;
end
a8=0.01*(a9-6);
%保存k时刻及以前的n个输出值y（k），y（k-1），...,y(k-n),以供模型运算
for i=1:n-1
yy(i,1)=yy(i+1,1);
end
yy(n,1)=y;
yyy=[yyy;y]; %保存个k时刻的nm个输出量以便绘图
%根据最小二乘公式，由y（k）计算G阵的各元素值g0,g1,...,gn
for i=1:n
X(1,i)=uu(i,1);
end
X(1,n+1)=1;
K=P*X'*inv(t1+X*P*X');
P=(eye(n+1)-K*X)*P/t1;
Q=Q+K*(y-X*Q);
%根据元素值g0,g1,...gn,求G阵
for j=1:m
for i=n:-1:j
i1=n-i+j;
G(i1,j)=Q(i,1);
end
end
%求nn维y0向量（y1为上一时刻的y0向量）
y0(1:nn-1,1)=y1(2:nn,1);y0(nn,1)=y1(nn,1);
y0=y0+(y-y1(1,1));
for i=1:n
y1(i,1)=y0(i,1)+u(1,1)*Q(n+1-i,1);
end
for i=n+1:nn
y1(i,1)=y1(n,1);
end
%根据y0,求n维f向量f(k+1),...f(k+n)
f(1:n,1)=y0(1:n,1);
%由当前k时刻的输出值y（k）和给定值yr，求k时刻以后的n个参考轨迹w（k+1），...w(k+n)
w=a*y+(1-a)*yr;
for i=2:n
w=[w;a^i*y+(1-a^i)*yr];
end
%计算k时刻及以后的m个控制增量Du（k），...Du（k+m）
u=inv(G'*G+t0*eye(m))*G'*(w-f);
%保存k时刻集以后的n个控制增量，以供模型运算
for i=1:n-1
uu(i,1)=uu(i+1,1);
end
uu(n,1)=u(1,1);
uuu=[uuu;u(1,1)]; %保存各k时刻的nm个控制增量以便绘图
%控制量限幅
if(u(1,1)>1)
u(1,1)=1;
end
if(u(1,1)<-1)
u(1,1)=-1;
end
end
%绘制给定值、输出值和控制增量曲线
if(d3==1)
%绘制去掉前100步的给定值、输出值和控制增量曲线
yyy1(1:(T-100),1)=yyy(101:T,1);
uuu1(1:(T-100),1)=yyy(101:T,1);
t1(1:(T-100),1)=t(101:T,1)-100;
yr01(1:(T-100),1)=yr0(101:T,1);
subplot(2,1,1):plot(t1,(yr01+1),t1,yyy1,'-');
axis([0,nm-100,0,2.5]);
xlabel('t'); ylabel('yr,y');
subplot(2,1,2):plot(t1,uuu1,'*');
axis([0,nm-100,-1.5,1.5]);
legend('yr','-';'y','*');
xlabel('t'); ylabel('u')

else
% 绘制完整的给定值、输出值和控制增量曲线
subplot(2,1,1):plot(t,(yr0+1),t,yyy);
axis([0,nm,0,2.5]);
xlabel('t'); ylabel('yr,y'); legend('yr,y')
subplot(2,1,2):plot(t,uuu);
axis([0,nm,-1.5,1.5]);
xlabel('t'); ylabel('u')
end

??? Output argument "u" (and maybe others) not assigned during call to "D:\Program Files\MATLAB71\toolbox\control\control\gensig.m (gensig)".

Error in ==> gensig at 21
ni = nargin;

---

杂多县19380226759： 用c语言编程计算π/2 - ？
夹须复方： #include<stdio.h> int main() { double i=1,y=1,t=1; for(;t>=1e-6;i++) {t*=i/(2*i+1);y+=t; } printf(＂Pi/2=%lf\nPi=%lf\n＂,y,2*y); getch(); }

杂多县19380226759： 如何将k - means算法与de算法结合用matlab实现 - ？
夹须复方： data=input('请输入样本数据矩阵:'); m=size(data,1); n=size(data,2); counter=0; k=input('请输入聚类数目:'); whilek>m disp('您输入的聚类数目过大,请输入正确的 k 值'); k=input('请输入聚类数目:'); end if k==1 disp('聚类数...

杂多县19380226759： MATLAB题目:用π/4=1 - 1/3+1/5 - 1/7+......公式求π的近似值,直到某一项的绝对值小于10^( - 6)为止. - ？
夹须复方： sum=0 n=1 while 1/(2*n-1)>1e-6 sum=sum+1/(2*n-1)*((-1)^(n-1)); n=n+1; end n sum*4 例如:do while (abs(j-i)>=1/1000) j=1 i=j+1/(2i+1) j=i enddo print(j) 扩展资料: 计算机语言中,正数的二进制符号位)为0,负数的二进制首位为1. 32位系...

杂多县19380226759： 一个数列前n项和为Sn=1 - 2+3 - 4----+( - 1)的n+1次方*n,则S17+S33+S50=? - ？
夹须复方： n为偶数,Sn=(1-2)+(3-4)+……+(n-1)*(-1)^n+n*(-1)^(n+1)=-n/2 n为奇数,Sn=1+(-2+3)+(-4+5)+……+(n-1)*(-1)^n+n*(-1)^(n+1)=(n+1)/2 S17=(17+1)/2=9 S33=(33+1)/2=17 S50=-50/2=-25 S17+S33+S50=1

杂多县19380226759： int a[]={1,2,3,4,5,6,7},*p=a int n,s=0; for(n=0;n<6;n++) s+=p[n++] printf(＂%d＂,s) - ？
夹须复方： s的值为1+3+5(即a[0]+a[2]+a[4]),所以为9

杂多县19380226759： 在Fibonacci问题中,若a(0)=1, a(1)=1, a(n+1)=a(n)+a(n - 1), 则a(n...？
夹须复方： #include<stdio.h> main() { int n,s; n=1; s=0; while(n<21) { s=s+n*(n+1)*(n+2); n++; } printf(＂S=%d＂,s);}

杂多县19380226759： lim趋向无穷,5n的3次方分之(n+1)x(n+2)x(n+3)为什么等于5分之一? - ？
夹须复方： 记原式=P, P=[(n+1)(n+2)(n+3).(n+n)/n^n]^(1/n) ={[(n+1)/n][(n+2)/n][(n+3)/n].[(n+n)/n]}^(1/n) =[(1+1/n)(1+2/n)(1+3/n).(1+n/n)]^(1/n) 取自然对数, lnP=(1/n)[ln(1+1/n)+ln(1+2/n)+ln(1+3/n)+.+ln(1+n/n)] 设f(x)=ln(1+x), 则P=[f(1/n)+f(2/n)+...+f(n/n)]/n, 当n→∞时, 应用分部积分法可求得 则当n→∞时,lnP=ln(4/e),即P=4/e.

杂多县19380226759： 各项均为正数的数列{an}中,a1=a,a2=b,且满足m+n=p+q的正整数m,n,p,q都有am+an/(1+am)(1+an)= - ？
夹须复方： (am+an)/【(1+am)(1+an)】=(ap+aq)/【(1+ap)(1+aq)】 所以(a1+an)/(1+a1)(1+an)=(a2+a[n-1])/(1+a2)(1+a[n-1]) 带入a=0.5,b=0.8 an=(2a[n-1]+1)/(an+2) 所以数列{1+an /1-an}是一个等比数列 an=(3^n -1)/(3^n +1)==================...

杂多县19380226759： 已知数列cn=2n^2+3n+1,求证:当n>=2时,数列{1/cn}的前n项和小于5/12 - ？
夹须复方：[答案] 1/cn=(2n+1)(n+1) 数列{1/cn}的前n项和=1/(2*3)+1/(3*5)+1/(4*7)+.+1/(2n+1)(n+1) 1/(3*5)+1/(4*7)+.+1/(2n+1)(n+1)