个人收集的matlab代码截止到0528.docx
- 文档编号:13133398
- 上传时间:2023-06-11
- 格式:DOCX
- 页数:113
- 大小:1.43MB
个人收集的matlab代码截止到0528.docx
《个人收集的matlab代码截止到0528.docx》由会员分享,可在线阅读,更多相关《个人收集的matlab代码截止到0528.docx(113页珍藏版)》请在冰点文库上搜索。
个人收集的matlab代码截止到0528
Rbf神经网络预测:
clc
closeall
clear
%ld=100;%thenumberoflearningdata
%x=rand(2,200);%generatedata(0--1)
x=[2.10000.2000
3.10000.2000
4.10000.2000
6.10000.2000
7.10000.2000
2.10001.0000
3.10001.0000
4.10001.0000
6.10001.0000
7.10001.0000
2.10001.8000
3.10001.8000
4.10001.8000
6.10001.8000
7.10001.8000
2.10002.6000
3.10002.6000
4.10002.6000
6.10002.6000
7.10002.6000
]%(x-0.5)*1.5*2%definexto-1.5--0.5;
[xn,xs]=mapminmax(x);
xn=xn';
%x1=x(1,:
);
%x2=x(2,:
);
f=[21.3000
22.2000
23.4000
21.3000
21.3000
17.5000
19.2000
25.4000
17.6000
17.5000
17.0000
23.7000
29.0000
22.1000
17.0000
25.1000
28.4000
29.2000
27.6000
25.1000
]%10+x1.^2.*cos(2*pi.*x1)+x2.^2.*cos(2*pi.*x2);
f=f';
%subplot(2,1,1)
%plot(x,f)
%holdon
%subplot(2,1,2)
plot(xn,f)
title('initialdata')
%generateneuralnetwork
net=newrb(xn,f)
%x1=1;
%x2=1;
%a=10+x1.^2.*cos(2*pi.*x1)+x2.^2.*cos(2*pi.*x2)
h=sim(net,xn)
%subplot(2,1,1)
plot(x,h)
title('forecastfunction')
%holdon
、////////////////////////////////////////////////////////////////////////////////////////////////////////
clear
clc
closeall
x=0:
20:
500;
y=[32003400360038004015401540154015401540154015401540154015401540154015401540154015401540154015401540154015];
f=polyfit(x,y,2);
x0=0:
20:
500;
y0=polyval(f,x0);
plot(x,y,'*',x0,y0);
xlabel('横坐标x')
ylabel('纵坐标y')
title('描点绘制曲线')
legend('描点坐标','绘制的曲线')
gridoff
f1=5;f2=100;
T=10;
f0=(f2+f1)/2;
f3=f2-f1;
t=0:
0.1:
10;
s=sin(2*pi*(f0+(f3*t)/(2*T)).*t);
plot(t,s)
functionK=f1(Nx,i,j,Yj,Y1,Y2,Y3)
K1=0;
K2=1.51;
K3=0.93;
K4=2.91;
ifj
K=K1;
else
ifYj(i) K=K2; end ifYj(i) K=K3; end ifYj(i)>Y2 K=K4; end End functionS=f2(Nx,i,j,Xw,Xj,Yw,Yj,dl,K,MN) R=0.01; h1=2074; Tf1=50.00000; ifNx(i)==j dt=abs(R-(Xj(j)^2+Yj(i)^2)^0.5); dV=(Xw(j)-Xw(j-1))*(Yw(i)-Yw(i-1)); Sc=pi/2/MN*R*h1*K*Tf1/(K+h1*dt)/dV; Sp=-pi/2/MN*R*h1*K/(K+h1*dt)/dV; else Sc=0; Sp=0; end S=[Sc;Sp]; functionHeat4%边界不规则 R=0.01;X=0.125;Y1=0.03;Y_1=-0.03;Y2=0.05;Y3=0.07; MN=25;N1=115;M1=20;M_1=20;M2=20;M3=20; K1=0;K2=1.51;K3=0.93;K4=2.91;h1=2309;Tf1=50.00000;h2=60;Tf2=18.00000; dl=pi/2/MN*R; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% fori=1: MN+1 Xw0(i)=R*sin((i-1)*pi/2/MN); end fori=1: 2*MN+1 Yw0(i)=-R*cos((i-1)*pi/2/MN); end Xw1=R: (X-R)/N1: X; Xw=[Xw0,Xw1(2: N1+1)]; Yw_1=Y_1: (-R-Y_1)/M_1: -R; Yw1=R: (Y1-R)/M1: Y1; Yw2=Y1: (Y2-Y1)/M2: Y2; Yw3=Y2: (Y3-Y2)/M3: Y3; Yw=[Yw_1,Yw0(2: 2*MN),Yw1,Yw2(2: M2+1),Yw3(2: M3+1)]; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% M=length(Yw)-1; N=length(Xw)-1; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Xj (1)=Xw (1); forj=2: N+1 Xj(j)=(Xw(j)+Xw(j-1))/2; end Xj(N+2)=Xw(N+1); Yj (1)=Yw (1); fori=2: M+1 Yj(i)=(Yw(i)+Yw(i-1))/2; end Yj(M+2)=Yw(M+1); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% fori=1: M+2 ifi<=M_1+1 Nx(i)=0; end ifi>M_1+1&i<=M_1+MN+1 Nx(i)=i-M_1; end ifi>M_1+MN+1&i Nx(i)=M_1+2*MN+3-i; end ifi>M_1+2*MN+1 Nx(i)=0; end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% forj=2: N+2 iii=1; fori=1: M+2 ifNx(i)==j Ny(j,iii)=i; iii=iii+1; end end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% T=ones(M+2,N+2)*20.00000; forii=1: 3000 ii fori=2: M+1 forj=2: N+1 Kp=f1(Nx,i,j,Yj,Y1,Y2,Y3); Kw=f1(Nx,i,j-1,Yj,Y1,Y2,Y3); Ke=f1(Nx,i,j+1,Yj,Y1,Y2,Y3); Ks=f1(Nx,i-1,j,Yj,Y1,Y2,Y3); Kn=f1(Nx,i+1,j,Yj,Y1,Y2,Y3); S=f2(Nx,i,j,Xw,Xj,Yw,Yj,dl,Kp,MN); Sc=S (1); Sp=S (2); if(Kp+Kw)==0 Aw(j)=0; else Aw(j)=(Yw(i)-Yw(i-1))*Kw*Kp/(Kw*(Xj(j)-Xj(j-1))/2+Kp*(Xj(j)-Xj(j-1))/2); end if(Kp+Ke)==0 Ae(j)=0; else Ae(j)=(Yw(i)-Yw(i-1))*Ke*Kp/(Ke*(Xj(j+1)-Xj(j))/2+Kp*(Xj(j+1)-Xj(j))/2); end if(Kp+Ks)==0 As(j)=0; else As(j)=(Xw(j)-Xw(j-1))*Ks*Kp/(Ks*(Yj(i)-Yj(i-1))/2+Kp*(Yj(i)-Yj(i-1))/2); end if(Kp+Kn)==0 An(j)=0; else An(j)=(Xw(j)-Xw(j-1))*Kn*Kp/(Kp*(Yj(i+1)-Yj(i))/2+Kn*(Yj(i+1)-Yj(i))/2); end Ap(j)=-Aw(j)-Ae(j)-As(j)-An(j)+Sp*(Yw(i)-Yw(i-1))*(Xw(j)-Xw(j-1)); B(j)=-Sc*(Yw(i)-Yw(i-1))*(Xw(j)-Xw(j-1))-As(j)*T(i-1,j)-An(j)*T(i+1,j); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%上下边界条件 ifi==2 Ap(j)=-Aw(j)-Ae(j)-An(j)+Sp*(Yw(i)-Yw(i-1))*(Xw(j)-Xw(j-1)); B(j)=-Sc*(Yw(i)-Yw(i-1))*(Xw(j)-Xw(j-1))-An(j)*T(i+1,j); end ifi==M+1 KK=K4;%2.7; Ap(j)=Ap(j)+An(j)*KK/((Yj(M+2)-Yj(M+1))*h2+KK); B(j)=-Sc*(Yw(i)-Yw(i-1))*(Xw(j)-Xw(j-1))-As(j)*T(i-1,j)-An(j)*Tf2*h2*(Yj(M+2)-Yj(M+1))/((Yj(M+2)-Yj(M+1))*h2+KK); end end %%%%%%%%%%%%%%%%%%%%%%%%左右边界条件 Ap (2)=Ap (2)+Aw (2); Ap(N+1)=Ap(N+1)+Ae(N+1); ifNx(i)==0 A=diag(Ap(2: N+1))+diag(Aw(3: N+1),-1)+diag(Ae(2: N),1); T1=A\B(2: N+1)'; T(i,2: N+1)=T1'; else A=diag(Ap(Nx(i): N+1))+diag(Aw(Nx(i)+1: N+1),-1)+diag(Ae(Nx(i): N),1); T2=A\B(Nx(i): N+1)'; T(i,Nx(i): N+1)=T2'; T(i,2: Nx(i)-1)=Tf1;%20; end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% forj=2: N+1 fori=2: M+1 Kp=f1(Nx,i,j,Yj,Y1,Y2,Y3); Kw=f1(Nx,i,j-1,Yj,Y1,Y2,Y3); Ke=f1(Nx,i,j+1,Yj,Y1,Y2,Y3); Ks=f1(Nx,i-1,j,Yj,Y1,Y2,Y3); Kn=f1(Nx,i+1,j,Yj,Y1,Y2,Y3); S=f2(Nx,i,j,Xw,Xj,Yw,Yj,dl,Kp,MN); Sc=S (1); Sp=S (2); if(Kp+Kw)==0 Aw(i)=0; else Aw(i)=(Yw(i)-Yw(i-1))*Kw*Kp/(Kw*(Xj(j)-Xj(j-1))/2+Kp*(Xj(j)-Xj(j-1))/2); end if(Kp+Ke)==0 Ae(i)=0; else Ae(i)=(Yw(i)-Yw(i-1))*Ke*Kp/(Ke*(Xj(j+1)-Xj(j))/2+Kp*(Xj(j+1)-Xj(j))/2); end if(Kp+Ks)==0 As(i)=0; else As(i)=(Xw(j)-Xw(j-1))*Ks*Kp/(Ks*(Yj(i)-Yj(i-1))/2+Kp*(Yj(i)-Yj(i-1))/2); end if(Kp+Kn)==0 An(i)=0; else An(i)=(Xw(j)-Xw(j-1))*Kn*Kp/(Kp*(Yj(i+1)-Yj(i))/2+Kn*(Yj(i+1)-Yj(i))/2); end Ap(i)=-Aw(i)-Ae(i)-As(i)-An(i)+Sp*(Yw(i)-Yw(i-1))*(Xw(j)-Xw(j-1)); B(i)=-Sc*(Yw(i)-Yw(i-1))*(Xw(j)-Xw(j-1))-Aw(i)*T(i,j-1)-Ae(i)*T(i,j+1); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%左右边界条件 ifj==2 Ap(i)=Ap(i)+Aw(i);%-Ae(i)-As(i)-An(i)+Sp*(Yw(i)-Yw(i-1))*(Xw(j)-Xw(j-1)); B(i)=-Sc*(Yw(i)-Yw(i-1))*(Xw(j)-Xw(j-1))-Ae(i)*T(i,j+1); end ifj==N+1 Ap(i)=Ap(i)+Ae(i);%-Aw(i)-As(i)-An(i)+Sp*(Yw(i)-Yw(i-1))*(Xw(j)-Xw(j-1)); B(i)=-Sc*(Yw(i)-Yw(i-1))*(Xw(j)-Xw(j-1))-Aw(i)*T(i,j-1); end end %%%%%%%%%%%%%%%%%%%%%%%%上下边界条件 Ap (2)=Ap (2)+As (2); KK=K4;%2.7; Ap(M+1)=Ap(M+1)+An(M+1)*KK/((Yj(M+2)-Yj(M+1))*h2+KK); B(M+1)=B(M+1)-An(M+1)*Tf2*h2*(Yj(M+2)-Yj(M+1))/((Yj(M+2)-Yj(M+1))*h2+KK); ifj<=length(Ny) A1=diag(Ap(2: Ny(j,1)))+diag(As(3: Ny(j,1)),-1)+diag(An(2: Ny(j,1)-1),1); T3=A1\B(2: Ny(j,1))'; A2=diag(Ap(Ny(j,2): M+1))+diag(As(Ny(j,2)+1: M+1),-1)+diag(An(Ny(j,2): M),1); T4=A2\B(Ny(j,2): M+1)'; T(2: Ny(j,1),j)=T3; T(Ny(j,2): M+1,j)=T4; else A3=diag(Ap(2: M+1))+diag(As(3: M+1),-1)+diag(An(2: M),1); T5=A3\B(2: M+1)'; T(2: M+1,j)=T5; end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% A=0 fori=2: N+1 A=A+T(M+1,i) end t=A/N end%迭代循环 %K=2.7 %fori=1: M+2 %forj=1: N+2 %S=f2(Nx,i,j,Xw,Xj,Yw,Yj,dl,K,MN); %KK(i,j)=S (1); %KK1(i,j)=S (2); %end %end %KK1 surf(Xj(2: N+1),Yj(2: M+1),T(2: M+1,2: N+1)) shadinginterp view(0,90); figure,contour(Xj(2: N+1),Yj(2: M+1),T(2: M+1,2: N+1)); mydir=uigetdir('选择一个目录'); ifmydir(end)~='\' mydir=[mydir,'\']; end; files=dir([mydir,'*.bmp']); y=length(files); mkdir('C: \DocumentsandSettings\Administrator\桌面\','huidutuzhuose'); forx=1: y I=imread([mydir,files(x).name]); I1=rgb2gray(I) B=I1; [m,n]=size(I1); [L,num]=bwlabel(I1); stats=regionprops(L,'Area'); allArea=[stats.Area]; M=max(allArea) allArea=allArea/(M/60); B=double(B); %set('B',[111]); forl=1: num fora=1: m forb=1: n ifL(a,b)==l B(a,b)=B(a,b)-1+allArea(l); elseB(a,b)=B(a,b); end end end end figure,imagesc(B) daspect([111]);%将图片还原为原来大小 colorbar F=getframe(gcf); imwrite(F.cdata,['C: \DocumentsandSettings\Administrator\桌面\huidutuzhuose\',files(x).name]); %closeall; end function[f]=dfft(x,n) %dfft: 计算序列的n点fft并绘制图形; %x: 输入序列向量; %n: 序列fft的点数; %f: x的fft; f=fft(x,n); subplot(2,1,1) stem(0: length(x)-1,x,'.') title('信号原型') axis([0length(x)-1min(x)max(x)]) subplot(2,1,2) stem(0: n-1,abs(f),'.') title('上图序列FFT后的图象') end 调用: a=1: 5 dfft(a,5) 自组织竞争神经网络: loadabs.txt; data=abs; data=data(: 2: 115); p1=data(1: 40,: ); t1=data(41: 60,: ); p1=p1'; t1=t1'; %DataNormalization [pn,ps]=mapminmax(p1); [tn,ts]=mapminmax(t1); Q=minmax(pn); %GenerateNeuralNetwork net=newc(Q,2,0.01) net=init(net) net.trainParam.epochs=300; %Training net=train(net,pn) %Texting a=sim(net,pn) 建模2012: %%通道3温度和风速 clc clear closeall x=[0.30.91.52.12.7]; y=[2.457.2]; Z=[2729293029 3029313230 2731315231]; ZN=[0.40.60.70.80.9 0.40.50.60.70.6 0.40.60.60.60.5]; %subplot(2,1,1);mesh(x,y,z);title('通道3温度分布'); %%温度分布 xi=0.3: 0.01: 2.7;yi=2: 0.01: 7.5; [XI,YI]=meshgrid(xi,yi); %ZI=INTERP2(x,y,Z,XI,YI,'cubic'); ZI=interp2(x,y,Z,XI,YI,'cubic'); figure mesh(XI,YI,ZI);title('热通道3温度分布'); xlabel('高度z') ylabel('离空调的距离y') %%风速分布 xi=0.3: 0.01: 2.7; yi=2: ... 0.01: 7.5; [XI,YI]=meshgrid(xi,yi); ZI=griddata(x,y,ZN,XI,YI,'cubic'); figure mesh(XI,YI,ZI);title('热通道3风速分布'); xlabel('高度z') ylabel('离空调的距离y') clc clear closeall a=[1.12.10.221.321.321.321.3 1.13.10.222.222.222.222.2 1.14.10.223.423.623.423.4 1.16.10.221.321.321.321.3 1.17.10.221.321.321.321.3 1.12.1117.517.517.517.5 1.13.1119.219.219.219.2 1.14.1125.426.425.425.4 1.16.1117.617.617.617.6 1.17.1117.517.517.517.5 1.12.11.817.017.017.017.0 1.13.11.823.723.824.323.7 1.14.11.829.031.529.029.0 1.16.11.822.122.122.122.1 1.17.11.817.017.017.017.0 1.12.12.625.125.125.125.1 1.13.12.628.428.929.628.4 1.14.12.629.231.129.7
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 个人 收集 matlab 代码 截止 0528