chap9.bas

介绍一些matlab在数学方面的应用

Attribute VB_Name = "chap9"
Public Sub PutOnClipboard(str As Variant)
    Dim MyDataObj As New DataObject
    MyDataObj.SetText Format(str)
    MyDataObj.PutInClipboard
End Sub
Sub i1()
CR = Chr(13) + Chr(10)
PutOnClipboard ("x=[0:15]'; y1=[]; y2=[]; lam1=[1,2,5,10];" + CR + "for i=1:length(lam1)" + CR + "y1=[y1,poisspdf(x,lam1(i))]; y2=[y2,poisscdf(x,lam1(i))];" + CR + "end" + CR + "plot(x,y1), figure; plot(x,y2)")
End Sub
Sub i2()
CR = Chr(13) + Chr(10)
PutOnClipboard ("x=[-5:.02:5]'; y1=[]; y2=[];" + CR + "mu1=[-1,0,0,0,1]; sig1=[1,0.1,1,10,1]; sig1=sqrt(sig1);" + CR + "for i=1:length(mu1)" + CR + "y1=[y1,normpdf(x,mu1(i),sig1(i))]; y2=[y2,normcdf(x,mu1(i),sig1(i))];" + CR + "end" + CR + "plot(x,y1), figure; plot(x,y2)")
End Sub
Sub i3()
CR = Chr(13) + Chr(10)
PutOnClipboard ("x=[-0.5:.02:5]'; y1=[]; y2=[];" + CR + "a1=[1,1,2,1,3]; lam1=[1,0.5,1,2,1];" + CR + "for i=1:length(a1)" + CR + "y1=[y1,gampdf(x,a1(i),lam1(i))]; y2=[y2,gamcdf(x,a1(i),lam1(i))];" + CR + "end" + CR + "plot(x,y1), figure; plot(x,y2)")
End Sub
Sub i4()
PutOnClipboard ("x=[-eps:-0.02:-0.05,0:0.02:5]; x=sort(x');")
End Sub
Sub i5()
CR = Chr(13) + Chr(10)
PutOnClipboard ("x=[-eps:-0.02:-0.05,0:0.02:2]; x=sort(x');" + CR + "k1=[1,2,3,4,5]; y1=[]; y2=[];" + CR + "for i=1:length(k1)" + CR + "y1=[y1,chi2pdf(x,k1(i))]; y2=[y2,chi2cdf(x,k1(i))];" + CR + "end" + CR + "plot(x,y1), figure; plot(x,y2)")
End Sub
Sub i6()
CR = Chr(13) + Chr(10)
PutOnClipboard ("x=[-5:0.02:5]'; k1=[1,2,5,10]; y1=[]; y2=[];" + CR + "for i=1:length(k1)" + CR + "y1=[y1,tpdf(x,k1(i))]; y2=[y2,tcdf(x,k1(i))];" + CR + "end" + CR + "plot(x,y1), figure; plot(x,y2)")
End Sub
Sub i7()
CR = Chr(13) + Chr(10)
PutOnClipboard ("x=[-eps:-0.02:-0.05,0:0.02:5]; x=sort(x');" + CR + "b1=[.5,1,3,5]; y1=[]; y2=[];" + CR + "for i=1:length(b1)" + CR + "y1=[y1,raylpdf(x,b1(i))]; y2=[y2,raylcdf(x,b1(i))];" + CR + "end" + CR + "plot(x,y1),[ figure; plot(x,y2)")
End Sub
Sub i8()
CR = Chr(13) + Chr(10)
PutOnClipboard ("x=[-eps:-0.02:-0.05,0:0.02:1]; x=sort(x');" + CR + "p1=[1 2 3 3 4]; q1=[1 1 1 2 1]; y1=[]; y2=[];" + CR + "for i=1:length(p1)" + CR + "y1=[y1,fpdf(x,p1(i),q1(i))]; y2=[y2,fcdf(x,p1(i),q1(i))];" + CR + "end" + CR + "plot(x,y1), figure; plot(x,y2)")
End Sub
Sub i9()
PutOnClipboard ("b=1; p1=raylcdf(0.2,b); p2=raylcdf(2,b); P1=p2-p1")
End Sub
Sub i10()
PutOnClipboard ("p1=raylcdf(1,b); P2=1-p1")
End Sub
Sub i11()
PutOnClipboard ("syms x y; f=x^2+x*y/3; P=int(int(f,x,0,1/2),y,0,1/2)")
End Sub
Sub i12()
CR = Chr(13) + Chr(10)
PutOnClipboard ("b=1; p=raylrnd(1,30000,1);" + CR + "xx=0:.1:4; yy=hist(p,xx); yy=yy/(30000*0.1);" + CR + "bar(xx,yy), y=raylpdf(xx,1); line(xx,y)")
End Sub
Sub i13()
CR = Chr(13) + Chr(10)
PutOnClipboard ("syms x; syms a lam positive" + CR + "p=lam^a*x^(a-1)/gamma(a)*exp(-lam*x);" + CR + " m=int(x*p,x,0,inf)")
End Sub
Sub i14()
PutOnClipboard ("s=simple(int((x-1/lam*a)^2*p,x,0,inf))")
End Sub
Sub i15()
PutOnClipboard ("p=normrnd(0.5,1.5,30000,1); [mean(p), var(p), std(p)]")
End Sub
Sub i16()
CR = Chr(13) + Chr(10)
PutOnClipboard ("p=normrnd(0.5,1.5,300,1); [mean(p), var(p), std(p)]")
End Sub
Sub i17()
PutOnClipboard ("[m,s]=raylstat(0.45)")
End Sub
Sub i18()
CR = Chr(13) + Chr(10)
PutOnClipboard ("syms x; syms a lam positive;" + CR + " p=lam^a*x^(a-1)/gamma(a)*exp(-lam*x);" + CR + "for n=1:5, m=int(x^n*p,x,0,inf), end")
End Sub
Sub i19()
PutOnClipboard ("syms n; m=simple(int((x)^n*p,x,0,inf))")
End Sub
Sub i20()
PutOnClipboard ("for n=1:7, s=simple(int((x-1/lam*a)^n*p,x,0,inf)), end")
End Sub
Sub i21()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[]; B=[]; p=normrnd(0.5,1.5,30000,1); n=1:5;" + CR + "for r=n," + CR + "A=[A, sum(p.^r)/length(p)]; B=[B,moment(p,r)];" + CR + "end" + CR + "A,B")
End Sub
Sub i22()
CR = Chr(13) + Chr(10)
PutOnClipboard ("syms x; A1=[]; B1=[];" + CR + " p=1/(sqrt(2*pi)*1.5)*exp(-(x-0.5)^2/(2*1.5^2));" + CR + "for i=1:5" + CR + "A1=[A1,vpa(int(x^i*p,x,-inf,inf),12)];" + CR + "B1=[B1,vpa(int((x-0.5)^i*p,x,-inf,inf),12)];" + CR + "end")
End Sub
Sub i23()
PutOnClipboard ("A1, B1")
End Sub
Sub i24()
PutOnClipboard ("p=randn(30000,4); cov(p)")
End Sub
Sub i25()
CR = Chr(13) + Chr(10)
PutOnClipboard ("mu1=[-1,2]; Sigma2=[1 1; 1 3]; % 输入均值向量和协方差矩阵" + CR + "[X,Y]=meshgrid(-3:0.1:1,-2:0.1:4); xy=[X(:) Y(:)]; % 产生网格数据并处理" + CR + "p=mvnpdf(xy,mu1,Sigma2); P=reshape(p,size(X)); % 求取联合概率密度" + CR + "surf(X,Y,P) % 绘制联合概率密度的三维网格图")
End Sub
Sub i26()
CR = Chr(13) + Chr(10)
PutOnClipboard ("Sigma2=diag(diag(Sigma2)); % 消除协方差矩阵的非对角元素" + CR + "p=mvnpdf(xy,mu1,Sigma2);" + CR + " P=reshape(p,size(X)); surf(X,Y,P)")
End Sub
Sub i27()
CR = Chr(13) + Chr(10)
PutOnClipboard ("mu1=[-1,2]; Sigma2=[1 1; 1 3];" + CR + "R1=mvnrnd(mu1,Sigma2,2000); plot(R1(:,1),R1(:,2),'o')" + CR + "Sigma2=diag(diag(Sigma2)); figure;" + CR + "R2=mvnrnd(mu1,Sigma2,2000); plot(R2(:,1),R2(:,2),'o')")
End Sub
Sub i28()
CR = Chr(13) + Chr(10)
PutOnClipboard ("p=gamrnd(1.5,3,30000,1);" + CR + " Pv=[0.9,0.92,0.95,0.98]; A=[];" + CR + "for i=1:length(Pv)" + CR + "[a,b]=gamfit(p,Pv(i)); A=[A; Pv(i),a(1),b(:,1)',a(2),b(:,2)'];" + CR + "end")
End Sub
Sub i29()
CR = Chr(13) + Chr(10)
PutOnClipboard ("num=[300,3000,30000,300000,3000000]; A=[];" + CR + "for i=1:length(num)" + CR + "p=gamrnd(1.5,3,num(i),1);" + CR + "[a,b]=gamfit(p,0.95); A=[A;num(i),a(1),b(:,1)',a(2),b(:,2)'];" + CR + "end")
End Sub
Sub i30()
CR = Chr(13) + Chr(10)
PutOnClipboard ("a=[1 -1.232 2.23 2 4,3.792]'; X=randn(120,6); " + CR + "y=X*a; a1=inv(X'*X)*X'*y")
End Sub
Sub i31()
PutOnClipboard ("[a,aint]=regress(y,X,0.02)")
End Sub
Sub i32()
PutOnClipboard ("yhat=y+sqrt(0.5)*randn(120,1); [a,aint]=regress(yhat,X,0.02)")
End Sub
Sub i33()
PutOnClipboard ("errorbar(1:6,a,aint(:,1)-a,aint(:,2)-a)")
End Sub
Sub i34()
CR = Chr(13) + Chr(10)
PutOnClipboard ("yhat=y+sqrt(0.1)*randn(120,1); [a,aint]=regress(yhat,X,0.02);" + CR + "errorbar(1:6,a,aint(:,1)-a,aint(:,2)-a)")
End Sub
Sub i35()
CR = Chr(13) + Chr(10)
PutOnClipboard ("f=inline('a(1)*exp(-a(2)*x)+a(3)*exp(-a(4)*x).*sin(a(5)*x)','a','x');" + CR + "x=0:0.1:10; y=f([0.12,0.213,0.54,0.17,1.23],x);" + CR + "[a,r,j]=nlinfit(x,y,f,[1;1;1;1;1]); a")
End Sub
Sub i36()
PutOnClipboard ("ci=nlparci(a,r,j)")
End Sub
Sub i37()
CR = Chr(13) + Chr(10)
PutOnClipboard ("y=f([0.12,0.213,0.54,0.17,1.23],x)+0.02*rand(size(x));" + CR + "[a,r,j]=nlinfit(x,y,f,[1;1;1;1;1]); a")
End Sub
Sub i38()
PutOnClipboard ("ci=nlparci(a,r,j)")
End Sub
Sub i39()
PutOnClipboard ("errorbar(1:5,a,ci(:,1)-a,ci(:,2)-a)")
End Sub
Sub i40()
CR = Chr(13) + Chr(10)
PutOnClipboard ("a=[1;1;1;1;1;1]';" + CR + "f=inline(['(a(1)*x(:,1).^3+a(2)).*sin(a(3)*x(:,2).*x(:,3))+',..." + CR + "'(a(4)*x(:,3).^3+a(5)*x(:,3)+a(6))'],'a','x');" + CR + "X=randn(120,4); y=f(a,X)+sqrt(0.2)*randn(120,1);")
End Sub
Sub i41()
PutOnClipboard ("[ahat,r,j]=nlinfit(X,y,f,[0;2;3;2;1;2]); ahat")
End Sub
Sub i42()
PutOnClipboard ("y1=f(ahat,X); plot([y y1])")
End Sub
Sub i43()
PutOnClipboard ("ci=nlparci(ahat,r,j); ci")
End Sub
Sub i44()
PutOnClipboard ("errorbar(1:6,ahat,ci(:,1)-ahat,ci(:,2)-ahat)")
End Sub
Sub i45()
CR = Chr(13) + Chr(10)
PutOnClipboard ("n=200; mu0=9.94; xbar=9.73; s=1.62;" + CR + " u=sqrt(n)*(mu0-xbar)/s")
End Sub
Sub i46()
CR = Chr(13) + Chr(10)
PutOnClipboard ("alpha=[0.01:0.01:0.05 0.07, 0.09];" + CR + " K=norminv(1-alpha/2,0,1); [alpha' K']")
End Sub
Sub i47()
PutOnClipboard ("abs(u)<K")
End Sub
Sub i48()
PutOnClipboard ("r=normrnd(1,2,400,1); [H,p,ci]=ztest(r,1,2,0.02)")
End Sub
Sub i49()
PutOnClipboard ("[H,p,ci]=ztest(r,0.5,2,0.02)")
End Sub
Sub i50()
PutOnClipboard ("[H,p,ci]=ttest(r,1,0.02)")
End Sub
Sub i51()
CR = Chr(13) + Chr(10)
PutOnClipboard ("X=[216,203,197,208,206,209,206,208,202,203,206,213,218,207,208,..." + CR + "202,194,203,213,211,193,213,208,208,204,206,204,206,208,209,..." + CR + "213,203,206,207,196,201,208,207,213,208,210,208,211,211,214,...")
End Sub
Sub i52()
CR = Chr(13) + Chr(10)
PutOnClipboard ("220,211,203,216,224,211,209,218,214,219,211,208,221,211,218,..." + CR + "218,190,219,211,208,199,214,207,207,214,206,217,214,201,212,..." + CR + "213,211,212,216,206,210,216,204,221,208,209,214,214,199,204,...")
End Sub
Sub i53()
CR = Chr(13) + Chr(10)
PutOnClipboard ("211,201,216,211,209,208,209,202,211,207,202,205,206,216,206,..." + CR + "213,206,207,200,198,200,202,203,208,216,206,222,213,209,219];" + CR + "[H,p]=jbtest(X,0.05)")
End Sub
Sub i54()
CR = Chr(13) + Chr(10)
PutOnClipboard ("[mu1,sig1,mu_ci,sig_ci]=normfit(X,0.05);" + CR + "mu=[mu1,mu_ci'],sig=[sig1, sig_ci']")
End Sub
Sub i55()
PutOnClipboard ("r=gamrnd(1,3,400,1); [H,p,c,d]=jbtest(r,0.05)")
End Sub
Sub i56()
PutOnClipboard ("r=gamrnd(1,3,400,1); alam=gamfit(r)")
End Sub
Sub i57()
PutOnClipboard ("[H0,p]=kstest(r,gamcdf(sort(r),alam(1),alam(2)),0.05)")
End Sub
Sub i58()
CR = Chr(13) + Chr(10)
PutOnClipboard ("A=[5,4,6,7,9; 8,6,4,4,3; 7,6,4,6,5; 7,3,5,6,7; 10,5,4,3,7; 8,6,3,5,6];" + CR + "mean(A)")
End Sub
Sub i59()
PutOnClipboard ("[p,tbl,stats]=anova1(A)")
End Sub
Sub i60()
CR = Chr(13) + Chr(10)
PutOnClipboard ("B=[23,15,26,13,21,25,20,21,16,18,21,17,16,24,27,14,17,19,20,24;" + CR + "28,22,25,19,26,30,26,26,20,28,19,24,19,25,29,17,21,18,26,23;" + CR + "18,10,12,22,13,15,21,22,14,12,23,25,19,13,22,16,12,23,22,19];" + CR + "anova2(B',5);")
End Sub
Sub i61()
CR = Chr(13) + Chr(10)
PutOnClipboard ("C=[];" + CR + "for i=1:3" + CR + "for j=1:4" + CR + "C(i,j)=mean(B(i,[1:5]+(j-1)*5));" + CR + "end, end" + CR + "C=[C; mean(C)]; C=[C mean(C')']")
End Sub
Sub i62()
CR = Chr(13) + Chr(10)
PutOnClipboard ("")
End Sub