Please fix my code

지원 안

2021 5 16

0 个回答

1 次查看（30 天）

0 个投票

to solve y'=t(y.^2-4)

I try like this.

but something was wrong... Amend pls......

clear all
[x,y]=meshgrid(linspace(0,3,100));
u=ones(100);
v=(y.^2-4)*(1/2).*(log(abs(2-y/2+y)).^(1/2));
quiver(x,y,u,v,0.8)
hold on
 
t=0, w=1;
dt=0.01;
N=1/dt;
 
for n=1:N
    t=t+dt;
    w=w+dt*(y.^2-4)*(1/2).*(log(abs(2-y/2+y)).^(1/2));
    st=[st t
end
plot(st,sw,'r')

what is wrong?

4 个评论
显示 2更早的评论隐藏 2更早的评论

지원 안 2021-5-17

clear all

[x,y]=meshgrid(linspace(0,3,100));

u=ones(100);

v=(y.^2-4)*(1/2).*(log(abs(2-y/2+y)).^(1/2));

quiver(x,y,u,v,0.8)

hold on

t=0, w=1;

dt=0.01;

N=1/dt;

st=t; sw=w;

for n=1:N

t=t+dt;

w=w+dt*(y.^2-4)*(1/2).*(log(abs(2-y/2+y)).^(1/2));

st=[st t]; sw=[sw w];

end

plot(st,sw,'r')

but something wrong... what's wrong?

DGM 2021-5-17

在 MATLAB Online 中打开

What's the error message say?

Error using horzcat
Dimensions of arrays being concatenated are not consistent.
    sw=[sw w];

so look at the size of sw and w

>> size(sw)
ans =
     1     1
>> size(w)
ans =
   100   100

why is sw the size it is? It's initialized as 1. You could initialize it as an empty vector [], but pay attention to the size of w.

w=1;
% ...
sw=w;

why is w the size it is?

[x,y]=meshgrid(linspace(0,3,100));
% ...
w=w+dt*(y.^2-4)*(1/2).*(log(abs(2-y/2+y)).^(1/2));

w is 2D. You're going to figure out how you want to plot these. Once you figure that out, you'll have to figure out how you need to store them. If you need to store all instances of them, you can concatenate them on dim3.

请先登录，再进行评论。

请先登录，再回答此问题。

Follow Question