Huber function matlab plotting with two intervals

2015 3 31
4 个回答
回答已采纳
更新时间：2017 11 15
4 次查看（30 天）

0 个投票

Hello,
I was trying to write this simple matlab function:
rho(eps)=0.5*eps^2 for eps<=k
rho(eps)=k*eps-0.5*k^2 for eps>k
for k=1.345
so my solution was:
k=1.345;
for eps=0:0.001:6
if eps<=k
f1=0.5*(eps^2);
else
f1=k*eps-0.5*(k^2);
end %if loop
end %for loop
plot(eps,f1)
There was almost empty graph and sometimes only linear function although this function is supposed to be nonlinear in the first part and linear in the other part.
and also i'm supposed to plot the derivative of this function after this is complete.

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

 采纳的回答

Torsten
Torsten 2015-3-31

0 个投票

k=1.345;
eps=zeros(6001);
f1=zeros(6001);
for i=1:6001
eps(i)=(i-1)*0.001;
if eps(i)<=k
f1(i)=0.5*(eps(i)^2);
else
f1(i)=k*eps(i)-0.5*(k^2);
end %if loop
end;
plot(eps,f1)
Best wishes
Torsten.

更多回答（3 个）

Kentman
Kentman 2015-3-31

0 个投票

Is there a way to plot the derivation of this function directly by matlab or we need to do it by hand then apply the result?
Torsten
Torsten 2015-3-31

0 个投票

Or more "MATLAB-like":
k=1.345;
eps=0:0.001:6;
f1=zeros(size(eps));
epsf=eps<=k;
f1(epsf)=0.5*eps(epsf).^2;
epsf=eps>k;
f1(epsf)=k*eps(epsf)-0.5*k^2;
plot(eps,f1);
Best wishes
Torsten.
Kentman
Kentman 2015-3-31

0 个投票

so the final matlab will be
clear all;close all;clc;
k=1.345; eps=zeros(6001); f1=zeros(6001); f2=zeros(6001); for i=1:6001 eps(i)=(i-1)*0.001; if eps(i)<=k f1(i)=0.5*(eps(i)^2); f2(i)=eps(i); else f1(i)=k*eps(i)-0.5*(k^2); f2(i)=k; end %if loop end; %for loop
plot(eps,f1,'k-') hold on plot(eps,f2,'b-')

类别

在 帮助中心 和 File Exchange 中查找有关 R Language 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by