Rectangular Function implementation in Matlab

Question

Rohan Repale 2013-3-20

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/68012-rectangular-function-implementation-in-matlab

评论： Vamshidhar Reddy 2022-11-15

采纳的回答： Azzi Abdelmalek

How do I implement this equation?

rect(x)

Thanks

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Azzi Abdelmalek 2013-3-20

4
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/68012-rectangular-function-implementation-in-matlab#answer_79386

编辑：Azzi Abdelmalek 2013-3-20

在 MATLAB Online 中打开

rect=@(x,a) ones(1,numel(x)).*(abs(x)<a/2) % a is the width of the pulse
x=-10:0.001:10
y=rect(x,2)   
plot(x,y)

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

Khushi Bhatti 2018-9-9

can someone please tell me how can i use a rectangular function or any function to visualize outliers in my data? i find the outliers in my data but i dnt know how to visualize them using a graph or plot.

Vamshidhar Reddy 2022-11-15

rectangular function do require symbolic toolbox

请先登录，再进行评论。

Answer 2

Johan Rincon Botia 2022-4-16

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/68012-rectangular-function-implementation-in-matlab#answer_944060

在 MATLAB Online 中打开

I did the implementation for rect function with period w

My function has 4 arguments rect (T,n,amp,xtranslation), where T is the period, n is the periods that you want to show, amp is the amplitude and xtranslation is the translation at x-axis, for example, i want to do this signal:

with any period, like T = 20s, and i wanna show this 3 times,

This is how i use the function

T = 20; %Period
amplitude = 5; %Amplitude of the signal
num = 3; %Numbers of periods that you want to show
xtranslation =-8; % Translation x-axis
[t,y] = rect(T,num,amplitude,xtranslation); %the function return the time vector and amplitud vector
plot(t,y)%use t and y to plot your rect signal

An here is the implementation of the function:

function [t,y] = rect(T_period,num,ampli,xtranslation)
    t = -T_period*num:0.01:T_period*num;
    y_0=0;
    y_1=1;
    if xtranslation>=T_period/2
        num = num+round((2*xtranslation/T_period))
    end
    for i=-num:num
        y_0=y_0 + 0.*(t>T_period/2+(T_period*(i-1)) & t<T_period*num);
        y_1=y_1 + 1.*(t>T_period*(i-1) & t<T_period/2+(T_period*(i-1)));
    end
    y= (y_0 + y_1)
    y = ampli*(y-min(y));
    t = t-xtranslation;
end