第十次上机实验

任务1:
输入一个字符串和一个正整数x,将该字符串中的后x个字符复制到另一个字符串y中,再对y串的内容前后倒置后存入数组z中并输出。
要求:用指针访问数组元素、用函数getx(char *c1)实现复制、用函数getr(char *c2)实现倒置。
运行示例
Enter a string: abcABCD
Enter an integer: 4
The new string is DCBA
----------------------分割线----------------------

代码如下:

#include<stdio.h>
#include<string.h>
#define N 20
void getx(char *c1,char *m,int x)
{
	int i;
	int j=0;
	int n=strlen(c1);
	for(i=n-x;i<n;i++)
	{
		*(m+j)=*(c1+i);
		j++;
	}
	*(m+x)='\0';
}
void getr(char *c2)
{
	int i;
	char temp;
	int n=strlen(c2);
	for(i=1;i<=n/2;i++)
	{
		temp=*(c2+i-1);
		*(c2+i-1)=*(c2+n-i);
		*(c2+n-i)=temp;
	}
	*(c2+n)='\0';
}
void main()
{
	int x;
	int n,i;
	char a[N],b[N],z[N];
	printf("Enter a string: ");
	gets(a);
	printf("Enter an integer: ");
	scanf_s("%d",&x);
	getx(a,b,x);
	n=strlen(b);
	getr(b);
	for(i=0;i<n;i++)
	{
		z[i]=b[i];
	}
	z[n]='\0';
	puts(z);
}
运行效果图如下:


任务2:
定义一维整形数组,对数组分别进行“由大到小”和"由小到大"排序并输出。
要求:用函数和指针实现排序
----------------------分割线----------------------

代码如下:

#include<stdio.h>
#include<string.h>
#define N 100
void paix1(int *p,int n)  
{  
    int t;  
    for(int i=0;i<n;i++)
	{
        for(int j=0;j<n-i-1;j++)  
        {  
            if(*(p+j)<*(p+j+1))  
            {  
                t=*(p+j);  
                *(p+j)=*(p+j+1);  
                *(p+j+1)=t;  
            }  
        }  
	}
}  
void paix2(int *p,int n)  
{  
    int t;  
        for(int i=0;i<n;i++)  
        for(int j=0;j<n-i-1;j++)  
        {  
            if(*(p+j)>*(p+j+1))  
            {  
                t=*(p+j);  
                *(p+j)=*(p+j+1);  
                *(p+j+1)=t;  
            }  
        }  
}  
void main()
{
	int a[N],*p;
	int n,i;
	printf("请输入数组的元素:");
	scanf("%d",&n);
	for(i=0;i<n;i++)
	{
		scanf("%d",&a[i]);
	}
	p=a;
	paix1(a,n);
	printf("由大到小为:");
	for(i=0;i<n;i++)
	{
		printf("%d ",a[i]);
	}
	printf("\n");
	paix2(a,n);
	printf("有小到小为:");
		for(i=0;i<n;i++)
	{
		printf("%d ",a[i]);
	}
	printf("\n");
	getchar();
}

运行效果图为:


任务3:
输入字符串s,将字符放入d数组中,最后输出d中的字符串。
要求:用函数和指针实现
运行示例
输入字符串:abc123edf456gh
输出字符串:abcedfgh

代码如下:

#include<stdio.h>  
#include<string.h>  
void main()  
{  
    int i,j=0,n;  
    char a[100],b[100];  
    gets(a);  
    n=strlen(a);  
    for(i=0;i<n;i++)  
    {  
        if(a[i]>='A'&&a[i]<='Z'&&a[i]!='\0'||a[i]>='a'&&a[i]<='z'&&a[i]!='\0')  
        {     
            *(b+j)=*(a+i);  
            j++;  
        }  
        else  
            continue;  
    }  
    *(b+j)=0;  
    puts(b);  
    getchar();  
}  

运行效果图为:


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
共轭梯度法是求解线性方程组的一种有效方法,MATLAB中提供了相关的函数可以使用。这里我们将介绍如何手动编写共轭梯度法程序。 首先,我们需要明确共轭梯度法的基本思路。共轭梯度法是一种迭代算法,其核心思想是通过利用前一次迭代得到的信息来加速收敛。具体来说,对于一个对称正定的系数矩阵A和一个右端向量b,共轭梯度法的基本步骤如下: 1. 初始化迭代:取一个初始向量x0和残差向量r0=b-Ax0,设置迭代次数k=0。 2. 计算搜索方向:对于第k次迭代,计算搜索方向p_k=r_k,如果k=0,则p_0=r_0;否则,p_k=r_k+beta_k*p_{k-1},其中beta_k是一个系数,定义为beta_k=(r_k'*r_k)/(r_{k-1}'*r_{k-1})。 3. 计算步长:沿着搜索方向p_k移动一个步长alpha_k,使得x_{k+1}=x_k+alpha_k*p_k。 4. 更新残差:计算新的残差r_{k+1}=b-Ax_{k+1}。 5. 检查停止条件:如果满足停止条件(如残差的范数小于某个阈值),则停止迭代;否则,令k=k+1,返回步骤2。 基于以上思路,我们可以编写共轭梯度法的MATLAB程序。下面是一个简单的例子: ```matlab function [x, k] = conjgrad(A, b, x0, tol, maxiter) % 输入参数: % A: 系数矩阵 % b: 右端向量 % x0: 初始向量 % tol: 残差的阈值 % maxiter: 最大迭代次数 % 输出参数: % x: 迭代结束时的解向量 % k: 实际迭代次数 r = b - A * x0; % 初始残差 p = r; % 初始搜索方向 x = x0; % 初始解向量 k = 0; % 初始迭代次数 while norm(r) > tol && k < maxiter alpha = (r'*r) / (p' * A * p); % 计算步长 x = x + alpha * p; % 更新解向量 r_new = r - alpha * A * p; % 计算新的残差 beta = (r_new' * r_new) / (r' * r);% 计算beta p = r_new + beta * p; % 计算新的搜索方向 r = r_new; % 更新残差 k = k + 1; % 更新迭代次数 end end ``` 这个程序定义了一个名为conjgrad的函数,输入参数包括系数矩阵A、右端向量b、初始向量x0、残差阈值tol和最大迭代次数maxiter。输出参数包括迭代结束时的解向量x和实际迭代次数k。 在函数内部,我们首先计算了初始残差r、初始搜索方向p和初始解向量x。然后进入迭代循环,直到满足停止条件为止。在每次迭代中,我们先计算了步长alpha,然后更新了解向量x。接着计算了新的残差r_new和新的搜索方向p,最后更新了残差和迭代次数。在迭代结束后,函数返回最终的解向量和实际迭代次数。 使用这个函数非常简单,只需要提供系数矩阵A、右端向量b和初始向量x0,以及其他参数即可。例如: ```matlab A = [4 -1 0; -1 4 -1; 0 -1 4]; b = [10; 10; 10]; x0 = [0; 0; 0]; tol = 1e-6; maxiter = 100; [x, k] = conjgrad(A, b, x0, tol, maxiter); disp(['迭代次数:', num2str(k)]); disp(['解向量:', num2str(x')]); ``` 这个例子是求解线性方程组Ax=b,其中系数矩阵A是一个对称正定矩阵。我们设置了初始向量x0=[0; 0; 0],残差阈值tol=1e-6和最大迭代次数maxiter=100,然后调用conjgrad函数求解。最后输出了迭代次数和解向量。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值