① 请教编程高手，如何产生从-5到5的一组正态分布的随机数
#include <stdio.h>
//产生任意均值与方差的正态分布的随机数
// double u μ 正颤空春态分布的均值
// double g σ2=g2 正态分布的方差
// double * r 指向随机数种子
// double grn1() 返回一个均值u 方差g2 正态分布 的随机数
double grn1(double u,double g,double * r)
{ int i,m;
double s,w,v,t;
s=65536.0; w=2053.0; v=13849.0;
t=0.0;
for (i=1; i<=12; i++)
{ *r=(*r)*w+v; m=(int)(*r/s);
*r=*r-m*s; t=t+(*r)/s;
}
t=u+g*(t-6.0);
return(t);
}
int main()
{
double u,g,r;
r=5.0;u=1.0;g=1.5;
printf("\n");
for(int i=0;i<=9;i++){
for(int j=0;j<茄耐=4;j++){
printf("%10.7lf ", grn1(u,g,&r));
}
printf("\n");
}
printf("\n");
}
运行结果:
1.2386322 -1.1779938 0.5128021 1.9047699 1.5916595
1.1672211 -0.7747955 1.8593597 1.6634369 1.7311859
0.6563568 -1.4673004 -0.0460358 -0.9860992 0.3062592
3.9247894 -0.5367584 3.5153656 -0.3250885 0.0356293
0.1912689 0.2355804 0.2623138 0.3652191 3.6380463
1.1745453 2.0684662 0.4135590 2.3035736 4.8322601
-0.9066315 0.1806488 2.1878510 2.2087250 1.8370209
-0.3335114 3.2908783 0.8039398 1.2994232 -1.1289215
4.1126556 0.6179047 0.4805756 0.7944183 -1.3468170
1.6506195 2.3804779 2.4365082 1.9124603 0.9020844
参考 常用算法程序集(c语言描述)第三版
第三章 随机数的产生
3.5 产生任意均值与方差的正态分布亏数的一个随机数
② 最小二乘法三次多项式曲线拟合 算法 C++ 实现，该怎么处理
/******************************************
//参考《常用算法程序集（C语言描述第三版)》
//最小二乘法
//x[n]y[n]已知输入
//n输入点个数
//a[m]返回m-1次拟合多项式的m个系数
//m拟合多项式的项数,即拟合多项式的最高次为m-1
//dt[3]dt[0]返回拟合多项式与各数据点误差的平方和，
dt[1]返回拟合多项式与各数据点误差的绝对值之和
dt[2]返回拟合多项式与各数据点误差的绝对值的最大值
//
//拟合多项式的输出
//Y(x)=a0+a1(x-X)+a2(x-X)^2+……am(x-X)^m
//其中X为已知点x的平均值
******************************************/
#include"math.h"
voidpir1(x,y,n,a,m,dt)
intn,m;
doublex[],y[],a[],dt[];
{
inti,j,k;
doublez,p,c,g,q,d1,d2,s[20],t[20],b[20];
for(i=0;i<=m-1;i++)a[i]=0.0;
if(m>n)m=n;
if(m>20)m=20;
z=0.0;
for(i=0;i<=n-1;i++)z=z+x[i]/(1.0*n);
b[0]=1.0;d1=1.0*n;p=0.0;c=0.0;
for(i=0;i<=n-1;i++)
{p=p+(x[i]-z);c=c+y[i];}
c=c/d1;p=p/d1;
a[0]=c*b[0];
if(m>1)
{t[1]=1.0;t[0]=-p;
d2=0.0;c=0.0;g=0.0;
for(i=0;i<=n-1;i++)
{q=x[i]-z-p;d2=d2+q*q;
c=c+y[i]*q;
g=g+(x[i]-z)*q*q;
}
c=c/d2;p=g/d2;q=d2/d1;
d1=d2;
a[1]=c*t[1];a[0]=c*t[0]+a[0];
}
for(j=2;j<=m-1;j++)
{s[j]=t[j-1];
s[j-1]=-p*t[j-1]+t[j-2];
if(j>=3)
for(k=j-2;k>=1;k--)
s[k]=-p*t[k]+t[k-1]-q*b[k];
s[0]=-p*t[0]-q*b[0];
d2=0.0;c=0.0;g=0.0;
for(i=0;i<=n-1;i++)
{q=s[j];
for(k=j-1;k>=0;k--)
q=q*(x[i]-z)+s[k];
d2=d2+q*q;c=c+y[i]*q;
g=g+(x[i]-z)*q*q;
}
c=c/d2;p=g/d2;q=d2/d1;
d1=d2;
a[j]=c*s[j];t[j]=s[j];
for(k=j-1;k>=0;k--)
{a[k]=c*s[k]+a[k];
b[k]=t[k];t[k]=s[k];
}
}
dt[0]=0.0;dt[1]=0.0;dt[2]=0.0;
for(i=0;i<=n-1;i++)
{q=a[m-1];
for(k=m-2;k>=0;k--)
q=a[k]+q*(x[i]-z);
p=q-y[i];
if(fabs(p)>dt[2])dt[2]=fabs(p);
dt[0]=dt[0]+p*p;
dt[1]=dt[1]+fabs(p);
}
return;
}
③ 用c语言怎么编写输入一个矩阵求其逆矩阵的程序
下面是实现Gauss-Jordan法实矩阵求逆。
#include
<stdlib.h>
#include
<math.h>
#include
<stdio.h>
int
brinv(double
a[],
int
n)
{
int
*is,*js,i,j,k,l,u,v;
double
d,p;
is=malloc(n*sizeof(int));
js=malloc(n*sizeof(int));
for
(k=0;
k<=n-1;
k++)
{
d=0.0;
for
(i=k;
i<=n-1;
i++)
for
(j=k;
j<=n-1;
j++)
{
l=i*n+j;
p=fabs(a[l]);
if
(p>d)
{
d=p;
is[k]=i;
js[k]=j;}
}
if
(d+1.0==1.0)
{
free(is);
free(js);
printf("err**not
inv
");
return(0);
}
if
(is[k]!=k)
for
(j=0;
j<=n-1;
j++)
{
u=k*n+j;
v=is[k]*n+j;
p=a[u];
a[u]=a[v];
a[v]=p;
}
if
(js[k]!=k)
for
(i=0;
i<=n-1;
i++)
{
u=i*n+k;
v=i*n+js[k];
p=a[u];
a[u]=a[v];
a[v]=p;
}
l=k*n+k;
a[l]=1.0/a[l];
for
(j=0;
j<=n-1;
j++)
if
(j!=k)
{
u=k*n+j;
a[u]=a[u]*a[l];}
for
(i=0;
i<=n-1;
i++)
if
(i!=k)
for
(j=0;
j<=n-1;
j++)
if
(j!=k)
{
u=i*n+j;
a[u]=a[u]-a[i*n+k]*a[k*n+j];
}
for
(i=0;
i<=n-1;
i++)
if
(i!=k)
{
u=i*n+k;
a[u]=-a[u]*a[l];}
}
for
(k=n-1;
k>=0;
k--)
{
if
(js[k]!=k)
for
(j=0;
j<=n-1;
j++)
{
u=k*n+j;
v=js[k]*n+j;
p=a[u];
a[u]=a[v];
a[v]=p;
}
if
(is[k]!=k)
for
(i=0;
i<=n-1;
i++)
{
u=i*n+k;
v=i*n+is[k];
p=a[u];
a[u]=a[v];
a[v]=p;
}
}
free(is);
free(js);
return(1);
}
void
brmul(double
a[],
double
b[],int
m,int
n,int
k,double
c[])
{
int
i,j,l,u;
for
(i=0;
i<=m-1;
i++)
for
(j=0;
j<=k-1;
j++)
{
u=i*k+j;
c[u]=0.0;
for
(l=0;
l<=n-1;
l++)
c[u]=c[u]+a[i*n+l]*b[l*k+j];
}
return;
}
int
main()
{
int
i,j;
static
double
a[4][4]={
{0.2368,0.2471,0.2568,1.2671},
{1.1161,0.1254,0.1397,0.1490},
{0.1582,1.1675,0.1768,0.1871},
{0.1968,0.2071,1.2168,0.2271}};
static
double
b[4][4],c[4][4];
for
(i=0;
i<=3;
i++)
for
(j=0;
j<=3;
j++)
b[i][j]=a[i][j];
i=brinv(a,4);
if
(i!=0)
{
printf("塌早MAT
A
IS:
");
for
(i=0;
i<=3;
i++)
{
for
(j=0;
j<=3;
j++)
printf("%13.7e
",b[i][j]);
printf("
");
}
printf("
");
printf("闷埋MAT
A-
IS:
");
for
(i=0;
i<=3;
i++)
{
for
(j=0;
j<=3;
j++)
printf("%13.7e
",a[i][j]);
printf("
");
}
printf("
");
printf("MAT
AA-
IS:
");
brmul(b,a,4,4,4,c);
for
(i=0;
i<=3;
i++)
{
for
(j=0;
j<=3;
j++)
printf("%13.7e
",c[i][j]);
printf("
");
}
}
}参考资料：蚂衫蚂C常用算法程序集－徐士良
④ 请为C++新手推荐几本教程好书吧。我不怕学着难，只要经典有价值的
你自己挑几本看吧，
入门推荐看
C Primer Plus中文版（第五版）
C程序设计（第三版）
学好C之后再去学c++，建议看
白话c++
你先把这些书复制下来，这些书都是经典中的经典了。都是非常好的。
C程序设计语言（第2版·新版）---C语言“倚天屠龙双剑”---Brian W.Kernighan“C语言之父”
C Primer Plus中文版（第五版）--------C语言“倚天屠龙双剑”---Stephen Prata
C程序设计（第三版）---------------------------谭浩强
C语言大全（第四版）---------------------------HERBERT SCHILDT
C语言接口与实现：创建可重用软件的技术-------------DAVID R.HANSON
C语言参考手册(原书第5版)--------------------------Samuel P.Harbison
C程序设计教程---------------------------------H.M.Deitel/P.J.Deitel
C陷阱与缺陷-----------------------------------Andrew Koenig
The C Programming Language by Kernighan & Ritchie Complement : The Development of the C Language
C Traps and Pitfalls by Andrew Koenig
Expert C Programming by Peter van der Linden
The Standard C Library by P J Plauger
Books on Linux & C :
Linux Application Development by Michael K. Johnson & Erik W. Troan
Linux device drivers by Alessandro Rubini
Developing Linux applications with GTK+ and GDK (ISBN 0-7357-0021-4, New Riders)
C/C++入门基础
Windows C 程序设计入门与提高
单片机C语言入门
C++ 入门基础教程
C语言常用算法源代码
C常用算法程序集
C++ 语言命令详解
Borland C++使用手册
C语言编程基础
C语言编程要点
C语言数值算法程序大全(第二版)
Borland C++4.0安装与使用
C++ Builder 5.0经典入门
非C程序员C++入门自学教程
C/C++指令字典
C语言常用函数手册
C++标准程序库自修教程与参考手册
C语言图形函数介绍
C语言数值计算初步
C++ 设计新思维
C++Builder基础进阶
C语言初学者入门讲座
C++完全参考手册
C++Builder英文帮助
《C++Builder学习大全》中文版
C++ Builder编程研究文档
C游戏编程从入门到精通
C++程序调试实用手册
Borland C++4.0从入门到精通
实用教程：
经典C语言教程
C++捷径教程（第三版）
Borland C++实用教材
C++程序设计语言题解
C++上机实践指导教程（第三版）
创世纪的C++ Builder教程
Borland C++程序设计教程
C++ BUILDER 实用培训教程
C++程序设计培训教程
C++编程金典(第三版)
C高级实用程序设计
嵌入式系统中C的开发
C++Builder自学培训教程
C++大学教程(第二版)
C++傻瓜教程
精通C++ Builder 5程序设计高级教程
C++程序设计语言经典题与实验指导
C++程序设计教程实验书
C++大学基础教程
C++面向对象到窗口程序设计
C语言核心编程
嵌入式系统的C程序设计
经典C教程
实用C语言详解
Borland C++深入编程
Borland C++4.0程序设计
C++ 面向对象多线程编程
编程实例
C++实例教程
C++30天极限教程
C++Builder程序设计范例
C语言程序基础及应用实例
C语言趣味程序百例精解
C++编程实例详解
C++ Builder网络开发实例
C++实践之路
21天学通C语言(第六版)
C语言学习300例
C++Builder6.0界面实例开发经典编程900例(C语言)
一个月挑战C++
Borland C++应用篇
自学C++半月通
边用边学C语言
C语言完美演绎
设计60系列C++游戏
C语言程序设计题典
C++物件模型
C++ Builder 5编程实例与技巧
C语言程序设计及应用实例
Borland C++ 高级开发实践
Borland C++ Builder6开发人员
经验技巧：C++编码规范与指导
C++核心编程技术
C++程序设计陷阱
C数值算法程序大全
C和C++代码精粹
C++编程-数据结构与程序设计方法
C++语言教程大全
C++Builder核心program
C高级编程技术
C语言常见问题集
C++ Builder 3 核心编程技术
Borland C++程序设计技术
用C实现面向对象
C专家编程
C++数据结构
C语言高级实用编程技巧
C程序设计与应用
C&C++深层探索
C++ 和面向对象的数值计算
用C++开发Web商用程序
C++代码设计与重用
C++程序员JAVA编程
Borland C++宝典
C语言最新编程技巧200例
C++语言学习经验集合
C语言编程宝典大全
C++ Builder高级开发人员技术指南
C++ 内存管理算法和实现
Turboc库函数源程序与参考大全
C++ Builder 5 开发人员指南
C++程序设计实验指导（钱能）
用C语言设计屏幕界面技术
C和C++代码锦囊
C语言学习和精华文摘
标准C++宝典
C_C++深层探索
Borland C++库函数详解
高质量C、C++编程指南
C++和面向对象的数值计算
C++、
C++程序设计语言（特别版)---c++八大金刚----Bjarne Stroustrup“C++之父”
C++ Primer (第3版)中文版----c++八大金刚---Stanley B.Lippman
C++ Primer (第4版)中文版----c++八大金刚---Stanley B.Lippman
C++标准程序库—自修教程与参考手册--c++八大金刚--Nicolai M.Josuttis
C++语言的设计和演化-----c++八大金刚----Bjarne Stroustrup“C++之父”
深度探索C++对象模型---c++八大金刚----Stanley B.Lippman
Essential C++中文版---c++八大金刚---Stanley B.Lippman
Effective C++中文版 2nd Edition-----c++八大金刚------Scott Meyers
More Effective C++中文版----c++八大金刚------Scott Meyers
C++编程思想（第2版） 第1卷：标准C++导引--------Bruce Eckel
C++编程思想（第2版）第2卷：实用编程技术 --------Bruce Eckel
C++程序设计--------------------------谭浩强
C++ 程序设计教程(第2版)--------------钱能
C++ Primer Plus(第五版)中文版---Stephen Prata
广博如四库全书The c++ programming language、c++ Primer
深奥如山重水复Inside the c++ object model
程序库大全The c++ standard libray
工程经验之积累Effective c++、More Effective c++、Exceptional c++
c++八大金刚：
1、Essentital c++---lippman---C++之父，旁枝暂略，主攻核心，轻薄短小，初学者
2、The c++ programming language----C++之父，技术权威，用词深峻，思想深远，c++网络全书代表，圣经。
3、c++ Primer----lippman---纵横书市十数年，c++最佳教本，c++网络全书代表。
4、Inside the c++ object model-----lippman----揭示c++底层，非常好，非常难。
5、Effective c++-----通过50个编程实例，展示专家经验，行文有趣，深处浅出。
6、More Effective c++－－－－通过35个编程实例，展示专家经验，行文有趣，深处浅出。
7、The c++ standard libray－－－c++标准库的网络全书。
8、设计模式：可复用面向对象软件的基础
⑤ C语言 求矩阵的逆
//源程序如下#include<stdio.h>
#include<conio.h>
#include<string.h>
#include<iostream.h>
#include<stdlib.h>
#include<math.h>
#define max 100void inputstyle(int *); //输入函数
void input(int **,int); //输入函数
long danx(int **,int);
int sgnx(int);
void martx(int **,int);int main(void)
{
int style=0,i=0;
int matrix[max][max],*p[max];
for(i=0;i<max;i++)*(p+i)=matrix[i]; //*(p+i)是指针，指向第i个字符串
char exit1=' ';
while(exit1!='E'&& exit1!='e'){ printf("求n阶矩阵的逆\n"); inputstyle(&style);
input(p,style);
printf("原矩阵为:\n");
for(i=0;i<style;i++){
for(int j=0;j<style;j++){
printf("%4d",matrix[i][j]);
}
printf("\n");
}
martx(p,style);
printf("\n");
printf("Exit=e Continue=Press any key\n");
cin>>exit1;
fflush(stdin);
printf("\n\n"); }
return(0);
} void input(int **p,int n){
for(int i=0;i<n;i++){
for(int j=0;j<n;j++){
printf("迟蚂输入矩阵(%d行,%d列)元素:",j+1,i+1);
*(*(p+j)+i)=0;
scanf("%d",*(p+j)+i);
fflush(stdin);
}
}
}void inputstyle(int *style){
do{
printf("码携埋输入矩阵n*n阶数n(0<n<%d):",max);
fflush(stdin);
scanf("隐告%d",style);
fflush(stdin);
}while(*style<=0 && *style>max);
}long danx(int **p,int n){
int i=0,j1=0,k1=0,j2=0,k2=0;
long sum=0;
int operate[max][max],*po[max];
for(i=0;i<max;i++)*(po+i)=operate[i]; if(n==1)return *(*(p+0)+0);
else{
for(i=0;i<n;i++){
for(j1=1,j2=0;j1<n;j1++,j2++){
k1=-1;k2=-1;
while(k2<n-1){
k1++;
k2++;
if(k1==i)k1++;
*(*(po+j2)+k2)=*(*(p+j1)+k1);
}
}
/*for(int i1=0;i1<n-1;i1++){
for(int h1=0;h1<n-1;h1++){
printf("(%d,%d)%d ",i1,h1,*(*(po+h1)+i1));
}
printf("\n");
}*/
sum+=*(*(p+0)+i) * sgnx(1+i+1) * danx(po,n-1);
}
return sum;
}
}int sgnx(int i){
if(i%2==0)return(1);
else return(-1);
}void martx(int **p,int n){
int i=0,j=0,j1=0,k1=0,j2=0,k2=0,num=0;
int tramform[max][max];
int operate[max][max],*po[max];
for(i=0;i<max;i++)*(po+i)=operate[i];
num=danx(p,n);
if(num==0)printf("矩阵不可逆\n");
else{
if(n==1)printf("矩阵的逆为: 1/%d\n",num);
else{
printf("矩阵的逆为: 系数 1/%d *\n",num);
for(i=0;i<n;i++){
for(j=0;j<n;j++){
j1=-1;j2=-1;
while(j2<n-1){
j1++;j2++;
if(j1==j)j1++; k1=-1;k2=-1;
while(k2<n-1){
k1++;
k2++;
if(k1==i)k1++;
*(*(po+j2)+k2)=*(*(p+j1)+k1);
}
}
tramform[i][j]=sgnx(2+i+j) * danx(po,n-1);
}
}
for(i=0;i<n;i++){
for(j=0;j<n;j++){
printf("%4d",tramform[i][j]);
}
printf("\n");
}
}
}
}
//运行结果//希望对你有帮助