tempdirObj = fidemo.fiTempdir（'fi_instrumentation_fixed_point_filter_demo'）;

copyfile（fullfile（matlabroot，'toolbox'，' fixedpoint '，'fidemos'，'+ fidemo'，... 
                  'fi_2nd_order_df2t_filter.m'），'。'，'f'）;

FIPREF_STATE = get（fipref）;
复位（fipref）

N = 256;
x = fi（零（N，1），1,16,15）;

F = fimath（'RoundingMethod'，'Floor'，...... 
           'OverflowAction'，'Wrap'，...... 
           'ProductMode'，'KeepLSB'，...... 
           'ProductWordLength'，40，...... 
           'SumMode'，'KeepLSB'，...... 
           'SumWordLength'，40）;

num = [0.0299545822080925 0.0599091644161849 0.0299545822080925];
den = [1 -1.4542435862515900 0.5740619150839550];

b = fi（num，1,16）; b.fimath = F;
a = fi（den，1,16）; a.fimath = F;

B = coder.Constant（b）;
A = coder.Constant（a）;

yisd = fi（零（N，1），1,16,15，'DataType'，'ScaledDouble'，'fimath'，F）;
zisd = fi（零（2,1），1,16,15，'DataType'，'ScaledDouble'，'fimath'，F）;

buildInstrumentedMex fi_2nd_order_df2t_filter  ... 
    -o  filter_scaled_double  ... 
    -args  {B，A，X，yisd，zisd}

t = linspace（0,1，N）;       ％时间矢量从0到1秒 
f1 = N / 2;                  ％啁啾的目标频率设置为Nyquist 
xchirp = sin（pi * f1 * t。^ 2）;  ％线性啁啾从0到Fs / 2 Hz，1秒 
x（:) = xchirp;             ％将啁啾投射到定点

ychirp = filter_scaled_double（x，yisd，zisd）;

CLF
情节（t，x，'c'，t，ychirp，'bo-'）
标题（'唧唧'）
图例（'输入'，'缩放双输出'）
图（GCF）; 的DrawNow;

showInstrumentationResults filter_scaled_double  -proposeFL

xstep = [ones（N / 2,1）;  -  ones（N / 2,1）];
x（:) = xstep;

ystep = filter_scaled_double（x，yisd，zisd）;

CLF
情节（t，x，'c'，t，ystep，'bo-'）
标题（'步骤'）
图例（'输入'，'缩放双输出'）
图（GCF）; 的DrawNow;

showInstrumentationResults filter_scaled_double  -proposeFL

yi = fi（零（N，1），1,16,14，'fimath'，F）;
zi = fi（零（2,1），1,16,14，'fimath'，F）;

buildInstrumentedMex fi_2nd_order_df2t_filter  ... 
    -o  filter_fixed_point  ... 
    -args  {B，A，x，yi，zi}

x（:) = xchirp;
[y，z] = filter_fixed_point（x，yi，zi）;
[ysd，zsd] = filter_scaled_double（x，yisd，zisd）;
err = double（y） -  double（ysd）;

CLF
subplot（211）; plot（t，x，'c'，t，ysd，'bo-'，t，y，'mx'）
xlabel（'Time（s）'）;
ylabel（'幅度'）
图例（'输入'，'缩放双输出'，'定点输出'）;
标题（'定点啁啾'）
subplot（212）; plot（t，err，'r'）; title（'Error'）; xlabel（'t'）; ylabel（'err'）;
图（GCF）; 的DrawNow;

showInstrumentationResults filter_fixed_point

clearInstrumentationResults filter_fixed_point

x（:) = xstep;
[y，z] = filter_fixed_point（x，yi，zi）;
[ysd，zsd] = filter_scaled_double（x，yisd，zisd）;
err = double（y） -  double（ysd）;

CLF
subplot（211）; plot（t，x，'c'，t，ysd，'bo-'，t，y，'mx'）
标题（'定点步骤'）;
图例（'输入'，'缩放双输出'，'定点输出'）
subplot（212）; plot（t，err，'r'）; title（'Error'）; xlabel（'t'）; ylabel（'err'）;
图（GCF）; 的DrawNow;

showInstrumentationResults filter_fixed_point

fipref（FIPREF_STATE）;
clearInstrumentationResults filter_fixed_point 
clearInstrumentationResults filter_scaled_double 
明确fi_2nd_order_df2t_filter_fixed_instrumented 
明确fi_2nd_order_df2t_filter_float_instrumented

tempdirObj.cleanUp;
％＃确定<* ASGLU>

生成定点 C 代码

此示例说明如何为简单的函数生成代码，该函数将两个输入值相乘并累加。这是可以嵌入在外部硬件中的代码的类型。函数是

function acc = mult_acc(x,a,acc)
acc = accumpos(acc,x*a); 

以下代码定义测试平台输入，设置必需的代码生成属性并生成代码。测试平台输入被指定为定点数。x 输入为随机数，a 为 0.9，累加器 acc 初始化为 0。coder.HardwareImplementation 对象指定影响生成代码的外部硬件的属性。这些示例指定一个 40 位累加器。coder.CodeConfig 对象具有直接影响代码生成的属性。codegen 命令将函数、配置对象作为输入参数并生成可嵌入的 C 代码。

x = fi(rand,true,16,15);
a = fi(0.9,true,16,15);
acc = fi(0,true,40,30);


%% 
hi = coder.HardwareImplementation;
hi. ProdHWDeviceType = 'Generic->Custom'
hi. TargetHWDeviceType = 'Generic->Custom'
hi.TargetBitPerLong = 40;
hi.ProdBitPerLong   = 40;

hc = coder.config('lib');
hc.HardwareImplementation = hi;
hc.GenerateReport         = true;

codegen mult_acc -config hc -args {x,a,acc}

生成的 C 代码是

/*
 * mult_acc.c
 *
 * Code generation for function 'mult_acc'
*/

/* Include files */
#include "rt_nonfinite.h"
#include "mult_acc.h"

/* Function Definitions */
void mult_acc(short x, short a, long *acc)
{
  *acc += x * a;
}

/* End of code generation (mult_acc.c) */

function y = mysum(x)
  y = 0;
  for n = 1:length(x)
    y = y + x(n);
  end
end

function y = mysum(x) %#codegen
  y = 0;  
  for n = 1:length(x)
    y = y + x(n);
  end
end

function y = mysum(x,T) %#codegen
  y = cast(0,'like',T.y);
  for n = 1:length(x)
    y(:) = y + x(n);
  end
end

function T = mytypes(dt)
  switch dt
    case 'double'
      T.x = double([]);
      T.y = double([]);
    case 'single'
      T.x = single([]);
      T.y = single([]);
    case 'fixed'
      T.x = fi([],true,16,15);
      T.y = fi([],true,16,13);
    case 'scaled'
      T.x = fi([],true,16,15,...
           'DataType','ScaledDouble');
      T.y = fi([],true,16,13,...
           'DataType','ScaledDouble');
  end
end