YD-602S 源代码

function c = reedcf(yd,lat,dsw); % REEDCF: computes daily mean cloud cover following Reed (1977). % c = REEDCF(yd,lat,dsw) com…

function yd=yearday(mon,day,leapyr) % YEARDAY: converts calender month and day into yearday. % yd = YEARDAY(mon,day,leapyr) co…

function sradav = clskswr(yd,lat) % CLSKSWR: computes clear sky insolation following Seckel&Beaudry (1973). % sradav = CLSKSWR…

function [x,k]=niudunfa(f1,f2,x0,ep) % 求解非线性方程的牛顿法，其中 % f1(x) --- 需要求根的函数。 % f2(x) --- f1(x)的导数。 % x0 --- 初始点。 % ep -…

%test 开始测试 count=0 load('-ascii', 'usps_test.txt');%载入测试集文件 for isamp=1:2000 for i=2:inNum+1 x(1,i-1)=usps_test(isa…

function Yd=DYDt2(t,Y,flag,G,ME) % flag 按ODE文件格式规定，必须是第三输入宗量。对它的赋值由ode45指令自动产生。 % 第4、5宗量是被传递的参数 switch flag case '' X=Y(1:2)…

#--------------------------------------------------- # Makefile.machine # This file defines machine dependencies #-------------…

%Pi-Sigma Type Fuzzy Neural Network clear all; close all; xite=0.85; alfa=0.05; p0=0.1*ones(9,1); p1=0.1*ones(9,1); p2=…

function Yd=DYDt2(t,Y,flag,G,ME) % flag 按ODE文件格式规定，必须是第三输入宗量。对它的赋值由ode45指令自动产生。 % 第4、5宗量是被传递的参数 switch flag case '' X=Y(1:2)…

function y = DEimpEuler1(f, h,a,b,y0) format long; N = (b-a)/h; y = zeros(N+1,1); x = a:h:b-h; y(1) = y0; var = findsym(f); …