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;