用c程序实现的算法

机器实验算法总结

程序设计算法与实现

一、数据结构

给定一段数组a[20]（为了便于以后的插入移动操作，这里我们给这个数组分配的空间大些）

初始化10个数，即a[0]~a[9]。原数组为a[0]~a[19]大小，我们只用了a[0]~a[9]，因此称这段区域为有效区域。当然删除有效区域的数据元素，有效区要变小。最后只输出有效区。

图示一

问题一：如何删除数组里一段连续的数据？

先定义数组a[10]，用于存放10个数字，由于是向前覆盖删除数据，所以这里可以不需要多定义空间。有效区域是a[0]~a[9]

上面是为这个数组分配了空间，下面需要初始化数组….

初始化数组有两种方法，一种是定义后直接赋值初始化，一种是用户输入用scanf接收后赋值初始化，第一种明显没有第二种灵活，建议使用第二种做法。代码片段如下

int a[10],i;//这里定义了数组a[10]和变量i for(i=0;i<10;i++) scanf(“%d”,&a[i]); 代码一 上述代码用一个for循环让用户自己输入数值进行a[10]的初始化。为检验初始化的结果，我们可以将a[0]~a[9]的初始化数值也用for循环配合printf输出出来，代码如下：

for(i=0;i<10;i++) printf(“%-5d”,a[i]);/*这里%-5d为格式化输出控制符，%d表示输出整型（int）]为先输出5个列空格再输出数字，%-5d为先输出数字后输出5个列空格。*/ 代码二

为了解决删除数组中连续的元素这个问题，可以先考虑如何删除数组里某个元素。我们知道，若要删除某个元素，

其做法就是让这个元素后面的元素依次向前移动，这样这个待删除的元素就会被他后面的元素所替代，这样依次下去，达到目的。

王晔2007年10月20日编写

C语言与程序的算法

第1章 C语言与程序的算法

总体要求：

了解C程序的一般结构及特点。

了解C程序的编辑、调试与运行等上机操作过程。

理解算法的概念及其特性，了解算法的描述方法。

相关知识点：

了解计算机软件、计算机语言及分类的知识。

熟悉Windows系统基础知识及操作。

学习重点：

C程序的结构、特点以及算法的概念。

C语言是20世纪70年代初出现的、国际上广泛流行的计算机高级语言，它既可以用来编写底层系统软件，也可以用来开发应用软件。目前，C语言在我国已经成为学习程序设计的入门语言。本章将首先向大家展示几个简单的C程序，并介绍C程序的特点，然后在Visual C++ 6.0中演示C程序的编辑、编译、调试和运行的操作过程，最后介绍程序算法的概念及其描述方法。

1.1 简单的C程序及特点

1.1.1 几个简单的C程序

【实例1-1】 显示Hello, World!

#include "stdio.h"

main()

{

} printf("Hello,World!\n");

本程序的作用是输出以下一行信息：

在该程序中，#include用来引用C语言的库函数文件。stdio.h是标准输入输出库函数文件。main

#include #include\

#define PMAX 100 #define PMIN -100 #define DMAX 100

#define DMIN -100 #define FMAX 100

int PFF[4]={0,12,24,48}; /*语言值的满幅值*/

int DFF[4]={0,16,32,64}; /*输入量D语言值特征点*/ int UFF[7]={0,15,30,45,60,75,90}; /*输出量U语言值特征点*/

/*采用了调整因子的规则表,大误差时偏重误差,小误差时偏重误差变化*/ /*a0=0.3,a1=0.55,a2=0.74,a3=0.89 */

int rule[7][7]={ //误差变化率 -3,-2,-1, 0, 1, 2, 3 // 误差 {-6,-6,-6,-5,-5,-5,-4,}, // -3 {-5,-4,-4,-3,-2,-2,-1,}, // -2

用C程序实现对二元关系性质的判定

ＩＴ技术

数字技术与应用

用Ｃ程序实现对二元关系性质的判定

岳晓红

（陇东学院信息工程学院　　甘肃庆阳　　７４５０００）

摘　要：本文利用Ｃ语言程序实现了对离散数学中二元关系的５种性质的判断。在计算机相关专业的教学中可以培养学生的理论和计算机操作相结合的学习能力。

关键词：Ｃ语言　　离散数学　　二元关系中图分类号：Ｏ１４１．１２文献标识码：Ａ文章编号：１００７－９４１６（２０１１）０２－００７１－０２

Ｔｏ　Ｊｕｄｇｅ　ｔｈｅ　Ｐｒｏｐｅｒｔｙ　ｏｆ　Ｂｉｎａｒｙ　Ｒｅａｌｔｉｏｎｓ　ｂｙ　Ｃ　ｐｒｏｇｒａｍ

Ｙｕｅ　Ｘｉａｏｈｏｎｇ

（Ｓｃｈｏｏｌ　ｏｆ　Ｉｎｆｏｒｍａｔｉｏｎ　Ｅｎｇｉｎｅｅｒｉｎｇ　　Ｑｉｎｇｙａｎｇ　ｇａｎｓｕ　　７４５０００）

Ａｂｓｔｒａｃｔ：Ｔｈｅ　ｐａｐｅｒ　ｐｒｅｓｅｎｔｓ　ｓｏｍｅ　ｗａｙｓ　ｔｏ　ｊｕｄｇｅ　ｔｈｅ　５　ｐｒｏｐｅｒｔｉｅｓ　ｏｆ　ｂｉｎａｒｙ　ｒｅｌａｔｉｏｎｓ　ｉｎ　ｔｈｅ　ｄｉｓｃｒｅｔｅ　ｍａｔｈｅｍａｔｉｃｓ，　ａｎｄ　ｔｈｅｙ　ｃａｎ　ｂｅ　ａｐｐｌｉｅｄ　ｔｏ　ｃｕｌｔｉｖａｔｅ　ｔｈｅ　ｓｔｕｄｅｎｔｓ’ａｂｉｌｉｔｉｅｓ　ｂｙ　ｉｎｔｅｇｒａｔｉｎｇ　ｔｈ

16点的FFT算法VERILOG程序实现

问题的提出 解决问题的思路与方法

基2时间抽取FFT算法基2时间抽取FFT算法的计算复杂度

基2时间抽取FFT算法流图规律基2频率抽取FFT算法

FFT算法的实际应用

问题的提出4点序列{2,3,3,2} DFT的计算复杂度

X [m] k 0

N 1

km x[k ]WN ,

m 0,1, N 1

0 0 0 0 X [0] 2WN 3WN 3WN 2WN 10 0 1 2 3 X [1] 2WN 3WN 3WN 2WN 1 j 0 2 4 6 X [2] 2WN 3WN 3WN 2WN 0 0 3 6 9 X [3] 2WN 3WN 3WN 2WN 1 j

如 何 提 高 DFT 的 运 算 效 率

复数加法 N(N-1)

复数乘法 N 2

?

解决问题的思路1. 将长序列DFT分解为短序列的DFTkm 2. 利用旋转因子 WN 的周期性、对称性、可约性。

旋转因子1)周期性

的性质

km WN

(k N )m k ( m N ) km WN WN WN

2) 对称性WNmk N 2

mk W N

km

一 需求分析

1．本程序演示的是用简单遗传算法随机一个种群，然后根据所给的交叉率，变异率，世代数计算最大适应度所在的代数

2．演示程序以用户和计算机的对话方式执行，即在计算机终端上显示“提示信息”之后，由用户在键盘上输入演示程序中规定的命令；相应的输入数据和运算结果显示在其后。 3．测试数据

输入初始变量后用y=100*(x1*x1-x2)*(x1*x2-x2)+(1-x1)*(1-x1)其中-2.048<=x1,x2<=2.048作适应度函数求最大适应度即为函数的最大值

二 概要设计 1．程序流程图

开始 Gen=0 编码 随机产生M个初始个体 满足终止条件? Y 输出结果 终止 N 计算群体中各个体适应度 从左至右依次执行遗传算子 pm j = 0 选择个体变异点 执行变异 变异后添入新群体中 pc j = 0 根据适应度选择复制个体 执行复制 复制的个体添入新群体中 j = j+1 j = 0 选择两个交叉个体 执行交叉 交叉后添入新群体中 j = j+2 j = j+1 N Y N Y Gen=Gen+1 N Y

2．类型定义

int popsize; //种群大小

int maxgeneration; //最大世代数 double pc; //交叉率

遗传算法的c语言程序

一 需求分析

1．本程序演示的是用简单遗传算法随机一个种群，然后根据所给的交叉率，变异率，世代数计算最大适应度所在的代数

2．演示程序以用户和计算机的对话方式执行，即在计算机终端上显示“提示信息”之后，由用户在键盘上输入演示程序中规定的命令；相应的输入数据和运算结果显示在其后。 3．测试数据

输入初始变量后用y=100*(x1*x1-x2)*(x1*x2-x2)+(1-x1)*(1-x1)其中-2.048<=x1,x2<=2.048作适应度函数求最大适应度即为函数的最大值

二 概要设计 1．程序流程图

2．类型定义

int popsize; //种群大小

int maxgeneration; //最大世代数 double pc; //交叉率 double pm; //变异率 struct individual

遗传算法的c语言程序

{

char chrom[chromlength+1]; double value;

double fitness; //适应度 };

int generation; //世代数 int best_index;

int worst_

#include\int E0=0; int E1=0; int E2=0; int Error0=0; int Error1=0; int Ppid=0;

int Pctr(int ch1,int ch2,int Kp,int Ki,int Kd) {

E0=ch1-ch2;//增量计算当前E0

Error0=E0-E1;//当前差值Eo是差值

Error1=E1-E2;//上一次差值上一次值也再上一次值得差值 Ppid=(int)(Ki*Error0+Kp*E0+Kd*(Error0-Error1)); E2=E1;//

E1=E0;//存储误差，用于下次计算

return Ppid;//Return }

/*====================================================================================================

增量式PID计算部分 //增量式PID算法核心部分！P控制误差D

==============

#include\int E0=0; int E1=0; int E2=0; int Error0=0; int Error1=0; int Ppid=0;

int Pctr(int ch1,int ch2,int Kp,int Ki,int Kd) {

E0=ch1-ch2;//增量计算当前E0

Error0=E0-E1;//当前差值Eo是差值

Error1=E1-E2;//上一次差值上一次值也再上一次值得差值 Ppid=(int)(Ki*Error0+Kp*E0+Kd*(Error0-Error1)); E2=E1;//

E1=E0;//存储误差，用于下次计算

return Ppid;//Return }

/*====================================================================================================

增量式PID计算部分 //增量式PID算法核心部分！P控制误差D

==============

分步傅里叶算法的MATLAB程序实现举例

模型：

线性部分： 非线性部分：2

i U1 U z

2 x

2 nU 0

2

2

n U x,z d z

n x

2其中：

d z exp z 或d z exp z

2

U z

=

i U2 x

2

两边同时对x变量作傅里叶变换

~

U z

=

i2

i

2

~

U

两边积分

x h21

~

U=x hix

~

U

x

2

i

2

dz

即

~

ln

U x h/2 ~

=

iU x

2

i

2

*

h2

最后有

~

~

U x h/2 =U x exp i

2h i * 2

2

再对x变量作作傅里叶逆变换

U x h/2 =FFT-1

FFT i2h U x exp 2 i *

2

U z

inU

两边积分

x h1 hx

U

=

xx

indz

当h 0时

ln

U x h U x

inh

最后有

U x h U x exp inh

2

折射率部分： n U x,

两边同时对x变量作傅里叶变换

n=FFT

~~

2

nd 2

x

2

d z * i n

2

~

2

FFT

n= 2

1 d z

再对x变量作作傅里叶逆变换

FFT U2

-1

n=FFT 2

1 d z