DSP总结之C语言环境

 

资料来源于:重剑无锋的BLOG 原文作者为:wxjia2004

1.C语言软件开发过程



   软件开发过程涉及编译器(compiler),汇编器(assembler),连接器(linker),归档器(archiver),建库器(library-build utility),运行支持库(run time support library),HEX转换器(hex conversion utility),交叉引用列表器(cross reference lister),绝对列表器(absolute lister)等。其大都设置既可通过命令,也可通过ccs的project\build options设置。

   1>.编译器:对C语言源代码进行编译,产生汇编语言源代码。

      一步编译:cl2000 –v28 [-options] filenames [object files] [-z [link_options]]

      分布编译:文法分析:  生成file.if的中间文件 ac2000 file.c

                代码优化:  生成file.opt的文件     opt2000 file.if

                代码生成:  生成file.asm的汇编文件 cg2000 file.opt

                汇编:      生成file.obj目标文件   asm2000 file.asm

   2>.汇编器:把汇编语言源文件翻译成机器语言目标文件,机器语言格式为公用目标格式(COFF)。具体命令格式如下:

      asm2000 Version [input file [object file [listing file]]] [options]

   3>.连接器:把多个目标文件组合成单个可执行目标模块。它一边创建可执行模块,一边完成重定位以及决定外部参考。连接器的输入是可重定位的目标文件和目标库文件.具体命令格式如下:

      汇编语言:lnk2000 [options] filename1,filename2……[filenamen]

      C语言:lnk2000 {-c|-cr} filenames [options] [-o name.out] [lnk.cmd] [-l libname,lib]

         注:        -c             运行时自动初始化变量

                     -cr            载入程序时自动初始化变量

                     options:       详见手册。

                     -o name.out    生成的输出文件名

                     lnk.cmd        连接的cmd文件名

                     -l libname.lib 连接的运行支持库名

         例:lnk2000 –c prog1 prog2 prog3 –o prog.out rts2800.lib

4>.归档器(Archiver)

      归档器Archiver可以用来对文档(Archive)或者库(library)中的文件进行分离和合并。这些文档或库可以是源文件库,也可以时目标文件库。归档器可以对库进行新建,添加,删除,替换,提取等操作,具体命令格式如下:

      ar2000 [-]command [options] libname [filename0,filename1,……filenamen]

1>     此命令针对TMS320F2812而言,其他类DSP命令以及格式略有不同,可详查TI的文档。此命令所在目录在\CCS_v3.1\C2000\cgtools\bin

2>     command:

@                使用CMD文件

a (add)         向指定文档中添加指定文件

d (delete)      删除指定文档中的指定文件

r (replace)    替换指定文档中的指定文件

t (table)        列出指定文档中的文件

x (extract)    提取指定文档中的指定文件

3>     options:

q(quiet)       屏蔽状态信息

s(symbol)    列出库中定义的全局符号(对命令a,r,d无效)

u                替换文件时同步更新修改日期

v(verbose)   提供详细的描叙

4>     libname 指定的文档名

5>     filename 文档中指定的文件名

   5>.建库器(library-build utility): 建立满足你自己要求的运行支持库.CCS中提供有rts库文件(\CCS_v3.1\C2000\cgtools\lib),并提供了对应的源程序文件rts.src(\CCS_v3.1\C2000\cgtools\lib)。如:rts2800.lib:C/C++运行支持库;rts2800_ml.lib C/C++大内存模式运行支持库;用户可以利用归档器和建库器对rts.src中的文件进行提取,修改等操作,并重新生成自己需要的库文件。

命令的具体格式如下:

      mk2000 –v28 [options] src_arch1 [-l objlib1] src_arch1 [-l objlib1]……

   6>.运行支持库(run time support library): 它包括C编译器所支持的ANSI标准运行支持函数、编译器公用程序函数、浮点运算函数和C编译器支持的I/O函数. CCS中提供有rts库文件(\CCS_v3.1\C2000\cgtools\lib),并提供了对应的源程序文件rts.src(\CCS_v3.1\C2000\cgtools\lib)。如:rts2800.lib:C/C++运行支持库;rts2800_ml.lib C/C++大内存模式运行支持库.库文件包括标准的C/C++运行支持库函数,浮点运算程序,系统启动程序_c_int00.

   7>.HEX转换器(hex conversion utility): 它把COFF目标文件转换成TI-Tagged、ASCII-hex、 Intel、 Motorola-S、或 Tektronix 等目标格式,可以把转换好的文件下载到EPROM编程器中.

      hex2000 [-options] filename

          -a    ASCII-HEX

          -i    Intel

          -t    Ti-Tagged

          -m    Motorola-S

          -x    Tektronix

   8>.交叉引用列表器(cross reference lister): 它用连接的目标文件产生参照列表文件,可显示符号及其定义,以及符号所在的源文件.

      xref2000 [options] [input filename [output filename]]

          options:  -lxx        每页显示xx行(十进制数)

                    -q(quiet)   不显示提示消息

          input filename        obj文件

          output filename       生成的xrf文件

   9>.绝对列表器(absolute lister): 它输入连接后的目标文件,输出.abs文件,通过汇编.abs文件可产生含有绝对地址的列表文件。如果没有绝对列表器,这些操作将需要冗长乏味的手工操作才能完成。

      abs2000 [-options] input file

          options:  -e      改变缺省的文件后缀名

                        -ea[.]asmext    asm->asmtxt

                        -ec[.]cext      c->ctxt

                        -eh[.]hext      h->htxt

                        -ep[.]cppext    cpp/cc/cxx ->ptxt

                    -fs     指定输出文件目录。如:–fs C:\ABSDIR

                    -q      不显示提示消息

          input file        .out文件

      例如:abs2000 –ea s –ec csr –eh hsr hello.out生成文件hello.s(hello.asm),hello.csr(hello.c),hello.hsr(hello.h).

2.CCS概述

      DSP芯片开发工具可分为两大类:

        代码生成工具:
         代码调试工具:C/汇编语言源码调试器,初学者工具DSK,软件模拟器(Simulator),评价模块EVM,软件开发系统SWDS,仿真器XDS.
   3.C编译器的优化:

      1>C编译器优化分为2类:

       C语言通用优化:简化表达式,数据流优化,删除公共子表达式和冗余分配,优化跳转,简化控制流,优化与循环有关的变量,将循环体内计算值不变的表达式移至循环体前,运行支持库函数的行内扩展。

        根据DSP芯片的特定优化:高效地使用寄存器,自动增量寄存器寻址方式,使用块重复,使用并行指令,使用延迟跳转。

      2>CCS的优化选项:



   4.GEL语言的使用:

      GEL(General Extension Language通用扩展语言)是一种类似于C语言的一种解释性语言,它可以创建GEL函数,以扩展CCS用途。

      GEL是C语言的一个子集,然而它不能声明主机变量,所有的变量必须在DSP程序中定义,存在于仿真/实际目标板中,唯一不在目标板上定义的标识符是GEL函数及其参数。

      GEL函数可在任何能键入C表达式的地方调用,既可以在任何可键入C表达式的对话框中调用,也可以在其他GEL函数中调用。但不支持递归。

可以将常用的GEL函数添加到CCS的GEL菜单下,此时需要使用menuitem关键词在GEL菜单下创建一个新的下拉菜单列表(一级菜单),再使用hotmenu,dialog和slider在该菜单项中添加新的菜单项(二级菜单)。

在CCS启动是自动执行GEL函数。SETUP CCS时设置环境时设置了自动执行GEL函数。自动运行StartUp()函数。这样要求每个工程建立时都载入GEL文件。

CCS提供了一系列嵌入GEL的函数。使用这些函数,用户可以控制仿真/实际目标板的状态,访问存储器,并在输出窗口中显示结果。

   5.Probe Point的应用

      Probe Point只是暂时中断程序的运行,更新与之相连接的窗口,然后自动运行以后的程序。与之相连接的窗口可在Probe Point设置中设置好,然后replace。

      可使用Probe Point从PC文件中导入数据,也可存储数据到PC文件中,对PC文件的格式只能使用两种:COFF程序文件(.out);CCS数据文件(.dat)。

      可利用Probe Point和Break Point配合显示图形和动画。显示图形使用View->Graph->Time/Frequency,设置好跳出的对话框,可看到一块内存的数据,可使用Probe Point将PC上的数据传给目标板,然后继续运行程序。然后可以创建断点,使图形窗口自动更新,使用Animate命令,使到达断点后更新窗口后程序自动继续执行。

1.  用户库的创建:

1>    DOS命令:使用归档器ar2000命令和汇编编译连接等命令进行手工创建。

2>    CCS直接创建,在CCS中建立工程对话框的时候有一个可选择建立lib文件的选项。可直接把工程中的asm,c等文件直接整合生成一个库文件(不包括CMD文件)。

DSP算法大全C语言版本 完整版,共407页,审阅过的。包含多种数字信号产生、处理、分析方式,并附参考代码。第六章FIR数宇滤波器的设计………………… 227 §6.1窗函数方法……………… .227 §6.2频域最小误差平方设计…… “·自“238 §6.3切比雪夫逼近方法……………………………………………242 第三篇随机数字信号处理 第一章经典谱佔计……………………………………….264 的周期图方法………………·…264 12功率谱估计的相关方法 271 第二章现代谱估计……… 280 §2.1求解一般托布利兹方程组的莱文森算法……………0 82.2求解对称正定方程组的乔里斯基算法 §2.3求解尤利沃克方程的莱文森德宾算法…………灬….28 §24计算ARMA横型的功率谱密度18 §2.5尤利沃克谱佔计算法…………………*…………………22 §2.6协方差谱估计算法 ……ts29 §27Burg谱估计算法 ◆鲁b+吾·合品品‘山点亠4+·日叶中·中‘甲争早导 §2.8最大似然谱估计算法…… 308 第三章时频分析 甲甲甲手曾鲁卧鲁哲雪 §3.1堆格纳( wigner)分布 中中中“由节“昏音山曲画 32离散小波变换 委甲■即 318 第四章随机信号的数字滤波……… §41维纳( Wiener)数字滤波……… …………量·自330 §4-2卡尔曼( Kalman)数字滤波… 会血中自4B44品西4垂4+中如甲吾卧d古 §4-3最小均方(LMS)自适应数字滤波…… ●·中自·自·中中中平看 °·341 §4.4归一化LMS自适应数字滤波……34 §45递推最小二乘(RLs)自适应数字滤波 348 第四篇数字图像处理 第一章图像基本运算 ……·352 §1,1图像读取、存储与显示…… 81.2图像旋转……………………………………………………………366 1.3图像灰度级直方图的计算……………………………………*……3 §1.4图像二值化的固定阀值法… ma·也d■血dp §1.5图像二值化的自适应阀值法 导b血 第二章图像增颯… §2.1图像直方图均衡…………………… w"376 §2.2中值滤波……38 §23图像锐化…… 鲁平t自d “*………382 §2.4图像平滑 P····383 第三章图像边缘检测……………………….356 831 Roberts算子边缘检测…46 §32拉普拉斯算子边缘检测……… “………“敌…388 83.3 Sobel算子边缘检测……; 83.4 Robinson算子边缘检测 小392 §35 Kirsch算子边缘检测……………………………………………394 §3.6 Prewitt算子边缘检测 争昏平辛辛平中萨 396 第四章图像细化…… 喜即香看d画命合b分bb画品目如画如bL晶 品品4甲。自·。·中 399 §4.1 Hilditch细化算法… ●·命··“““…“399 §4.2 Pavlidis细化算法……404 §43 Rosenfeld细化算法………………… 0 第五篇人工神经网络 第一章神经网络模型 t……416 §1.1多层感知器神经网络 ·416 §1.2离散 Hopfield神经网络…………“………………………425 §13连续 Hopfield神经网络………………49434 §I4 Tank-Hopfield线性规划神经网络 ……437 参考文献 44命↓◆命啡4每◆普““女4古“4b中d●·4·面···4·= ,··442 第一篇常用数字信号的产生 第一章数字信号的产生 §1.1均匀分布的随机数 -、功能 产生(a,b)区间上均匀分布的随机数。 方法简介 均匀分布的概率密度函数为 ≤x≤b fCx) 其它 通常用U(b表示均匀分布的均值为“士方差为2 产生均匀分布随机数的方法如下 首先,由给定的初值x0,用混合同余法 T-1+c)(mod M) y; /M 产生(0,1)区间上的随机数y。其中a=2045c=1M=2然后,通过变换z=a+ (b…a)y产生(a,b)区间上的随杋数x 三、使用说明 1.子函数语句 double uniform (a, b, seed) 2形参说明 双精度实型变量。给定区间的下限 b—双精度实型变量。给定区间的上限。 seed——长整型指针变量。*seed为随机数的种子。 四、子函效程序(文件名: uniform,c) double uniform(a,bseed) long int seed double t; 2045爷( seed-#seed-(操Seed/1948576)*1048576; t=(“seed)/1048576.0; t=a+(b-a)头t; return(t)i 五、例题 产生50个0到1之间均匀分布的随机数。 主函数程序(文件名: uniform,m) 杜 include〃 stdio.h includ iform ain( doable a, b,x; int 1,J; g Int s double uniform (double, double, long int *) a-0.0;b=10;s=13579; for(i-0;i<10;i-+) for(j=0:j<5++) =unite printf (". 7fM,x>; intf("\n")y 运行结果 0.48263550.98959450.72067070.77158260.8864250 0.73916340.58915140.81457810.81212620.7979975 0.9483t 0.39095970.51266860.40730760.9440937 0.67162610.47535710.10517980.09266470.4993505 0.13187120.47657490.59566690,13878150,8082657 C.90332890.30759240.02644250.07497020.3141527 0.44230840.5207319 89678960.93463230.3230572 0.65192320.18290330.03722860.13245580.8721647 0,5768661 0.6912775 0.66249660.80548000.2066078 0.51299000.0645466099776740.43443300.4154520 §1.2正态分布的随机数 功能 产生正态分布N(,a2)的随机数。 二,方法简介 正态分布的概率密度函数为 √2πa 通常用、(p,d2)表示。式中p是均值,a2是方差。正态分布也称为高斯分布 产生正态分布随机数的方法如下 设r1 r为(0,1)上n个相互独立的均匀分布的随机数,由于E(n)=1, 根据中心极限定理可知,当n充分大时 r 的分布近似于正态分布N(,1)。通常取n=12,此时有 r一6 最后,再通过变换y=H+x,便可得到均值为八方差为m2的正态分布随机数y 三、使用说明 1.子函数语句 double gauss(mean, sigma, seed) 2.形参说明 mean—双精度实型变量。正态分布的均值k sigma—-双精度实型变量。正态分布的均方差o see 长整型指针变量。餐seed为随机数的种子。 四、子函数程序(文件名: gauss,c) include Uniform,ci double gauss(mean, sigma, s) double mean, sigma; long int *si i int i; double x, y; double uniform() for(x=0,=0;<12i++) x+= uniform(0,0,1.0,) X=x一6.0: y〓mean十x并sgm丑 return (y)f 五、例题 产生50个均值为0方差为1的正态分布的随机数。 主函数程序(文件名 tgauss.m) #include"stdio. hm #include"gauss, c main() t int i,j s long int s# double x, mean, sigma; double gauss(double, double, long int *) mean=0.0; sigma=1.05=13579 for(i=0<10i++) for(j=0<5i十+) x=gauss(mean, sigma, &s) printf(".7",x); printi("\n") 运行结果 2.8997211-0,90885730,2041950-0.2572155-0.8516827 0.79969980,9866190.04313851.919498702543507 0.36892511.2145863 ,05370901.70509531.6925945 0.49287221.9956684.-0.59806631.29232980.1707630 0.5213604-0.40513420.8358479-0.54450801.6452045 0.5338917-0.8120403-0.3886852-0.25463680.4690113 0.4013348-0.1117687-0.9708830.650224713179646 0.53624150.74646191.3275318-0.40414241.8053455 0.8525982-0,24906731.68234440.945543304819355 1.1704273-0.172575002068348-1.9993710.8360157 §1.3指数分布的随机数 功能 产生指数分布的随机数。 方法简介 1.产生随机变量的迆变换法 定理设F(x)是任一连续的分布函数如果~U(0,1)且=F-1(a),那么η~F(x)。 证明由于t~U(0,1),则有 P(≤x)=P(F-(x)≤x)=P(x≤F(x))=F(x) 所以,7~F(x),定理证毕 此定理给出了从均匀分布随机数到给定分布F(x)的随机数的变换。根据该变换可产 生分布函数为F(x)的机效x,其算法可用下列两个步骤实现 (1)产生均匀分布的随机数g,即4~U(0,1)g(2)计算x=F-(t) 2.产生數分布随机的方法 指数分布的概率密度函数为 x≥0 f(x) 0 其它 其分布函效为 FCr)= ,其它 指数分布的均值为,方为P2 根播上述的逆交换法,产生指数分布随机数的方法为 (1)产生均匀分布的随机数M,即w~U(0,1);(2)计算x=-ln(v) 、使用说明 1.子函微语句 double exponent(beta, s) 2.彩参说明 beta——双精皮实型变量。指数分布的均值。 s—长整型指针变量。¥s为随机数的种子。 四、子函数程序(文件名: exponent.[) include "math. h" Include uniform. c uble exponent(bcta, double beta; long int *s i double u,x double uniform) u=uniform(0. 0,1-0,s)+ ta i logt return(x) 五、例题 产生50个均值为2、方差为4的指数分布的随机数。 主函数程序(文件名; exponent n) toinclude stdio. h f include exponentc int i,j; long int s; ouble x 4a; double exponent()+ bea=2.仍;s=t3579; for(i=0;i<10;i++) {for(j=0;j5;j+) i x=exponent(beta, &s printf(".7f", x> 运行结果 45698710,02092010,6551459051862310,2411175 0.6044725 1.0581442 0.4101700 0.4161992 0.451299 0.20000171,87830141.33625131.79637310.1150597 0,7961070 .4873781450416854.75753501.3888938
几个c程序,希望有用 include <math.h> /* In the VMS compiler, M_PI is not defined in math.h */ #ifdef vms #define M_PI 3.14159265358979323846 #endif #ifndef KR #define M_PI 3.14159265358979323846 #endif #define FALSE 0 #define TRUE 1 #define BIG 1e10 #define SMALL 1e-10 #define ORDER5 1e-5 #define ORDER4 1e-4 #define ABS(x) ((x) >= 0 ? (x) : -(x)) #define MIN(a,b) ((a) <= (b) ? (a) : (b)) #define MAX(a,b) ((a) >= (b) ? (a) : (b)) /* Table of constant values */ static long neg_i1 = -1; static long pos_i1 = 1; static long pos_i2 = 2; static long pos_i4 = 4; static long pos_i5 = 5; static float neg_f1 = -1.0; static float pos_f2 = 2.0; static double pos_d10 = 10.0; typedef struct { float r, i; } complex; typedef struct { double r, i; } doublecomplex; /* .................End of include file.................*/ /* SPBFCT 11/14/85 */ /* GENERATES (I1)!/(I1-I2)!=I1*(I1-1)*...*(I1-I2+1). */ /* NOTE: 0!=1 AND SPBFCT(I,I)=SPBFCT(I,I-1)=I!. */ #ifndef KR double spbfct(long *i1, long *i2) #else double spbfct(i1, i2) long *i1, *i2; #endif { /* Local variables */ long i; double ret_val; ret_val = 0.0; if (*i1 < 0 || *i2 < 0 || *i2 > *i1) { return(ret_val); } ret_val = 1.0; for (i = *i1 ; i >= (*i1 - *i2 + 1) ; --i) { ret_val *= (double) i; } return(ret_val); } /* spbfct */ /* SPBILN 11/13/85 */ /* CONVERTS ANALOG H(S) TO DIGITAL H(Z) VIA BILINEAR TRANSFORM */ /* ANALOG TRANSFER FUNCTION DIGITAL TRANSFER FUNCTION */ /* D(L)*S**L+.....+D(0) B(0)+......+B(L)*Z**-L */ /* H(S)=--------------------- H(Z)=---------------------- */ /* C(L)*S**L+.....+C(0) 1+.......+A(L)*Z**-L */ /* H(S) IS ASSUMED TO BE PRE-SCALED AND PRE-WARPED */ /* LN SPECIFIES THE LENGTH OF THE COEFFICIENT ARRAYS */ /* FILTER ORDER L IS COMPUTED INTERNALLY */ /* WORK IS AN INTERNAL ARRAY (2D) SIZED TO MATCH COEF ARRAYS */ /* IERROR=0 NO ERRORS DETECTED IN TRANSFORMATION */ /* 1 ALL ZERO TRANSFER FUNCTION */ /* 2 INVALID TRANSFER FUNCTION; Y(K) COEF=0 */ #ifndef KR void spbiln(float *d, float *c, long *ln, float *b, float *a, float *work, long *error) #else void spbiln(d, c, ln, b, a, work, error) long *ln, *error; float *d, *c, *b, *a, *work; #endif { /* Local variables */ long i, j, l, zero_func, work_dim1; float scale, tmp; double atmp; zero_func = TRUE; for (i = *ln ; i >= 0 && zero_func ; --i) { if (c[i] != 0.0 || d[i] != 0.0) { zero_func = FALSE; } } if ( zero_func ) { *error = 1; return; } work_dim1 = *ln + 1; l = i + 1; for (j = 0 ; j <= l ; ++j) { work[j * work_dim1] = 1.0; } tmp = 1.0; for (i = 1 ; i <= l ; ++i) { tmp = tmp * (float) (l - i + 1) / (float) i; work[i] = tmp; } for (i = 1 ; i <= l ; ++i) { for (j = 1 ; j <= l ; ++j) { work[i + j * work_dim1] = work[i + (j - 1) * work_dim1] - work[i - 1 + j * work_dim1] - work[i - 1 + (j - 1) * work_dim1]; } } for (i = l ; i >= 0 ; --i) { b[i] = 0.0; atmp = 0.0; for (j = 0 ; j <= l ; ++j) { b[i] += work[i + j * work_dim1] * d[j]; atmp += (double) work[i + j * work_dim1] * c[j]; } scale = (double) atmp; if (i != 0) { a[i - 1] = (double) atmp; } } if (scale == 0.0) { *error = 2; return; } b[0] /= scale; for (i = 1 ; i <= l ; ++i) { b[i] /= scale; a[i - 1] /= scale; } if (l < *ln) { for (i = l + 1 ; i <= *ln ; ++i) { b[i] = 0.0; a[i - 1] = 0.0; } } *error = 0; return; } /* spbiln */ /* SPBSSL 11/13/85 */ /* GENERATES ANALOG FILTER COEFFICIENTS FOR LTH ORDER */ /* NORMALIZED LOWPASS BESSEL FILTER */ /* COEFFICIENTS ARE RETURNED IN ARRAYS D AND C */ /* LN SPECIFIES ARRAY SIZE (LN>=L) */ /* WSCL CONTROLS FREQUENCY SCALING SUCH THAT RESPONSE AT 1 RAD/SEC */ /* IS EQUAL TO THAT OF UNSCALED H(S) AT WSCL RAD/SEC */ /* IERROR=0 NO ERRORS DETECTED */ /* 1 INVALID FILTER ORDER (L<=0 OR L>=LN) */ /* 2 INVALID SCALE PARAMETER (WSCL<=0) */ #ifndef KR void spbssl(long *l, float *wscl, long *ln, float *d, float *c, long *error) #else void spbssl(l, wscl, ln, d, c, error) long *l, *ln, *error; float *d, *c, *wscl; #endif { /* Builtin functions */ long pow_ii(); double pow_ri(); /* Local variables */ double spbfct(); long k, ll, tmp_int1, tmp_int2; if (*l <= 0 || *l > *ln) { *error = 1; return; } if (*wscl <= 0.0) { *error = 2; return; } for (ll = 0 ; ll <= *ln ; ++ll) { d[ll] = 0.0; c[ll] = 0.0; } for (k = 0 ; k <= *l ; ++k) { tmp_int1 = (*l * 2) - k; tmp_int2 = *l - k; c[k] = pow_ri(wscl, &k) * spbfct(&tmp_int1,&tmp_int1) / (pow_ii(&pos_i2, &tmp_int2) * spbfct(&k,&k) * spbfct(&tmp_int2,&tmp_int2)); } d[0] = c[0]; *error = 0; return; } /* spbssl */ /* SPBWCF 11/13/85 */ /* GENERATES KTH SECTION COEFFICIENTS FOR LTH ORDER NORMALIZED */ /* LOWPASS BUTTERWORTH FILTER */ /* SECOND ORDER SECTIONS: K<=(L+1)/2 */ /* ODD ORDER L: FINAL SECTION WILL CONTAIN 1ST ORDER POLE */ /* LN DEFINES COEFFICIENT ARRAY SIZE */ /* ANALOG COEFFICIENTS ARE RETURNED IN D AND C */ /* IERROR=0 NO ERRORS DETECTED */ /* 1 INVALID FILTER ORDER L */ /* 2 INVALID SECTION NUMBER K */ #ifndef KR void spbwcf(long *l, long *k, long *ln, float *d, float *c, long *error) #else void spbwcf(l, k, ln, d, c, error) long *l, *k, *ln, *error; float *d, *c; #endif { /* Local variables */ long i; float tmp; if (*l <= 0) { *error = 1; return; } if (*k <= 0 || *k > ((*l + 1) / 2)) { *error = 2; return; } d[0] = 1.0; c[0] = 1.0; for (i = 1 ; i <= *ln ; ++i) { d[i] = 0.0; c[i] = 0.0; } tmp = (float) *k - (*l + 1.0) / 2.0; if (tmp == 0.0) { c[1] = 1.0; } else { c[1] = (float) (-2.0 * cos((double) ((*k * 2) + *l - 1) * M_PI / (double) (*l * 2))); c[2] = 1.0; } *error = 0; return; } /* spbwcf */ 还有其它的几个例子
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值