matlab 复现文献中的图，用到求解未知数

帮别人做的文献中的程序，原文见
[1]. 江岳文与陈晓榕, 风电波动成本分摊方法. 电力自动化设备, 2020.
要求根据以下原理作图

求等效出力曲线图，如下

想到了两种求解思路，第一种求解思路和文献中的图有出入，编程思路是先求出l和pw的值，满足公式（1），然后讲求得的和相除得到比例因子，然后将l中所有的纵坐标除以比例因子，即得到要求得的pw1
代码如下

第二种思路是严格按照公式（1）和公式（2）中的原来编程，设未知数x，迭代求得pw1的和表达式，然后利用matlab来方程的解，然后将解带入到pw1中
，代码实现如下

clear;clc
x=xlsread('wind2.xls');
pw=x(:,1);
l=x(:,2);
syms x
% k=sum(pw);
pw1=[];
pw1=[pw1;x];


for i=2:length(pw)
   
       q=(((l(i)-l(i-1))/l(i-1))*pw1(i-1))+pw1(i-1);
       pw1=[pw1;q];
    
    
i
end
res=solve(sum(pw1)-sum(pw));
x=res;
pw1=eval(pw1)

建设使用的时候用第二段代码，严格按照文献中的原理编程，缺点是计算缓慢，数据量大的时候求未知数可能有误差