为什么我的程序没报错，图像却出不来？求大佬解惑！

我的代码如下：

close all;

clear all;

L=1; EI=5; T=2; rho=2; M=6;

nt=80000; nx=20; tf=15; Ttr=50;

dx=L/nx; dt=tf/nt;

x = linspace(0,L,nx+1); t_tr = linspace(0,tf,Ttr+1);

w=zeros(nx+1,nt+1);

w3D_free=zeros(Ttr+1,nx+1);

d=zeros(nt+1,1);

for j=1:nt+1

d(j) = 1 + sin(2*pi*(j-1)*dt) + cos(3*pi*(j-1)*dt) ;

end

Br=90;Bl=-90;k=100;k1=1;k2=4;r=1;m=1;

w=zeros(nx+1,nt+1); w3D_model=zeros(Ttr+1,nx+1);

u_model=zeros(nt+1,1);u3D_model=zeros(Ttr+1,1);

ua=zeros(nt+1,1); d_estimate=zeros(nt+1,1);

d1=zeros(nt+1,1);

dbar=max(abs(d));

d_estimate(1