⑴ 求高手给出一维热传导方程解的图形的matlab代码
function main
%热传导方程的差分解法
lda=0.000015;%lad即为a^2
l=50;
h=0.05;%空间步长
alpha=0.5;
tao=alpha*h^2/lda; %时间步长
N=l/h;%空间步数
M=100;%时间步数,由于你条件没有限制,所以我随便取了个值
for i=1:N+1
u(1,i)=fai((i-1)*h);
end
for k=1:M
u(k,1)=g1(k*tao);
u(k,N+1)=g2(k*tao);
end
for k=1:M
for i=2:N
u(k+1,i)=alpha*u(k,i+1)+(1-2*alpha)*u(k,i)+alpha*u(k,i-1);
end
plot([0:h:l],u(k+1,:));
hold on; pause(0.05);
end
function u=fai(x)
u=100;
function u=g1(x)
u=0;
function u=g2(x)
u=0;
就这个拉,你试试看行不,不行的话在联系我吧!!
⑵ C语言编程(一维下料问题):某厂有一批长度为7.4m的钢管原材料(数量充分多)
先看有几种下料方法 按照你给闹逗的长度排列 可碧凯能的方悔弯唤法有(数字表示切割成零件的个数 我们排除掉那些留下余料还可以切割的情况 否则就会太复杂)
2 0 0
1 2 0
1 1 1
1 0 2
0 3 0
0 2 1
0 1 3
0 0 4
下面是程序
model:
sets:
method/1..8/:n;
parts/1..3/:;
link(method,parts):m;
endsets
data:
m=
2 0 0
1 2 0
1 1 1
1 0 2
0 3 0
0 2 1
0 1 3
0 0 4;
enddata
min=@sum(method:n);
@for(parts(j):@sum(method(i):m(i,j)*n(i))>=200);
@for(method:@gin(n));
end
⑶ C语言一维数组求平均值
C语言一维数组求平均值的源代码如下:
#include#define N 10
typedef int TYPE;
TYPE sum(TYPE a[],int n)
{
TYPE s=0;
int i;
for(i=0;ia[j]
{
temp=a[i];
a[i]=a[j];
a[j]=temp;
}
int main()
{
(3)一维下料问题代码扩展阅读
数组的初始化
1、在定义数组时对数组元素赋以初值,如:static int a[4]={1,2,3,4}。给一部分元素赋初值,如static int b[10]={1,2}。
2、对static数组不赋初值,系统会对所有数组元素自动赋以0值。
3、若对全部数组元素赋初值时,可以不指定数组长度,如:int c[]={1,2,3,4,5}。
4、若被定义数组长度与提供初值的个数不相同,则数组长度不能忽略。