可以这样改:
function Yd=power_tra(t,YY)
%发电机转子摇摆曲线方程
%t一定是标量形式的自变量
%YY必须是列向量
global y0 Tj Pt E U x1 %在函数中定义全局变量传递参数
%发电机转子摇摆曲线微分方程
Yd=[(YY(2)-1)*y0;(Pt-(E*U/x1)*sin(YY(1)))/Tj];
end
%{
%主程序,另存为main.m,再运行,或copy到命令行运行
%在主程序中定义全局变量传递参数
global y0 Tj Pt E U x1
y0=2*pi*50;Tj=11.28;Pt=1;E=1.47;U=1;
%系统转移电抗
x1=2.82;%故障时
%指定解算微分方程的时间区间
tspan=[0.0 0.3];
%给定初值向量
y1=[31.54*pi/180;1]%故障时
%求解微分方程
[t,YY]=ode45('power_tra',tspan,y1);
%输出求解结果
x=YY(:,1);
y=YY(:,2);
%曲线绘制
plot(t,x*180/pi);
xlabel('t/s');
ylabel('delta/deg');
grid on
%}
dim tit
tit=request("tit")
tit=replace(tit,"'","''")
tit=trim(tit)
tit=cstr(tit)
sql="select * from blog_Content where log_Title like '%"+ tit +"%'"