Plotting a Piecewise function (2024)

Direkter Link zu diesem Kommentar

https://de.mathworks.com/matlabcentral/answers/89871-plotting-a-piecewise-function#comment_173634

t = (1:56)

while t>=1 & t<=28

Ws

t=48+3.6*t+.6363*t^2+.00963*t^3

end

while all(t>=1 & t<=28)

t = 1:56;

tt(1:28) = 48+3.6*t+.6363*t^2+.00963*t^3;

tt(29:56) = -1004+65.8*t;

function graph;W (1,56)

function graph

N=1000;

n=(56-1)/N;

x=[1:n:56];

for i=1:length(x)

F(i)=W(x(i));

end

plot(x,F)

xlabel('Time')

ylabel('Body Weight')

end

function F = W(t)

if t>=1 && t<=28

Ws

F=48+3.6*t+.6363*t^2+.00963*t^3

end

if 28<t && t<=56

Wr

F=-1004+65.8*t

end

function Ws

disp('t is greater than or equal to 1 and t is less than or equal to 28')

end

function Wr

disp ('t is greater than 28 and less than or equal to 56')

end

end

j = 1:56;

t = zeros(1, length(j));

for i = 1:56

if (i >= 1 && i < 29)

disp('t is greater than or equal to 1 and t is less than or equal to 28')

t(i) = 48 + 3.6 * i + .6363 * i^2 + .00963 * i^3;

else if (i > 28 && i <= 56)

disp ('t is greater than 28 and less than or equal to 56')

t(i) = -1004 + 65.8 .* i;

end

end

end

N = 1000;

x = min(j):((max(j) - min(j)) / N):56;

F = interp1(j, t, x);

plot(x, F), xlabel('Time'), ylabel('Body Weight')