请大家帮我看看这2个函数图像怎么画啊？？

xigang

我是新手，弱弱向大家请教几个问题。
最近在写一篇桩纵向振动理论方面的论文（岩土工程），推导了一个半解析解，结果函数图像画不出来，主要问题：
一是，函数里有一个Fourier逆变换，从负无穷到正无穷积分，我用MATLAB函数ifourier()对其求解，结果表达式里仍有ifourier符号，等于是答案求解不出来。
二是，函数里有一个tg()函数，我不知道如何处理tg()出现无穷的情况。 听说要编程，很茫然，特来向大家请教。

有不少公式，不好弄上来，我把问题写成WORD文档，麻烦大家看看，小弟万分感激

xigang

>> syms k1 E p C v rg w1 t1 l Tc T1 A M real
>> E=3*10^10;
>> p=2450;
>> C=sqrt(E/p);
>> l=10;
>> Tc=l/C;
>> v=0.2;
>> rg=0.2475;
>> T=2*10^(-3);
>> T1=T/Tc;
>> A=0.3847;
>> k1=w1/(E*Tc^2/p-v^2*rg^2*w1^2)^(1/2);
>> Hv=w1*tan(k1*l)/(E*A*Tc*k1);
>> ezplot(w1,Hv1,[0,100])
>> g1=i*w1*tan(k1*l)*T1*(1+exp(-(i*T1*w1)/2))*exp(i*w1*t1)/(k1*(4*pi^2-T1^2*w1^2))
>> g=ifourier(g1,w1,t1)
g =

-246290604621824/5*i*6^(1/2)*(ifourier(tan(200*w1/(40000-10*w1^2)^(1/2))*exp(i*w1*t1)*(40000-10*w1^2)^(1/2)/(-136124286778060735+1688849860263936*w1^2),w1,t1)+ifourier(tan(200*w1/(40000-10*w1^2)^(1/2))*exp(i*w1*t1)*(40000-10*w1^2)^(1/2)/(-136124286778060735+1688849860263936*w1^2)*exp(-1/7*i*6^(1/2)*w1),w1,t1))
g为桩顶速度时域速度直，可求出的结果里面仍有ifourier符号，图像也画不出来，肯定大家指点

xigang

后来问了一些老师，说是tan()函数问题，可我怎么看就没有什么问题，是不是两边出现　无穷的情况？？
恳请高手，帮我实现一下