matlab将符号推导表达式转换为函数

Matlab将符号推导表达式转换为函数

平时使用Matlab进行符号推导时，常常遇到这样一个问题，即如何能将推导出的符号表达式转换为函数，无论是求值还是绘图都比较方便。否则，当推导出的符号表达式很复杂时，如果手动录入，很麻烦也很容易出错。最终，想到了这样一个方法，即先将符号表达式转换为字符串，再由字符串转换为函数。
例如，分析一个二阶低通滤波器

利用matlab推导其传递函数：

syms R1 R2 C1 C2 s;
ZC1=1/(s*C1);
ZC2=1/(s*C2);
Zseries1=ZC2+R2;
Zparallel1=1/(1/Zseries1+1/ZC1);
g=Zparallel1/(R1+Zparallel1)*ZC2/Zseries1;

完成后，如果我们想绘制不同参数下的Bode图，则需要将其转换为函数（当然，利用tf函数去推导，再用bode命令也可以实现）。利用之前提到的思路，将符号表达式转换为字符串，再转换为函数：

g=char(g); %转换为字符串
g=replace(g,'*','.*'); %为了使转换后的函数具备向量计算能力
g=replace(g,'/','./');
g=replace(g,'^','.^');
g=eval(['@(R1,R2,C1,C2,s) ',g]);

此时，g转换为函数：

可以利用该函数绘制不同参数下的Bode图：

f=logspace(1,5,1000);
w=1j*2*pi*f;
R1=100;
C1=100e-9;
C2=200e-9;
R2_list=[50,100,150,200];
fig1=figure(1);
clf(fig1);
num=length(R2_list);
p=[];
h=cell(1,num);
for i=1:1:num
    R2=R2_list(i);
    p=[p,semilogx(f,abs(g(R1,R2,C1,C2,w)),'LineWidth',3)];
    h{i}=['R2=',num2str(R2),'\Omega'];
    hold on;
end
legend(p,h);
xlabel('Frequency/Hz');
ylabel('Magnitude');
grid minor;