入门基础-动态规划

最新推荐文章于 2022-05-27 19:47:35 发布
最新推荐文章于 2022-05-27 19:47:35 发布
阅读量166 收藏 1
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.youkuaiyun.com/qq_52155014/article/details/119106932
版权

https://www.luogu.com.cn/problem/P1060

P1060 [NOIP2006 普及组] 开心的金明

一道典型的背包问题

思路：先说二维的dp[i][j],i表示0---i件物品是花费j元的情况，也就是说先从第一层（i=1）的情况都更新，然后再更新（i=2）的情况，i=2的情况是在i=1的情况下推倒出来的，因此我们就可以得到

一维转移方程：f[j]=max(f[j],f[j-v[i]]+v[i]*w[i])

dfs:思路，选或者不选，每个都枚举一遍，要注意最后意见物品的选择（如果我们选择最后一件物品，那么number会出现number>m的情况，因此我们要在判断位置之前要更新我们的答案）以及花钱的边界。

dfs：

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
ll value[3000005],spend[3000005],n,m;
int ans=0;
void dfs(int number,int money,int jiazhi)
{
	if(money>0)ans=max(jiazhi,ans);
	if(number>m)return;
	dfs(number+1,money-spend[number],jiazhi+value[number]);//选 
	dfs(number+1,money,jiazhi);//不选 
}
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		cin>>spend[i]>>value[i];
		value[i]=spend[i]*value[i];
	}
	dfs(1,n,0);
	cout<<ans<<endl;
	return 0;
}

二维代码：转移方程 dp[i][j]=max(dp[i-1][j],dp[i-1][j-spen[i]]+value[i]);

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
ll dp[30][3000005],n,m,ans,value[3000005],spend[30000005];
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		cin>>spend[i]>>value[i];
		value[i]=spend[i]*value[i];
	}
	for(int i=1;i<=m;i++)
	{
		for(int j=1;j<=n;j++)
		{
			if(j>=spend[i]) 
			{
				dp[i][j]=max(dp[i-1][j],dp[i-1][j-spend[i]]+value[i]);
			}
			else dp[i][j]=dp[i-1][j];
		}
	}
	cout<<dp[m][n]<<endl;
	return 0;
}

一维代码：一维转移方程：f[j]=max(f[j],f[j-v[i]]+v[i]*w[i])

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
ll dp[3000005],n,m,ans,value[3000005],spend[30000005];
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>n>>m;
	for(int i=1;i<=m;i++)
	{
		cin>>spend[i]>>value[i];
		value[i]=spend[i]*value[i];
	}
	for(int i=1;i<=m;i++)
	{
		for(int j=n;j>=spend[i];j--)
		{
		dp[j]=max(dp[j],dp[j-spend[i]]+value[i]);
		}
	}
	cout<<dp[n]<<endl;
	return 0;
}

https://www.luogu.com.cn/problem/P1164

P1164 小A点菜

dfs：普通的dfs的时间复杂度是：2^100,显然如果都枚举完会超时，那么我们就需要通过剪枝来优化我们的复杂度，就能过了。

dfs+剪枝代码：

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
ll ans=0,n,m,spend[1005],dp[105][100005];
ll dfs(int i,int m)
{
	if(i>n||m<0)return 0;//边界 
	if(dp[i][m])return dp[i][m];//剪枝 
	if(m==spend[i])dp[i][m]+=1;//如果剩余的钱刚好够钱点第i个那么我们就让总的方案书+1 
	dp[i][m]+=dfs(i+1,m-spend[i]);//点第i个套餐 
	dp[i][m]+=dfs(i+1,m);//不点第i个套餐 
	//总方案数=点第i个套餐的方案+不点第i个套餐的方案，i从1开始，因此前面每个i的情况都枚举完了。 
	return dp[i][m];
}
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>n>>m;
	for(int i=1;i<=n;i++)cin>>spend[i];
	dp[n][m]=dfs(1,m);
	cout<<dp[n][m];
	return 0;
}

二维dp：转移方程：if(j>=spend[i]) dp[i][j]=dp[i-1][j]+dp[i-1][j-spend[i]];

                            else                     dp[i][j]=dp[i-1][j];

代码：

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
ll dp[105][10005],n,m,spend[10005]; 
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>n>>m;
	for(int i=1;i<=n;i++)cin>>spend[i];
	for(int i=0;i<=n;i++)dp[i][0]=1;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			if(j<spend[i])dp[i][j]=dp[i-1][j];
			else dp[i][j]=dp[i-1][j]+dp[i-1][j-spend[i]];
		}
	}
	cout<<dp[n][m]<<endl;
	return 0;
}

一维：dp[0]=1, dp[j]=dp[j]+dp[j-value[i]];

一维代码：

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
ll dp[10005],n,m,spend[10005]; 
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>n>>m;
	for(int i=1;i<=n;i++)cin>>spend[i];
	dp[0]=1;
	for(int i=1;i<=n;i++)
	{
		for(int j=m;j>=spend[i];j--)
		{
			dp[j]=dp[j]+dp[j-spend[i]];
		}
	}
	cout<<dp[m]<<endl;
	return 0;
}

https://www.luogu.com.cn/problem/P1115

P1115 最大子段和

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
int f[211111],n;
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		int x;
		cin>>x;
		f[i]=max(f[i-1]+x,x);//dp
	}
	int maxx=-999999;
	for(int i=1;i<=n;i++)
	{
		maxx=max(f[i],maxx);
	}
	cout<<maxx<<endl;
	return 0;
}

https://www.luogu.com.cn/problem/P1734

P1734 最大约数和

dp题

一维转移方程：dp[j]=max{dp[j],dp[j-i]+yue[i]};

如果用二维来写，那么需要注意的是在第二个循环里，我们必须对每一个数都更新，因为在下一层状态可能需要到本层的状态，这里我被卡了好久...

二维转移方程：if（j<i）dp[i][j]=dp[i-1][j];

                         else      dp[i][j]=max{dp[i-1][j],dp[i-1][j-i]+yue[i]};

一维代码：

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
int yueshu[1005],dp[1005],n;
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		int temp=0;
		for(int j=1;j<=i/2;j++)
		{
			if(i%j==0)
			{
				temp+=j;
			}
		}
		yueshu[i]=temp;
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=n;j>=i;j--)
		{
			dp[j]=max(dp[j],dp[j-i]+yueshu[i]);
		}
	}
	cout<<dp[n]<<endl;
	return 0;
}

二维代码：

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
int yueshu[1005],dp[1005][1005],n;
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		int temp=0;
		for(int j=1;j<=i/2;j++)
		{
			if(i%j==0)
			{
				temp+=j;
			}
		}
		yueshu[i]=temp;
	}
	for(int i=1;i<=n;i++)
	{
	    for(int j=1;j<=n;j++)
	    {
	    	if(j<i)dp[i][j]=dp[i-1][j];
	        else dp[i][j]=max(dp[i-1][j],dp[i-1][j-i]+yueshu[i]);
	    }
	}
	cout<<dp[n][n];
	return 0;
}

https://www.luogu.com.cn/problem/P1002

P1002 [NOIP2002 普及组] 过河卒

画图，转移方程：dp[i][j]=max(d[i][j],dp[i-1][j]+dp[i][j-1]);

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
ll dp[25][25],pd[25][25],mx,my,zx,zy;
int dx[]={0,1,1,-1,-1,-2,-2,2,2};
int dy[]={0,2,-2,2,-2,1,-1,1,-1};
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>zx>>zy>>mx>>my;
	zx=zx+2;zy=zy+2;mx=mx+2;my=my+2;
	for(int i=0;i<=8;i++)
	{
		int xx=mx+dx[i];
		int yy=my+dy[i];
		pd[xx][yy]=1;
	}
	dp[2][2]=1;
	for(int i=2;i<=zx;i++)
	{
		for(int j=2;j<=zy;j++)
		{
			if(pd[i][j])continue;
			else dp[i][j]=max(dp[i][j],dp[i-1][j]+dp[i][j-1]);
		}
	}
	cout<<dp[zx][zy]<<endl;
	return 0;
}

https://www.luogu.com.cn/problem/P1616

P1616 疯狂的采药

转移方程：dp[j]=max(dp[j],dp[j-spend[i]]+value[i]);

                        答案=max（选择，不选择）；

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
ll value[20000005],spend[20000005],dp[20000005],n,t;
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>t>>n;
	for(int i=1;i<=n;i++)
	{
		cin>>spend[i]>>value[i];
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=t;j++)
		{
			if(j>=spend[i])dp[j]=max(dp[j],dp[j-spend[i]]+value[i]);
		}
	}
	cout<<dp[t];
	return 0;
}

https://www.luogu.com.cn/problem/P1048

P1048 [NOIP2005 普及组] 采药

转移方程：dp[j]=max(dp[j],dp[j-spend[i]+value[i]);

#include<bits/stdc++.h>
#define ll long long
#define rep(a,b,c) for(int a=b;a<=c;a++)
#define per(a,b,c) for(int a=b;a>=c;a--)
#define hh 0x3f3f3f3f
using namespace std;
inline int read()
{
	ll s=0,f=1;
	char c=getchar();
	while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9'){s=(s<<3)+(s<<1)+c-'0';c=getchar();}
	return s*f;
}
ll dp[100000],n,t,value[100000],spend[100000];
int main()
{
	ios::sync_with_stdio(0);cout.tie(0);cin.tie(0);
	cin>>t>>n;
	for(int i=1;i<=n;i++)cin>>spend[i]>>value[i];
	for(int i=1;i<=n;i++)
	{
		for(int j=t;j>=spend[i];j--)
		{
			dp[j]=max(dp[j],dp[j-spend[i]]+value[i]);
		}
	}
	cout<<dp[t];

	return 0;
}

夜深人静写算法（二）- 动态规划入门
英雄哪里出来
09-05 2万+
全网独家《动态规划》从入门到精通
算法学习笔记（8.2）-动态规划入门进阶
2301_76874968的博客
07-09 978
算法学习笔记（8.2）-动态规划入门进阶
一些常用define
zawedx的博客
02-15 815
#define rep(a,b,c) for (int a=b;a&lt;=c;a++) #define per(a,b,c) for (int a=b;a&gt;=c;a--) #define go(u) for (int o=ft[u],v;v=E[o].t,o;o=E[o].n) // forEdge #define fi first #define se second // pair #d...
面向对象荣誉课学习笔记
自强不息 厚德载物
06-24 157
正则表达式 #include <bits/stdc++.h> #define rep(i, a, b) for(int i = (a); i < (b); i++) #define _for(i, a, b) for(int i = (a); i <= (b); i++) using namespace std; int main() { string s; regex r("([-+]*[0-9]*x\\^2)*([-+]*[0-9]*x)*([-+]*[
22/5/27
aidaqiudeaichao的博客
05-27 113
1，robot instructions； 1，robot instructions 题意：给出n条指令，和终点edx，edy，初始位置在(0,0);每条指令包含xi,yi，执行指令操作是使当前坐标x+=xi, y+=yi,问对于执行1~n次指令，可以到达终点的指令数量分别有多少；n<=40;（特别说明：时间限制和内存限制翻倍，分别是4s和512mb）； 思路：肯定是要枚举每个点的状态，选还是不选，这样的时间复杂度是O(2^n),n最大是40，这个复杂度肯定要t掉的，所以考虑优化； 对于枚举点
22/5/13
aidaqiudeaichao的博客
05-13 346
一，训练赛补题：1，Searching for Soulmates；2，Walking Home；二，cf round 789 div2补题：1，Tokitsukaze and Strange Inequality；三，acwing 提高课：acwing 188.武士风度的牛； 一，训练赛补题： 1，Searching for Soulmates 题意：给出两个数a，b，有如下操作：a*2，a/2，或者a+1，来使a==b，问最小的操作数； 思路：首先乘和除一定不会交替出现，因为会相互抵消掉的；所.
Splay平衡树 (LCT基础)
Free Loop~~~跳动的恒星
04-04 381
#include using namespace std; #define REP(i,j,k) for (int i = (j), _end_ = (k); i <= _end_; ++i) #define DREP(i,j,k) for (int i = (j), _start_ = (k); i >= _start_; --i) #define debug(...) fprintf(st
算法题-动态规划算法40道实例入门到熟练
07-03
根据给定文件的信息，我们可以提炼出以下几个重要的动态规划知识点及相关示例： ## 动态规划基础知识及解题...以上是针对这四个基础入门题目的详细解析，希望能帮助读者更好地理解和掌握动态规划的基本思想与技巧。
算法-动态规划-斐波那契模型
09-04
斐波那契数列是动态规划经典入门模型，它不仅是一个数学上的概念，更在算法学习中扮演了重要的角色。 在动态规划的学习过程中，斐波那契数列是最佳的切入点之一。斐波那契数列模型通常是指一个序列，其中每个数字是...
算法基础入门：90分钟搞定动态规划
06-17
动态规划（Dynamic programming，简称DP）很多人都觉得是比较难以理解和掌握的一种算法，为了应付面试更多的时候程序员会选择直接死记硬背斐波那楔数列或者背包问题的源码，其实只要认真学习、彻底理解，动态规划并...
用快排思路解决逆序对问题
独钓寒江雪
10-06 2271
#include #include #define rep(i, n) for(i=0;i<n;i++) #define repi(i,l,n) for(i=l;i<=r;i++) using namespace std; void quicksort(int l, int r); int a[100]; int b[100]; int hit; int main(){ int i, k;
C++ 中define的使用详解
说文科技，做有态度的研究。
04-01 4505
C++中define的使用详解 1.代码如下 #include<cstdio> #include<iostream> #define max(a,b) ((a)>(b)? (a): (b)) using namespace std; int main(){ int a = 2; int b = 2; int n = max(a,b); printf("n ...
#define新用法
里弗杰克
12-09 2295
#definerep(I,n) for(i=0;i Intmain() {          Int I,j,n,a[100];          Scanf(“%d”,&n); rep(I,n)scanf(“%d”,&a[i]); rep(j,n)printf(“%d ”,&a[j]); printf(“\n”); return0; }
AtCoder Regular Contest 123
qq_52155014的博客
08-09 261
https://atcoder.jp/contests/arc123/tasks/arc123_a A - Arithmetic Sequence 题意:给你a,b,c三个数，每次操作你可以让这三个数中的一个+1，直到a-b=b-c时，结束操作。问你怎么操作使得你的操作数最少。 思路：移项：a+c-2b=0；为了方便我们令a+b=A，2b=B；因此有A-B=0；每次操作我们可以使A+1或者让B-2，直到总和=0，这时我们的操作数是最少的。 代码： #include<bits/stdc++.
AtCoder Beginner Contest 215
qq_52155014的博客
08-22 230
传送门：A - Your First Judge 题意：给你一个字符串s，如果这个字符串完全符合“Hello,World!”就输出AC，否则就输出WA。 注意：不能有任何其他多出来的字符，比如：“Hello,World!Hello,World!”也是输出WA。 代码： #include<bits/stdc++.h> #define ll long long #define hh 0x3f3f3f3f const int maxn=1e5+5; using namespace std;
thymeleaf-2.0.15.jar中文文档.zip
05-04
# 压缩文件中包含： 中文文档 jar包下载地址 Maven依赖 Gradle依赖 源代码下载地址 # 本文件关键字： jar中文文档.zip,java,jar包,Maven,第三方jar包,组件,开源组件,第三方组件,Gradle,中文API文档,手册,开发手册,使用手册,参考手册 # 使用方法： 解压最外层zip，再解压其中的zip包，双击 【index.html】 文件，即可用浏览器打开、进行查看。 # 特殊说明： ·本文档为人性化翻译，精心制作，请放心使用。 ·只翻译了该翻译的内容，如：注释、说明、描述、用法讲解 等； ·不该翻译的内容保持原样，如：类名、方法名、包名、类型、关键字、代码 等。 # 温馨提示： （1）为了防止解压后路径太长导致浏览器无法打开，推荐在解压时选择“解压到当前文件夹”（放心，自带文件夹，文件不会散落一地）； （2）有时，一套Java组件会有多个jar，所以在下载前，请仔细阅读本篇描述，以确保这就是你需要的文件；
javassist-3.25.0-GA.jar中文-英文对照文档.zip
05-04
# 压缩文件中包含： 中文-英文对照文档 jar包下载地址 Maven依赖 Gradle依赖 源代码下载地址 # 本文件关键字： jar中文-英文对照文档.zip,java,jar包,Maven,第三方jar包,组件,开源组件,第三方组件,Gradle,中文API文档,手册,开发手册,使用手册,参考手册 # 使用方法： 解压最外层zip，再解压其中的zip包，双击 【index.html】 文件，即可用浏览器打开、进行查看。 # 特殊说明： ·本文档为人性化翻译，精心制作，请放心使用。 ·只翻译了该翻译的内容，如：注释、说明、描述、用法讲解 等； ·不该翻译的内容保持原样，如：类名、方法名、包名、类型、关键字、代码 等。 # 温馨提示： （1）为了防止解压后路径太长导致浏览器无法打开，推荐在解压时选择“解压到当前文件夹”（放心，自带文件夹，文件不会散落一地）； （2）有时，一套Java组件会有多个jar，所以在下载前，请仔细阅读本篇描述，以确保这就是你需要的文件；
“计及源荷双重不确定性的虚拟电厂微网日前随机优化调度模型设计与实现（基于蒙特卡洛与概率距离法处理不确定性变量）”
最新发布
05-04
内容概要：本文详细介绍了利用MATLAB进行虚拟电厂日前随机优化调度的方法。首先，通过蒙特卡洛模拟生成大量光伏出力和负荷的随机场景，考虑到光伏出力和负荷预测的不确定性。接着，采用概率距离削减法将大量场景缩减为少数典型场景，确保计算效率的同时保留重要特征。然后，构建随机优化模型，将燃气轮机、储能、弃光等成本纳入目标函数，并设置严格的功率平衡和设备出力限制作为约束条件。最终，使用CPLEX求解器求解该混合整数规划问题，得到最优调度方案。实验结果显示，相比传统确定性优化，随机优化显著提高了系统的鲁棒性和经济效益。 适合人群：从事电力系统研究、虚拟电厂调度、智能电网优化的研究人员和技术人员。 使用场景及目标：适用于需要处理光伏出力和负荷预测不确定性的虚拟电厂调度场景，旨在提高调度方案的稳定性和经济性，降低运营风险。 其他说明：文中提供了详细的MATLAB代码片段，帮助读者理解和实现具体的算法步骤。此外，还给出了若干实用技巧和注意事项，如场景削减时的概率归一化、储能SOC约束的跨场景耦合等。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值