
data0 = load('Competition_test')%data0是个struct结构,包含X域
A = data0.X%把X域赋给A
B = A(:,:,1)%取A的第一页给B
% size(B)%如果前面的“%”去掉,会看到B的大小,可能是100*64
C = zeros(6400,1)%生成全0矩阵,
for i = 1:100
for j = 1:64
C(i*j) = B(i,j,1) % 提取脑电信号
end
end
ode45 是解微分方程的函数,你在帮助里的说明书点ode45就能查到。
我这给你写一段,假设有微分方程
你要新建一个函数 写出这个微分方程,然后在主程序里调用这个函数,然后加入初始值,就是0 1 1。 然后用ode解。
这是定义的函数的样子:
function dy = rigid(t,y)
dy = zeros(3,1) % a column vector
dy(1) = y(2) * y(3)
dy(2) = -y(1) * y(3)
dy(3) = -0.51 * y(1) * y(2)
然后在主函数里,也就是另一个M文件里,解微分方程,
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5])
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options)
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.');%这个是绘制图形
最后那句就是画图,出来的图形就是y的值,y1 y2 y3 .如下图:
希望解决了你的疑惑。谢谢。
欢迎分享,转载请注明来源:内存溢出
微信扫一扫
支付宝扫一扫
评论列表(0条)