亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频

? 歡迎來到蟲蟲下載站！ | ?? 資源下載 ?? 資源專輯 ?? 關于我們
? 蟲蟲下載站

?? func.cpp

?? c語言編寫數學常用算法
?? CPP
字號: 
#include <values.h>
#include <math.h>
#include <stdio.h>
#include "func.h"


DOUBLE algo::cal(DOUBLE x) { // 基類的基本算法
	return yfactor*calculate(xfactor*(x-xshift))+addconst;
}

algo * algo::clone() // 克隆自己，必須被繼承子類改寫
{
	return new algo(*this);
}

algo * algo::mul(DOUBLE a)	// 乘a
{
	algo * alg;
	alg = this;
	if(refnum>1) { refnum--;
		alg = clone();	// 如引用數大于1，則產生新的對象
	}
	alg->yfactor *= a;
	alg->addconst *= a;
	return alg;
}

algo * algo::add(DOUBLE a)	// 加a
{
	algo * alg;
	alg = this;
	if(refnum>1) { refnum--;
		alg = clone();	// 如引用數大于1，則產生新的對象
	}
	alg->addconst += a;
	return alg;
}

algo * algo::neg() // 取負
{
	algo * alg;
	alg = this;
	if(refnum>1) { refnum--;
		alg = clone();	// 如引用數大于1，則產生新的對象
	}
	alg->addconst = -alg->addconst;
	alg->yfactor = -alg->yfactor;
	return alg;
}

algo * algo::setxfactor(DOUBLE x)		// 設置x軸因子
{
	algo * alg;
	alg = this;
	if(refnum>1) { refnum--;
		alg = clone();	// 如引用數大于1，則產生新的對象
	}
	alg->xfactor = x;
	return alg;
}

algo * algo::xroom(DOUBLE x)	// 將xfactor擴大x倍
{
	algo * alg;
	alg = setxfactor(x*xfactor);
	if(x!=0)alg->xshift/=x;
	return alg;
}

algo * algo::setxshift(DOUBLE x) // 設置xshift的值
{
	algo * alg;
	alg = this;
	if(refnum>1) { refnum--;
		alg = clone();	// 如引用數大于1，則產生新的對象
	}
	alg->xshift = x;
	return alg;
}

algo * algo::xshiftas(DOUBLE x)	// 從當前開始平移x
{
	return setxshift(xshift+x);
}

algo * algojoin::clone() // 克隆自己
{
	return new algojoin(*this);
}

DOUBLE algojoin::calculate(DOUBLE x)	// 實施結合算法
{
	if(leftalgo==0 || rightalgo==0)
		throw TMESSAGE("empty algo pointor!");
	DOUBLE a,b;
	if(met == ccom)		// 復合函數
		return leftalgo->cal(rightalgo->cal(x));
	a = leftalgo->cal(x);
	b = rightalgo->cal(x);
	if(met == cadd)	// 返回各結合運算值
		return a+b;
	else if(met == csub)
		return a-b;
	else if(met == cmul)
		return a*b;
	else if(met == cdiv)
		return a/b;
	else if(met == cpow)
		return pow(a,b);
	return 0.0;
}

algo * algofun::clone() // 克隆自己
{
	return new algofun(*this);
}

DOUBLE algofun::calculate(DOUBLE x)	// 實施函數算法
{
	if(f)
		return f(x);
	return 0.0;
}

func::func()// 缺省構造函數，產生函數f(x)=x;
{
	alg = new algo();
}

func::func(func & fn) // 拷貝構造函數
{
	alg = fn.alg;
	alg->refnum++;
}

func::func(DOUBLE (*fun)(DOUBLE))	// 函數指針的構造函數
{
	alg = new algofun(fun);
}

func::func(DOUBLE a) // 常函數構造函數
{
	alg = new algo(a);
}

func::func(algo * a) // 算法的構造函數
{
	if(a)
		alg = a;
}

func::func(algo& a) // 另一種算法構造函數
{
	alg = a.clone();	// 克隆出一個同樣的算法變量并將指針賦予alg.
}

func& func::operator=(func& fn)	// 賦值運算符
{
	if(this == &fn) return (*this); // 如果等于自己，干脆啥也別做
	if(alg) {
		alg->refnum--;	// 原算法引用數減1
		if(!alg->refnum)	// 如結果為0則刪除原算法
			delete alg;
	}
	alg = fn.alg;	// 將fn的算法移過來
	if(alg) alg->refnum++; // 引用數增加
	return (*this);
}

func& func::operator=(DOUBLE (*fn)(DOUBLE)) // 用函數指針的賦值運算符
{
	if(alg) {
		alg->refnum--;	// 原算法引用數減1
		if(!alg->refnum)	// 如結果為0則刪除原算法
			delete alg;
	}
	alg = new algofun(fn);
	return (*this);
}

func& func::operator=(DOUBLE a) // 常函數的賦值運算符
{
	if(alg) {
		alg->refnum--;	// 原算法引用數減1
		if(!alg->refnum)	// 如結果為0則刪除原算法
			delete alg;
	}
	alg = new algo(a);
	return (*this);
}

func& func::operator+=(func& fn) // 自身加一個函數
{
	algo * a = new algojoin(alg, fn.alg, cadd);
	alg->refnum--;	// 因為聯合算法對兩個算法都加了引用數，因此再減回來
	alg = a;	// 將新指針賦給alg
	return (*this);
}

func& func::operator+=(DOUBLE (*f)(DOUBLE)) // 自身加一個函數指針
{
	func fn(f);	// 將函數指針包裝為函數類
	operator+=(fn);	// 實施加操作
	return (*this);
}

func func::operator+(func& fn)	// 相加產生新函數
{
	algo * a = new algojoin(alg, fn.alg, cadd);
	func f(a);
	return f;
}

func func::operator+(DOUBLE a)	// 與常數相加產生新函數
{
	func f(*this);
	f += a;
	return f;
}

func func::operator+(DOUBLE (*f)(DOUBLE)) // 加一個函數指針產生新函數
{
	func ff(*this);
	ff += f;
	return ff;
}

func& func::neg() // 自身取負
{
	alg=alg->neg();
	return (*this);
}

func func::operator-() // 產生負函數
{
	func f(*this);
	f.neg();
	return f;
}

func& func::operator-=(func& fn) // 自身減一個函數
{
	algo * a = new algojoin(alg, fn.alg, csub);
	alg->refnum--;	// 因為聯合算法對兩個算法都加了引用數，因此再減回來
	alg = a;	// 將新指針賦給alg
	return (*this);
}

func& func::operator-=(DOUBLE (*f)(DOUBLE)) // 自身減一個函數指針
{
	func fn(f);	// 將函數指針包裝為函數類
	operator-=(fn);	// 實施減操作
	return (*this);
}

func func::operator-(func& fn)	// 相減產生新函數
{
	algo * a = new algojoin(alg, fn.alg, csub);
	func f(a);
	return f;
}

func func::operator-(DOUBLE a)	// 與常數相減產生新函數
{
	func f(*this);
	f -= a;
	return f;
}

func operator-(DOUBLE a, func& f) // 常數減函數
{
	return (-f)+a;
}

func func::operator-(DOUBLE (*f)(DOUBLE)) // 減一個函數指針產生新函數
{
	func ff(*this);
	ff -= f;
	return ff;
}

func operator-(DOUBLE (*f)(DOUBLE),func& fn) // 函數指針減函數
{
	func ff(f);
	ff -= fn;
	return ff;
}

func& func::operator*=(func& fn) // 自身乘一個函數
{
	algo * a = new algojoin(alg, fn.alg, cmul);
	alg->refnum--;	// 因為聯合算法對兩個算法都加了引用數，因此再減回來
	alg = a;	// 將新指針賦給alg
	return (*this);
}

func& func::operator*=(DOUBLE (*f)(DOUBLE)) // 自身乘一個函數指針
{
	func fn(f);	// 將函數指針包裝為函數類
	operator*=(fn);	// 實施乘操作
	return (*this);
}

func func::operator*(func& fn)	// 相乘產生新函數
{
	algo * a = new algojoin(alg, fn.alg, cmul);
	func f(a);
	return f;
}

func func::operator*(DOUBLE a)	// 與常數相乘產生新函數
{
	func f(*this);
	f *= a;
	return f;
}

func func::operator*(DOUBLE (*f)(DOUBLE)) // 乘一個函數指針產生新函數
{
	func ff(*this);
	ff *= f;
	return ff;
}

func& func::operator/=(func& fn) // 自身除以一個函數
{
	algo * a = new algojoin(alg, fn.alg, cdiv);
	alg->refnum--;	// 因為聯合算法對兩個算法都加了引用數，因此再減回來
	alg = a;	// 將新指針賦給alg
	return (*this);
}

func& func::operator/=(DOUBLE (*f)(DOUBLE)) // 自身除以一個函數指針
{
	func fn(f);	// 將函數指針包裝為函數類
	operator/=(fn);	// 實施除法操作
	return (*this);
}

func func::operator/(func& fn)	// 相除產生新函數
{
	algo * a = new algojoin(alg, fn.alg, cdiv);
	func f(a);
	return f;
}

func func::operator/(DOUBLE a)	// 與常數相除產生新函數
{
	func f(*this);
	f /= a;
	return f;
}

func operator/(DOUBLE a, func& f) // 常數除以函數
{
	func ff(a);
	return ff/f;
}

func func::operator/(DOUBLE (*f)(DOUBLE)) // 除以一個函數指針產生新函數
{
	func ff(*this);
	ff /= f;
	return ff;
}

func operator/(DOUBLE (*f)(DOUBLE),func& fn) // 函數指針除以函數
{
	func ff(f);
	ff /= fn;
	return ff;
}

void func::setxfactor(DOUBLE a)	// 設置x因子為a
{
	alg = alg->setxfactor(a);
}

void func::xroom(DOUBLE a)	  // x方向擴大a倍
{
	alg = alg->xroom(a);
}

func& func::power(func& f)	// 函數的f次乘冪，函數自身改變
{
	algo * a = new algojoin(alg, f.alg, cpow);
	alg->refnum--;	// 因為聯合算法對兩個算法都加了引用數，因此再減回來
	alg = a;	// 將新指針賦給alg
	return (*this);
}

func func::operator^(func& fn)	// f次乘冪，產生新函數
{
	algo * a = new algojoin(alg, fn.alg, cpow);
	func f(a);
	return f;
}

func& func::power(DOUBLE a)	// 函數的a次冪，函數自身改變
{
	func f(a);
	return power(f);
}

func func::operator^(DOUBLE a)  // 函數的a次冪，產生新函數，原函數不變
{
	func f(a);
	func ff(*this);
	return ff.power(f);
}

func func::operator()(func & fn)	// 復合函數，產生新的函數
{
	algo * a = new algojoin(alg, fn.alg, ccom);
	func f(a);
	return f;
}

DOUBLE algopoly::calculate(DOUBLE x)	// 實施函數算法
{
	size_t i;
	DOUBLE u;
	size_t n=data.rownum-1;
	u=data(n,0);
	for (i=n; i>0; i--)
		u=u*x+data(i-1,0);
	return u;
}

algo * algopoly::clone()	// 克隆自己
{
	return new algopoly(*this);
}

DOUBLE algogass::calculate(DOUBLE x)	// 實施函數算法
{
	return gass(0,1,x);
}

algo * algogass::clone() // 克隆自己
{
	return new algogass(*this);
}

func::func(DOUBLE a, DOUBLE d)	// 產生正態分布函數a為均值，d為標準差
{
	alg = new algogass(a,d);
}

func::func(matrix& m, funckind kind)	// 構造數值相關函數，
				// 如kind = polyfunc, 則m為nX1矩陣，是n-1階多項式系數，
				// 其中m(0,0)為常數項，m(n-1,0)為n-1次項。
{
	if(kind == polyfunc)
		alg = new algopoly(m);	// 產生多項式算法
	else if(kind == enter2func)
		alg = new algoenter2(m);	// 產生一元全區間不等距插值
}

DOUBLE algoenter2::calculate(DOUBLE x)	// 實施函數算法
{
	size_t n = data.rownum;
	size_t i,j,k,m;
	double z,s;
	z=0.0;
	if (n<1) return(z);
	if (n==1) { z=data(0,1); return z;}
	if (n==2) {
		z = data(0,1)*(x-data(1,0))-data(1,1)*(x-data(0,0))/
			(data(0,0)-data(1,0));
		  return(z);
	}
	i=0;
	while ((i<n) && (data(i,0)<x)) i++; // ((x[i]<t)&&(i<n)) i=i+1;
	if (i<4) k=0;
	else k=i-4;
	m=i+3;
	if (m>n-1) m=n-1;
	for (i=k;i<=m;i++)
	{	s=1.0;
		for (j=k;j<=m;j++)
			if (j!=i) s *= (x-data(j,0))/(data(i,0)-data(j,0)); //s=s*(t-x[j])/(x[i]-x[j]);
		z=z+s*data(i,1); // y[i];
	}
	return(z);
}

algo * algoenter2::clone()	// 克隆自己
{
	return new algoenter2(*this);
}

DOUBLE algoenter::calculate(DOUBLE x)	// 實施函數算法
{
// double eelgr(x0,h,n,y,t)
	size_t n=data.rownum;
//  double x0,h,t,y[];
	size_t i,j,k,m;
	DOUBLE z,s,xi,xj,p,q;
	z=0.0;
	if (n<1) return z;
	if (n==1) { z=data(0); return(z);}
	if (n==2)
	{	z=(data(1)*(x-x0)-data(0)*(x-x0-h))/h;
		return z;
	}
	if (x>x0)
	{	p=(x-x0)/h; q = floor(p); i=q;
		if (p>q) i++;
	}
	else i=0;
	k=i-4;
	if (i<4) k=0;
	else k=i-4;
	m=i+3;
	if (m>n-1) m=n-1;
	for (i=k;i<=m;i++)
	{	s=1.0; xi=x0+i*h;
		for (j=k; j<=m; j++)
			if (j!=i)
			{	xj=x0+j*h;
				s*=(x-xj)/(xi-xj);
			}
			z=z+s*data(i);
		}
	return z;
}

algo * algoenter::clone()	// 克隆自己
{
	return new algoenter(*this);
}

func::func(matrix& m, DOUBLE x0, DOUBLE h) // 構造等距插值函數
			// 其中m是nX1階矩陣，代表n個y值，x0是起始點，h是步長（采樣間隔）
{
	alg = new algoenter(m,x0,h);
}

void func::setxshift(DOUBLE a)	// 設置函數沿x軸平移a
{
	alg = alg->setxshift(a);
}

void func::shiftxas(DOUBLE a) // 函數沿x軸右移a
{
	alg = alg->xshiftas(a);
}

DOUBLE algo::integ(DOUBLE a, DOUBLE b, DOUBLE eps) // 在(a,b)區間積分
			// 采用勒讓德高斯求積法
{
	 size_t m,i,j;
	 DOUBLE s,p,ep,h,aa,bb,w,x,g;
	 static DOUBLE t[5]={-0.9061798459,-0.5384693101,0.0,
								 0.5384693101,0.9061798459};
	 static DOUBLE c[5]={0.2369268851,0.4786286705,0.5688888889,
								0.4786286705,0.2369268851};
	 m=1;
	 h=b-a; s=fabs(0.001*h);
	 p=1.0e+35; ep=eps+1.0;
	 while ((ep>=eps)&&(fabs(h)>s))
		{ g=0.0;
		  for (i=1;i<=m;i++)
			 { aa=a+(i-1.0)*h; bb=a+i*h;
				w=0.0;
				for (j=0;j<=4;j++)
				  { x=((bb-aa)*t[j]+(bb+aa))/2.0;
					 w=w+cal(x)*c[j];  //flrgsf(x)*c[j];
				  }
				g=g+w;
			 }
		  g *= h/2.0;
		  ep=fabs(g-p)/(1.0+fabs(g));
		  p=g; m=m+1; h=(b-a)/m;
		}
	 return g;
}

DOUBLE func::integ(DOUBLE a, DOUBLE b, DOUBLE eps)
		// 從a到b計算函數的定積分
{
	return alg->integ(a,b,eps);
}

DOUBLE func::singleroot(DOUBLE x, DOUBLE eps)
	// 計算函數在x附近的一個單根
{
	return getroot(*alg,x,eps);
}

DOUBLE algoinverse::calculate(DOUBLE x)	// 實施反函數算法
{
	algo * a = al->clone();	// 克隆一個一樣的函數算法變量
	a->add(-x); // 形成函數f(t)-x
	DOUBLE t;
	t = getroot(*a);	// 求f(t)-x=0的根
	delete a;
	return t;
}

algo * algoinverse::clone() // 克隆自己
{
	return new algoinverse(*this);
}

func func::inverse()	// 產生反函數
{
	algo * a = new algoinverse(alg);
	func f;
	f.alg = a;
	return f;
}

?? 快捷鍵說明

復制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
国产欧美一区二区精品性色| 免费看黄色91| 日本vs亚洲vs韩国一区三区 | 色婷婷综合中文久久一本| 欧美日韩精品一区视频| 国产精品美女久久久久久2018| 日韩福利视频导航| 成人av影院在线| 日韩精品专区在线影院重磅| 一区二区三区久久久| 成人18视频在线播放| 欧美成人在线直播| 亚洲国产精品久久久久秋霞影院 | 国产成人精品亚洲日本在线桃色| 欧美精品视频www在线观看| 国产精品午夜免费| 精品一区二区免费在线观看| 69堂成人精品免费视频| 亚洲情趣在线观看| 成人激情开心网| 久久免费的精品国产v∧| 国产盗摄一区二区三区| 欧美日韩亚洲国产综合| 亚洲三级小视频| 国产99久久久国产精品潘金| 久久综合久久久久88| 久久精品国产精品青草| 欧美一区二区精品在线| 亚洲第一久久影院| 欧美色图天堂网| 一区二区高清免费观看影视大全| 成人性色生活片| 欧美激情综合网| 国产a级毛片一区| 亚洲国产高清在线观看视频| 国产精品 日产精品 欧美精品| 欧美xxxxx裸体时装秀| 精品一区二区三区免费播放| 精品乱人伦小说| 狠狠色丁香九九婷婷综合五月| 日韩三级在线观看| 久草中文综合在线| 久久综合九色综合97婷婷女人| 国产一区二区伦理片| 国产日产精品1区| 国产成人激情av| 中文字幕永久在线不卡| 日本道免费精品一区二区三区| 亚洲综合色成人| 这里是久久伊人| 极品少妇一区二区三区精品视频 | 国产精品情趣视频| 99视频精品免费视频| 亚洲人吸女人奶水| 欧美日韩国产综合一区二区三区| 亚洲成a人v欧美综合天堂下载 | 中文字幕电影一区| jlzzjlzz欧美大全| 亚洲网友自拍偷拍| 日韩无一区二区| 成人夜色视频网站在线观看| 亚洲乱码中文字幕| 日韩一区二区在线观看视频| 国产精品一区在线| 亚洲精品中文在线影院| 欧美一二区视频| 不卡欧美aaaaa| 日韩av高清在线观看| 久久久久久免费毛片精品| 91小视频在线免费看| 亚洲v日本v欧美v久久精品| 亚洲精品一线二线三线| 99久久精品国产精品久久| 午夜私人影院久久久久| 久久免费视频色| 欧美日韩视频在线一区二区| 韩日av一区二区| 亚洲图片欧美视频| 国产亚洲制服色| 欧美三级日韩三级国产三级| 国产精品一区在线观看你懂的| 伊人性伊人情综合网| 337p日本欧洲亚洲大胆色噜噜| 91麻豆精品在线观看| 国内成人精品2018免费看| 亚洲激情图片qvod| 久久久精品国产免费观看同学| 欧美优质美女网站| 成人国产免费视频| 国产一区二区主播在线| 五月天精品一区二区三区| 国产精品日日摸夜夜摸av| 日韩欧美中文字幕制服| 精品日韩在线观看| 在线观看日韩电影| 9人人澡人人爽人人精品| 激情国产一区二区| 丝瓜av网站精品一区二区| 最新热久久免费视频| 久久久99精品免费观看不卡| 欧美一区二区视频免费观看| 91国偷自产一区二区开放时间| 国产大陆亚洲精品国产| 捆绑变态av一区二区三区| 亚洲福利国产精品| 国产精品久久久久久久久久免费看| 91精品国产麻豆国产自产在线| 91国产视频在线观看| a亚洲天堂av| 成人性生交大合| 国产成人高清在线| 国产69精品久久久久777| 九一久久久久久| 精品影视av免费| 老司机精品视频线观看86| 日韩av不卡在线观看| 日日夜夜精品免费视频| 亚洲综合色自拍一区| 亚洲一区精品在线| 一区二区三区在线播放| 亚洲免费色视频| 一区二区视频在线| 亚洲激情图片qvod| 洋洋成人永久网站入口| 一区二区三区四区蜜桃 | 99国内精品久久| 北条麻妃一区二区三区| av中文字幕亚洲| 91久久免费观看| 欧美日本乱大交xxxxx| 91麻豆精品国产91久久久资源速度| 欧美日韩国产高清一区二区三区 | 狠狠色综合色综合网络| 国产伦精品一区二区三区免费迷| 国产剧情一区二区三区| 成人黄色在线看| 欧美午夜电影在线播放| 在线播放视频一区| 欧美精品一区二区不卡 | 国产成人午夜精品影院观看视频| 国产成a人亚洲| 一本一本大道香蕉久在线精品| 精品视频一区二区三区免费| 日韩亚洲欧美在线| 国产欧美日韩不卡免费| 一区二区三区久久久| 青青草91视频| 成人免费电影视频| 欧美性生活久久| 精品国产一区久久| 亚洲视频1区2区| 精品一区二区免费在线观看| 99r国产精品| 欧美一区二区三级| 国产精品理论片在线观看| 亚洲午夜三级在线| 国产成人免费xxxxxxxx| 欧洲亚洲精品在线| 久久久久久久综合色一本| 一区二区激情视频| 国产精品资源在线看| 欧美色图一区二区三区| 国产日产欧美一区二区三区 | 中文字幕一区二区三区色视频 | 91一区在线观看| 在线观看91av| 中文字幕日韩一区| 男女性色大片免费观看一区二区| 成人性视频网站| 欧美一区二区在线免费观看| 国产精品三级av在线播放| 人人狠狠综合久久亚洲| 99热这里都是精品| 久久久久久免费毛片精品| 亚洲高清久久久| 91丨九色丨国产丨porny| 日韩一区二区三区电影| 亚洲黄色小视频| 亚洲美女淫视频| 亚洲欧美日韩久久| 国产一区二区h| 欧美精品乱人伦久久久久久| 国产精品久久久久久久久免费相片| 日产国产欧美视频一区精品| 91视频一区二区| 国产亚洲污的网站| 理论电影国产精品| 91精品欧美一区二区三区综合在| 亚洲精选视频在线| 国产69精品久久久久777| 欧美本精品男人aⅴ天堂| 日本不卡一区二区| 欧美群妇大交群中文字幕| 夜夜爽夜夜爽精品视频| 99国产精品久久久久| 国产精品视频第一区| 丁香婷婷综合色啪| 国产精品你懂的在线欣赏| 国产黄人亚洲片| 国产女主播一区|