嘿，我碰到过类似的情况——CFtool默认的函数库确实不包含square这种生成方波的内置函数，直接输入肯定会报错。给你几个实用的解决思路：

• 自定义函数导入CFtool
  MATLAB的CFtool支持导入用户自定义函数，按这个步骤来就行：

  1. 先写一个符合CFtool格式的自定义函数脚本，把square函数包进去，自变量t放前面，拟合参数f放后面：

     function y = custom_square_filter(t, f)
         y = square(2*pi*f*t) + 1;
     end
     

  2. 把这个脚本保存到MATLAB当前工作目录（或者添加到MATLAB搜索路径里）。
  3. 打开CFtool，选择「Custom Equation」，然后点击「Import」，找到你刚保存的.m文件，导入后就能用这个函数拟合数据了。

• 用阶跃函数组合替换方波表达式
  如果不想写自定义函数，也可以把square转换成CFtool能识别的基础表达式。因为square(x)是周期2π的方波，square(2πft)+1其实就是取值为2（高电平）和0（低电平）的方波，用取模和阶跃逻辑就能表示：

  y = 2*(mod(2*f*t, 1) < 0.5);
  

  把这个直接输入到CFtool的自定义方程框里，记得把f设为拟合参数，t设为自变量就行。

• 先滤波再拟合（适合先处理数据的场景）
  要是你只是想用这个函数先滤波再做后续拟合，不如先在MATLAB里手动生成滤波序列，处理完原始数据再导入CFtool：

  % 假设t是你的时间向量，f可以先给个初始猜测值（比如从FFT结果里找）
  filtered_signal = square(2*pi*f*t) + 1;
  % 之后把filtered_signal和你的响应数据一起导入CFtool分析
  

另外提醒一句：方波是不连续函数，CFtool拟合这类函数时容易收敛困难，建议先通过FFT分析原始数据的主频率，把这个值作为f的初始猜测，拟合成功率会高很多～

内容的提问来源于stack exchange，提问作者Rolly