Code error with parameter

2020 1 4
1 个回答
回答已采纳
更新时间：2020 1 4
1 次查看（30 天）

0 个投票

I'm trying to adapt my code below
rg = 0.01;
rl = 0.01;
s1 = [-1:0.000001:5];
x = ((s1-rg-rl)+sqrt((rg+rl-s1).^2+4.*s1.*rg))./(2.*s1);

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

 采纳的回答

Star Strider
Star Strider 2020-1-4

0 个投票

I am not certain what you want, however it most likely depends on how you define ‘t’.
Example —
rg = 0.01;
rl = 0.01;
s1 = [-1:0.000001:5];
x = ((s1-rg-rl)+sqrt((rg+rl-s1).^2+4.*s1.*rg))./(2.*s1);
t1 = 0:numel(s1)-1; % Define ‘t1’, Interval = 1 sec
t2 = (0:numel(s1)-1)*5; % Define ‘t2’, Interval = 5 sec
figure
plot(t1, x)
hold on
plot(t2, x)
hold off
grid
legend('\Delta t = 1 s', '\Delta t = 5 s', 'Location','E')

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

Ellen Brown
Ellen Brown 2020-1-4
I'm trying to get it to say e.g. s1 is drawn from an exponential distribution with mean 1 and s2 is drawn from an exponential distribution with mean 2, x will be calculated used the two states 's1' and 's2' , and which one it is calculated one from changes periodically every t time units.
Star Strider
Star Strider 2020-1-4
Are you intending to simulate a Poisson process?
If so, see if the poisspdf function willl do what you want.
Ellen Brown
Ellen Brown 2020-1-4
sorry I'm not explaining it clearly,
my s is simulated from a poisson:
U = rand(1,1000); % generate 1000 uniform (0,1) random variables
% use inverse transform method to generate 1000 exp rv with rate 1, mean
% shifted to -1
s = -log(1-U)-1;
but i want s to switch between two independent states s1 and s2, were s1 and s2 are simulated from:
s1 = -log(1-U)-1;
s2 = -log(1-U)-1;
I want to calculate x as s changes between two states, “s1” and “s2”, periodically every t time units, where x is:
x = ((s-rg-rl)+sqrt((rg+rl-s).^2+4.*s.*rg))./(2.*s);
Star Strider
Star Strider 2020-1-4
It appears that ‘s1’ and ‘s2’ are identical:
s1 = -log(1-U)-1;
s2 = -log(1-U)-1;
I’m still very much lost.
See if poissinv will do what you want.
Ellen Brown
Ellen Brown 2020-1-4
编辑：Ellen Brown 2020-1-4
ah typo, s2 is meant to be:
s2 = -log(1-U)+1;
I just want x to change which s it uses periodically, i.e. for first 5s it uses s1, then the next 5s it uses s2, then s1 again, e.t.c.
Star Strider
Star Strider 2020-1-4
I still have no clear idea of what you want.
The only thing I can think of is to use a loop:
rg = 0.01;
rl = 0.01;
U = rand(1,1000);
s1 = -log(1-U)-1;
s2 = -log(1-U)+1;
x = @(s) ((s-rg-rl)+sqrt((rg+rl-s).^2+4.*s.*rg))./(2.*s);
xs = zeros(size(U));
for k = 1: numel(U)
s_val = s1(k)*(rem(k,2)==0) + s2(k)*(rem(k,2)==1); % Even-Odd Values of ‘k’ Select ‘s’
xs(k) = x(s_val);
end
I am now completely out of ideas.
Ellen Brown
Ellen Brown 2020-1-4
That's exactly what I want! Thank you so much! I'm being a bit cheeky now but how would I get a vector of those s_val values? I'd like to plot xs against the s_val values
Star Strider
Star Strider 2020-1-4
As always, my pleasure!
Not cheeky at all! That’s a perfectly legitimate question.
To store the ‘s_val’ values (without re-writing the rest of the loop):
xs = zeros(size(U));
s_val_vct = xs;
for k = 1: numel(U)
s_val = s1(k)*(rem(k,2)==0) + s2(k)*(rem(k,2)==1); % Even-Odd Values of ‘k’ Select ‘s’
xs(k) = x(s_val);
s_val_vct(k) = s_val;
end
Then to plot them:
figure
plot(s_val_vct, xs, 'p')
grid
xlabel('s')
ylabel('x')
It’s best to plot them with markers, unlless you want to use the sort function to sort ‘xs’ and ‘s_val_vct’ and then plot with lines, because ‘s_val’ alternates between poisitive and negative values, and that would create a confusing and disordered plot. (I like plotting with pentagrams. Choose the marker you want.)

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Graphics Object Properties 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by