本文目录一览:
- 1、用C语言编写流星雨程序
- 2、C语言潮流程序中宏定义#define Y_(i,j)(*(*(Y+i)-i+j))是什么意思?
- 3、经典C语言程序例子
- 4、FORTRAN编程语言与C语言
- 5、为什么我编的电力系统潮流计算C语言PQ分解法不收敛
用C语言编写流星雨程序
数字流星雨代码:
//流星雨.cpp:Defines the entry point for the console application.
//
///////////////////////////////////////////////////
//程序名称:数字流星雨
//最后修改:2006-10-15
///////////////////////////////////////////////////
#includelt;windows.hgt;
#includelt;time.hgt;
#includelt;stdlib.hgt;
#include"graphics.h"
#includelt;conio.hgt;
#includelt;math.hgt;
/***********************宏定义**********************/
#define PI 3.1415926//圆周率
#define WIDTH 200//屏幕宽度,流星出生区域
#define HEIGHT 150//屏幕高度,流星出生区域
#define V 20//流星速度,单次移动的像素数
#define LENGTH 20//流星字符数
#define DELAY 30//延时
#define NUM 45//流星个数
/******************定义流星结构体*******************/
struct meteor
{
int x0;
int y0;
int yh;
char str[LENGTH];
}me[NUM]={0};
/*********************函数声明**********************/
char AsciiRand();
void Move(char*p);
void InitMeteor(struct meteor*me);
int color(int y,int y0,int yh);
void Meteors(struct meteor me[]);
/***********************主函数**********************/
///int main(void)
int _tmain(int argc,_TCHAR*argv[]){
char c='\0';//接收键盘输入的变量
initgraph(WIDTH,HEIGHT);//初始化WIDTH*HEIGHT的绘图窗口
HWND hwnd=GetHWnd();//获得窗口句柄
SetWindowText(hwnd,"Gavin Liu数字流星雨");//修改窗口名称
ShowWindow(hwnd,SW_SHOWMAXIMIZED);//最大化显示窗口
MessageBox(hwnd,TEXT("点击【确定】开始演示流星雨效果,Esc键退出"),TEXT("提示"),MB_OK|MB_ICONWARNING);//弹出提示
srand((unsigned)time(NULL));//设置随机种子
for(int i=0;ilt;NUM;i++){//对NUM个流星体初始化
InitMeteor(melt;igt;);
}
while(c!=27){
BeginBatchDraw();//开始批量绘图
Meteors(me);//绘制一帧动画
FlushBatchDraw();//执行未完成的绘制任务
Sleep(DELAY);//延时
cleardevice();//清屏
for(int i=0;ilt;NUM;i++){
melt;igt;.yh+=V;
Move(melt;igt;.str);
if(melt;igt;.yhgt;HEIGHT+LENGTH*V){
InitMeteor(melt;igt;);
}
}
if(kbhit()){
c=getch();
}
}
EndBatchDraw();//结束批量绘图
closegraph();//结束绘图环境
return 0;
}
/***********************函数体**********************/
char AsciiRand(){//产生随机可见ASCII码
return((char)(rand()%(126-33)+33));
}
void Move(char*p){//字符后移,可以使显示时字符相对屏幕位置不变
char*pt=p+LENGTH;
while(ptgt;p){
*(--pt)=*(pt-1);
}
*p=AsciiRand();
}
void InitMeteor(struct meteor*me){//对一颗流星初始化
me-gt;x0=rand()%WIDTH;
me-gt;yh=me-gt;y0=rand()%HEIGHT;
for(int i=0;ilt;LENGTH;i++)
{
*(me-gt;str+i)=AsciiRand();
}
}
int color(int y,int y0,int yh){//确定流星的颜色
int color;
//出生点之前的流星体置成黑色
if(ylt;y0){
color=0;
}
//流星颜色自头至尾按照余弦函数递减
else{
//尾迹消失
color=(int)(255*cos((yh-y)*PI/(2*LENGTH*V)));
}
return color;
}
//打印一帧流星的画面
void Meteors(struct meteor me[]){
//设置格式:背景透明,字符高度,字体粗细,字体
setbkmode(TRANSPARENT);
setfont(12,12,"宋体");
//开始打印一帧图像
int y;
for(int n=0;nlt;NUM;n++){
for(int j=0;jlt;LENGTH;j++){
//流星中第j个字符的纵坐标
y=me[n].yh-j*V;
//设置颜色,流星的头部是白色的
setcolor(RGB(255*(0==j),color(y,me[n].y0,me[n].yh),255*(0==j)));
//打印字符
outtextxy(me[n].x0,y,me[n].str[j]);
}
}
}
扩展资料:
include用法:
#include命令预处理命令的一种,预处理命令可以将别的源代码内容插入到所指定的位置;可以标识出只有在特定条件下才会被编译的某一段程序代码;可以定义类似标识符功能的宏,在编译时,预处理器会用别的文本取代该宏。
插入头文件的内容
#include命令告诉预处理器将指定头文件的内容插入到预处理器命令的相应位置。有两种方式可以指定插入头文件:
1、#includelt;文件名gt;
2、#include"文件名"
C语言潮流程序中宏定义#define Y_(i,j)(*(*(Y+i)-i+j))是什么意思?
#define Y_(i,j) (*(*(Y+i)-i+j))
你把Y_(i,j)当成一个子函数,i , j 则为参数,
当你调用Y_( i , j )时,把参数i, j带入(*(*(Y+i)-i+j))这个公式中计算,计算结果即为返回值
经典C语言程序例子
题目01:在一个已知的字符串中查找最长单词,假定字符串中只含字母和空格,空格用来分隔不同的单词。
直接编译,程序执行结果如下图所示:
题目02:编写一个int string_len(char *s),返回字符串s的字符长度(不包括\0)。
直接编译,程序执行结果如下图所示:
扩展资料:
C语言是一门通用计算机编程语言,应用广泛。C语言的设计目标是提供一种能以简易的方式编译、处理低级存储器、产生少量的机器码以及不需要任何运行环境支持便能运行的编程语言。
尽管C语言提供了许多低级处理的功能,但仍然保持着良好跨平台的特性,以一个标准规格写出的C语言程序可在许多电脑平台上进行编译,甚至包含一些嵌入式处理器(单片机或称MCU)以及超级电脑等作业平台。
FORTRAN编程语言与C语言
fortran语言中可通过定义子程序从而返回一系列变量值,而c语言的函数一般只能返回一个变量值。
Fortran语言的最大特性是接近数学公式的自然描述,在计算机里具有很高的执行效率。易学,语法严谨。可以直接对矩阵和复数进行运算,这一点类似MATLAB。自诞生以来广泛地应用于数值计算领域,积累了大量高效而可靠的源程序。很多专用的大型数值运算计算机针对Fortran做了优化。广泛地应用于并行计算和高性能计算领域。
C语言是一个有结构化程序设计、具有变量作用域(variable scope)以及递归功能的过程式语言。C语言传递参数均是以值传递(pass by value),另外也可以传递指针(a pointer passed by value)。不同的变量类型可以用结构体(struct)组合在一起。只有32个保留字(reserved keywords),使变量、函数命名有更多弹性。部份的变量类型可以转换,例如整型和字符型变量。通过指针(pointer),C语言可以容易的对存储器进行低级控制。预编译处理(preprocessor)让C语言的编译更具有弹性。
为什么我编的电力系统潮流计算C语言PQ分解法不收敛
//PQ 分解法源程序
#include "math.h"
#include "stdio.h"
#include "stdlib.h"
#define M 14 //预计支路数
#define N 14 //预计节点数
//
struct linestr //支路信息
{
int i,j; //和该支路相连的两个节点
float r,x,g,b; //该支路的参数,kb为非标准变比
float pij,qij,pji,qji; //支路潮流—计算数据
}line[M];
struct pointstr //节点信息
{
int sign; //节点类型,0—平衡节点,1—PQ节点,2—PV节点
float voltage,cta,p,q; //节点电压幅值和相位角,节点注入有功、无功
}pointpri[N],pointaft[N]; //pointpri用于存放节点的初始数据,pointaft用于存放计算后的节点数据
float G[N][N]={0},B[N][N]={0}; //节点导纳阵
//
float B1[N][N]={0},B2[N][N]={0};//修正方程系数矩阵
float detp[N]={0},detct[N]={0},detq[N]={0},detv[N]={0};//依次放修正方程的修正量
float pn[N]={0},ctn[N]={0},qn[N]={0},vn[N]={0}; //依次放修正量对应的节点号,即变量对应的节点号
int nline,npoint; //准确的支路数和节点数
int npq,npv,nbal; //PQ节点数、PV节点数、平衡节点数
void creat()
{
int i;
npq=0;npv=0;nbal=0;
printf("\n请输入总支路数:");scanf("%d",nline);
for(i=1;i=nline;i++)
{
printf("请输入%d号支路数据:\n",i);
printf("和该支路相连的节点号i,j(请用逗号隔开):");scanf("%d,%d",(line[i].i),(line[i].j));
printf("该支路的电阻和电抗R,X(请用逗号隔开):");scanf("%f,%f",(line[i].r),(line[i].x));
printf("该支路的对地电导和电纳G,B(请用逗号隔开):");scanf("%f,%f",(line[i].g),(line[i].b));
line[i].pij=0;line[i].qij=0;line[i].pji=0;line[i].qji=0;
//形成节点导纳阵
G[line[i].i][line[i].j]=-line[i].r/((line[i].r*line[i].r+line[i].x*line[i].x));
G[line[i].j][line[i].i]=G[line[i].i][line[i].j];
B[line[i].i][line[i].j]=line[i].x/((line[i].r*line[i].r+line[i].x*line[i].x));
B[line[i].j][line[i].i]=B[line[i].i][line[i].j];
G[line[i].i][line[i].i]+=line[i].r/((line[i].r*line[i].r+line[i].x*line[i].x));
G[line[i].i][line[i].i]+=line[i].g;
B[line[i].i][line[i].i]+=(-line[i].x)/((line[i].r*line[i].r+line[i].x*line[i].x));
B[line[i].i][line[i].i]+=line[i].b;
G[line[i].j][line[i].j]+=line[i].r/(line[i].r*line[i].r+line[i].x*line[i].x);
G[line[i].j][line[i].j]+=line[i].g;
B[line[i].j][line[i].j]+=(-line[i].x)/(line[i].r*line[i].r+line[i].x*line[i].x);
B[line[i].j][line[i].j]+=line[i].b;
}
printf("\n请输入总节点数:");scanf("%d",npoint);
for(i=1;i=npoint;i++)
{
printf("请输入节点%d的数据\n",i);
printf("选择节点类型(1-平衡节点,2-PQ节点,3-PV节点):");
scanf("%d",(pointpri[i].sign));
if(pointpri[i].sign==1)
{
printf("请输入该节点的电压幅值和相位角(用逗号隔开):");
scanf("%f,%f",(pointpri[i].voltage),(pointpri[i].cta));
pointpri[i].p=0;pointpri[i].q=0;
nbal++;
}
else if(pointpri[i].sign==2)
{
printf("请输入该节点的注入有功和无功(用逗号隔开):");
scanf("%f,%f",(pointpri[i].p),(pointpri[i].q));
pointpri[i].voltage=1;pointpri[i].cta=0;
npq++;
}
else if(pointpri[i].sign==3)
{
printf("请输入该节点的有功注入和电压幅值(用逗号隔开):");
scanf("%f,%f",(pointpri[i].p),(pointpri[i].voltage));
pointpri[i].q=0;pointpri[i].cta=0;
npv++;
}
else
{
printf("选择无效,请重新输入\n");
i--;
}
if(nbal1)
{
printf("出现多个平衡节点,请重新输入\n");
i--;nbal--;
}
}
}
void view()
{
int i,j;
for(i=1;i=npoint;i++)
{
for(j=1;j=npoint;j++)
printf("%f+j%f ,",G[i][j],B[i][j]);
printf("\n");
}
for(i=1;i=npoint;i++)
{
printf("节点%d参数:",i);
if(pointpri[i].sign==1)printf("平衡节点 v=%f,theta=%f\n",pointpri[i].voltage,pointpri[i].cta);
else if(pointpri[i].sign==2)printf("PQ节点 P=%f,Q=%f\n",pointpri[i].p,pointpri[i].q);
else printf("PV节点 P=%f,v=%f\n",pointpri[i].p,pointpri[i].voltage);
}
}
void modify()
{
}
void run()
{
int i,j,k,m,n;
float E,e;//E为题目所要求的迭代精度,由键盘输入,e为每一步迭代的实际计算精度,用以和E比较
int numPQ[N];//PQ节点号
int numPV[N];//PV节点号
int numBAL; //平衡节点号
int np_th[N],nq_v[N];//迭代方程中的有功-电压相位角方程、无功-电压幅值方程分别对应的节点号
float B1_fac[N][N]={0},B2_fac[N][N]={0},a;
j=1;k=1;m=1;n=1;
for(i=1;i=npoint;i++){B1_fac[i][i]=1;B2_fac[i][i]=1;}
for(i=1;i=npoint;i++)
{//将数据存入一个新的结构体数组,保证在迭代过程中即使数据有改变,仍能保留住原始数据
pointaft[i].sign=pointpri[i].sign;
pointaft[i].voltage=pointpri[i].voltage;
pointaft[i].cta=pointpri[i].cta;
pointaft[i].p=pointpri[i].p;
pointaft[i].q=pointpri[i].q;
//统计PV节点和PQ节点序号,并存入临时数组
if(pointpri[i].sign==1)numBAL=i;
else if(pointpri[i].sign==2){numPQ[j++]=i;np_th[m++]=i;nq_v[n++]=i;}
else {numPV[k++]=i;np_th[m++]=i;}
}
//形成有功-相位角的迭代方程系数矩阵B1
for(i=1;i=npoint-1;i++)
for(j=1;j=npoint-1;j++)
B1[i][j]=B[np_th[i]][np_th[j]];
//形成因子表
for(k=1;k=npoint-1;k++)
for(i=k+1;i=npoint-1;i++)
{ a=B1[i][k];
for(j=1;j=npoint-1;j++)
{
B1[i][j]-=a*B1[k][j]/B1[k][k];
B1_fac[i][j]-=a*B1_fac[k][j]/B1[k][k];
}
}
//形成无功-电压幅值的迭代方程的系数矩阵B2
for(i=1;i=npq;i++)
for(j=1;j=npq;j++)
B2[i][j]=B[nq_v[i]][nq_v[j]];
//形成因子表
for(k=1;k=npq;k++)
for(i=k+1;i=npq;i++)
{ a=B2[i][k];
for(j=1;j=npq;j++)
{
B2[i][j]-=a*B2[k][j]/B2[k][k];
B2_fac[i][j]-=a*B2_fac[k][j]/B2[k][k];
}
}
printf("请输入迭代精度:");scanf("%f",E);
k=-1;
do
{
e=0;
k++;
//有功迭代
//首先计算有功不平衡量
for(i=1;i=npoint-1;i++)
{
detp[i]=pointaft[np_th[i]].p/pointaft[np_th[i]].voltage;
for(j=1;j=npoint;j++)
detp[i]-=pointaft[j].voltage*(G[np_th[i]][j]*cos(pointaft[np_th[i]].cta-pointaft[j].cta)+B[np_th[i]][j]*sin(pointaft[np_th[i]].cta-pointaft[j].cta));
printf("%f,",detp[i]);
}
for(i=1;i=npoint-1;i++)
{
detct[i]=0;
for(j=1;j=npoint-1;j++)
detct[i]+=B1_fac[i][j]*detp[j];
}
for(i=npoint-1;i=1;i--)
{
for(j=npoint-1;ji;j--)
detct[i]-=B1[i][j]*detct[j];
detct[i]/=B1[i][i];
}
for(i=1;i=npoint-1;i++)detct[i]/=(-pointaft[np_th[i]].voltage);//求出相位角的修正量
for(i=1;i=npoint-1;i++){pointaft[np_th[i]].cta+=detct[i];e+=detct[i]*detct[i];}//修正角度,并计算精度
printf("第%d次迭代后的相位角:",k);
for(i=1;i=npoint;i++)printf("%f,",pointaft[i].cta);
printf("\n");
//无功迭代
//首先计算无功不平衡量
for(i=1;i=npq;i++)
{
detq[i]=pointaft[nq_v[i]].q/pointaft[nq_v[i]].voltage;
for(j=1;j=npoint;j++)
detq[i]-=pointaft[j].voltage*(G[nq_v[i]][j]*sin(pointaft[nq_v[i]].cta-pointaft[j].cta)-B[nq_v[i]][j]*cos(pointaft[nq_v[i]].cta-pointaft[j].cta));
printf("%f,",detq[i]);
}
for(i=1;i=npq;i++)
{
detv[i]=0;
for(j=1;j=npq;j++)
detv[i]+=B2_fac[i][j]*detq[j];
}
for(i=npq;i=1;i--)
{
for(j=npq;ji;j--)
detv[i]-=B2[i][j]*detv[j];
detv[i]/=B2[i][i];
}
for(i=1;i=npq;i++)detv[i]/=(-1);//求出电压的修正量
for(i=1;i=npq;i++){pointaft[nq_v[i]].voltage+=detv[i];e+=detv[i]*detv[i];}//修正电压
printf("第%d次迭代后的电压幅值为:",k);
for(i=1;i=npoint;i++)printf("%f,",pointaft[i].voltage);
printf("\n");
}while(eEk100);
//输出结果
printf("在满足精度E=%f的情况下,计算结果为:",E);
for(i=1;i=npoint;i++)
printf("v%d=%f,cta%d=%f\n",i,pointaft[i].voltage,i,pointaft[i].cta);
//计算平衡节点注入功率
pointaft[numBAL].p=0;pointaft[numBAL].q=0;
for(i=1;i=npoint;i++)
{
pointaft[numBAL].p+=pointaft[i].voltage*(G[numBAL][i]*cos(pointaft[numBAL].cta-pointaft[i].cta)+B[numBAL][i]*sin(pointaft[numBAL].cta-pointaft[i].cta));
pointaft[numBAL].q+=pointaft[i].voltage*(G[numBAL][i]*sin(pointaft[numBAL].cta-pointaft[i].cta)-B[numBAL][i]*cos(pointaft[numBAL].cta-pointaft[i].cta));
}
pointaft[numBAL].p*=pointaft[numBAL].voltage;
pointaft[numBAL].q*=pointaft[numBAL].voltage;
printf("平衡节点注入功率为:%f+j%f\n",pointaft[numBAL].p,pointaft[numBAL].q);
//计算各支路潮流
for(k=1;k=nline;k++)
{
float aa,bb,cc;
i=line[k].i;j=line[k].j;
aa=pointaft[i].voltage*cos(pointaft[i].cta)-pointaft[j].voltage*cos(pointaft[j].cta);
bb=pointaft[i].voltage*sin(pointaft[i].cta)-pointaft[j].voltage*sin(pointaft[j].cta);
cc=line[k].r*line[k].r+line[k].x*line[k].x;
line[k].pij=(line[k].r*(pointaft[i].voltage*cos(pointaft[i].cta)*aa+pointaft[i].voltage*sin(pointaft[i].cta*bb))-line[k].x*(pointaft[i].voltage*sin(pointaft[i].cta*aa-pointaft[i].voltage*cos(pointaft[i].cta)*bb)))/cc;
line[k].qij=(line[k].x*(pointaft[i].voltage*cos(pointaft[i].cta)*aa+pointaft[i].voltage*sin(pointaft[i].cta*bb))+line[k].r*(pointaft[i].voltage*sin(pointaft[i].cta*aa-pointaft[i].voltage*cos(pointaft[i].cta)*bb)))/cc;
line[k].pji=(line[k].r*(pointaft[j].voltage*cos(pointaft[j].cta)*(-aa)+pointaft[j].voltage*sin(pointaft[j].cta*(-bb)))-line[k].x*(pointaft[j].voltage*sin(pointaft[j].cta*(-aa)-pointaft[j].voltage*cos(pointaft[j].cta)*(-bb))))/cc;
line[k].qji=(line[k].x*(pointaft[j].voltage*cos(pointaft[j].cta)*(-aa)+pointaft[j].voltage*sin(pointaft[j].cta*(-bb)))+line[k].r*(pointaft[j].voltage*sin(pointaft[j].cta*(-aa)-pointaft[j].voltage*cos(pointaft[j].cta)*(-bb))))/cc;
printf("由节点%d流向节点%d的功率为:%f+j%f\n",i,j,line[k].pij,line[k].qij);
printf("由节点%d流向节点%d的功率为:%f+j%f\n",j,i,line[k].pji,line[k].qji);
}
}
main ()
{
int a;
printf("\n");
printf("***************电力系统潮流计算**************\n");
printf("\n");
printf(" 设计者: 李 冰 设计:2014 年 5 月\n");
printf("\n");
printf("\n");
printf(" 注:各参数用标么值表示,角度均以弧度制表示\n");
printf("\n");
printf("\n");
do
{
printf("请选择:1-建立电网;2-查看数据;3-修改数据;4-运行计算;0-退出程序:");//由键盘输入选项
scanf("%d",a);
if (a==0) //如果选中0,就表示退出程序
{
exit(0);
}
else if (a==1) //如果选中1,就跳去执行建立电网,即creat处
{
creat();
}
else if (a==2) //如果选中2,就跳去执行查看数据,即view处
{
view();
}
else if (a==3) //如果选中1,就跳去执行修改数据,即modify处
{
modify();
}
else if (a==4) //如果选中1,就跳去执行运行计算,即run处
{
run();
}
else //如果以上都不是,则表示无效输入,重新要求输入选择
{
printf("\n无效输入,请重新选择\n");
}
}while(1);
}
不知道能不能帮到你