我想讲单值计算"D=150" 改为矩阵计算"D =1:1:300". 但是如果直接修改的话,后面的计算就不正确了,无法显示全部结果. 所以想请问下应该使用什么code将计算修改.
- clc;
- clear all;
- Profit_S_B = [];
- Profit_S_Q = [];
- Profit_S_R =[];
- Q1 = 100;
- q_01 = 50;
- P_r = 400;
- c = 20;
- g = 50;
- x = norminv(g/(g+c),100,30);
- theta = 0.8351;
- Ps_B = ((1-theta)*P_r*g+P_r*c-c*c)/(g+c);
- Ps_Q = (P_r*c-c*c)/(g+c);
- Ps_R = (theta*P_r*c-c*c)/(g+c);
- D =150; % Here if I want change D into For D =1:1:300
- S_N0 = min(x,D);
- I_N0 = max(0,(D-x));
- S_L0 = min((D-S_N0),(Q1+q_01-x));
- I_L0 = max(0,(D-Q1-q_01));
- Profit_S_B = Ps_B*x + (1-theta)*P_r*S_N0 + (P_r-c)*S_L0 - I_L0*g-(1-theta)*P_r*x
- Profit_S_Q =Ps_Q*x + (P_r-c).*S_L0-I_L0*g
- Profit_S_R =Ps_R*x+(1-theta)*P_r*S_N0+(P_r-c).*S_L0-I_L0*g