**************************
VeryGood!
**************************
解：
mian()
{printf(“**************************”)
printf(“\n”)
printf(“VeryGood!\n”)
printf(“\n”)
printf(“**************************”)
}
1.6 编写一个程序，输入a、b、c三个值，输出其中最大值。
解：
mian()
{int a,b,c,max
printf(“请输入三个数a,b,c:\n”)
scanf(“%d,%d,%d”,&a,&b,&c)
max=a
if(max<b)
max=b
if(max<c)
max=c
printf(“最大数为：“%d”,max)
}
第三章
3.6写出以下程序运行的结果。
main()
{char c1=’a’,c2=’b’,c3=’c’,c4=’\101’,c5=’\116’
printf(“a%cb%c\tc%c\tabc\n”,c1,c2,c3)
printf(“\t\b%c %c”,c4,c5)
}
解：
aaㄩbbㄩㄩㄩccㄩㄩㄩㄩㄩㄩabc
AㄩN
3.7 要将"China"译成密码，译码规律是：用原来字母后面的第4个字母代替原来的字母．例如，字母"A"后面第4个字母是"E"．"E"代替"A"。因此，"China"应译为"Glmre"。请编一程序，用赋初值的方法使cl、c2、c3、c4、c5五个变量的值分别为，’C’、’h’、’i’、’n’、’a’，经过运算，使c1、c2、c3、c4、c5分别变为’G’、’l’、’m’、’r’、’e’，并输出。
解：
#include <stdio.h>
main()
{ char c1=’C’,c2=’h’,c3=’i’,c4=’n’,c5=’a’
c1+=4
c2+=4
c3+=4
c4+=4
c5+=4
printf("密码是%c%c%c%c%c\n",c1,c2,c3,c4,c5)
}
运行结果：
密码是Glmre
3.9求下面算术表达式的值。
（1）x+a%3*(int)(x+y)%2/4
设x=2.5,a=7,y=4.7
（2）(float)(a+b)/2+(int)x%(int)y
设a=2,b=3,x=3.5,y=2.5
（1）2.5
（2）3.5
3.10写出程序运行的结果。
main()
{int i,j,m,n
i=8
j=10
m=++i
n=j++
printf(“%d,%d,%d,%d”,i,j,m,n)
}
解：
9，11，9，10
3.12 写出下面表达式运算后a的值，设原来a=12。设a和n都已定义为整型变量。
（1）a+=a （2） a-=2 （3） a*=2+3 （4）a/=a+a
（5） a%=(n%=2)，n的值等于5
（6）a+=a-=a*=a
解：
(1) 24(2) 10(3) 60 (4) 0(5) 0 (6) 0
第四章
4．4若a=3，b=4，c=5，x=1.2，y=2.4，z=-3.6，u=51274，n=128765，c1=’a’，c2=’b’。想得到以下输出格式和结果，请写出程序（包括定义变量类型和设计输出）。
a=_3_ _b=_4_ _c=_5
x=1.200000,y=2.400000,z=-3.600000
x+y=_3.600_ _y+z=-1.20_ _z+x=-2.40
c1=ˊaˊ_or_97(ASCII)
c2=ˊbˊ_or_98(ASCII)
main()
{int a=3,b=4,c=5
long int u=51274,n=128765
float x=1.2,y=2.4,z=3.6
char c1=’a’,c2=’b’
printf("a=%2d b=%2d c=%2d\n",a,b,c)
printf("x=%f,y=%f,z=%f\n",x,y,z)
printf("x+y=%5.2f y+z=%5.2f z+x=%5.2f\n",x+y,y+z,z+x)
printf("u=%6ld n=%9ld\n",u,n)
printf("%s %s %d%s\n","c1=’a’","or",c1,"(ASCII)")
printf("%s %s %d%s\n","c2=’a’","or",c2,"(ASCII)")
}
4．7用scanf下面的函数输入数据，使a=3，b=7，x=8.5，y=71.82，c1=ˊAˊ，c2=ˊaˊ，问在键盘上如何输入？
main()
{
int a，b；float x，y；char c1c2；
scanf（"a=%d_b=%d"，&a，&b）；
scanf（"_x=%f_y=%e",&x，&y）；
scanf（"_c1=%c_c2=%c"，&c1，&c2）；
}
a=3_b=7
_x=8.5_y=71.82
_c1=A_c2=a
4．8设圆半径r=1.5，圆柱高h=3，求圆周长、圆面积、圆球表面积、圆球体积、圆柱体积。用scanf输入数据，输出计算结果，输出时要求文字说明，取小数点后两位数字。请编程序。
main()
{float r,h,C1,Sa,Sb,Va,Vb
scanf("%f,%f",&r,&h)
C1=2*3.14*r
Sa=3.14*r*r
Sb=4*Sa
Va=4*3.14*r*r*r/3
Vb=Sa*h
printf("C1=%.2f\n",C1)
printf("Sa=%.2f\nSb=%.2f\nVa=%.2f\nVb=%.2f\n",Sa,Sb,Va,Vb)
}
4．9输入一个华氏温度，要求输出摄氏温度。公式为
c=5(F-32)/9
输出要求有文字说明，取位2小数。
main()
{float F,c
scanf("%f",&F)
c=5*(F-32)/9
printf("c=%.2f",c)
}
4．10编程序，用getchar函数读入两个字符给c1、c2，然后分别用函数和函数输出这两个字符。并思考以下问题：（1）变量c1、c2应定义为字符型或整形？抑二者皆可？（2）要求输出c1和c2值的ASCII码，应如何处理？用putchar函数还是printf函数？（3）整形变量与字符变量是否在任何情况下都可以互相代替？如：
char c1，c2；
与
int c1，c2；
是否无条件的等价？
#include"stdio.h"
main()
{char c1,c2
c1=getchar()c2=getchar()
putchar(c1)putchar(’\n’)putchar(c2)putchar(’\n’)
}
#include"stdio.h"
main()
{char c1,c2
c1=getchar()c2=getchar()
printf("c1=%d c2=%d\n",c1,c2)
printf("c1=%c c2=%c\n",c1,c2)
}
第五章
5.1 什么是算术运算？什么是关系运算？什么是逻辑运算？
解：略。
5.2 C语言中如何表示“真”和“假”？系统如何判断一个量的“真”和“假”？
解：设有一个逻辑表达式，若其结果为“真”，则以1表示；若其结果为“假”，则以0表示。但是判断一个逻辑量的值时，以0代表“真”，以非0代表“假”。例如3&&5的值为“真”，系统给出3&&5的值为1。
5.3 写出下面各逻辑表达式的值。设a=3，b=4，c=5。
(1) a+b>c&&b==c
(2) a||b+c&&b-c
(3) !(a>b)&&!c||1
(4) !(x=a)&&(y=b)&&0
(5) !(a+b)+c-1&&b+c/2
解：
(1) 0
(2) 1
(3) 1
(4) 0
(5)1
5.4 有3个整数a、b、c，由键盘输入，输出其中最大的数。
解：
方法一
#include <stdio.h>
main()
{ int a,b,c
printf("请输入3个整数：")
scanf("%d,%d,%d",&a,&b,&c)
if(a<b)
if(b<c) printf("max=%d\n",c)
else printf("max=%d\n",b)
else if(a<c) printf("max=%d\n",c)
else printf("max=%d\n",a)
}
方法二：使用条件表达式，可以使程序更加简明、清晰。
程序如下：
#include <stdio.h>
main()
{ int a,b,c,temp,max
printf("请输入3个整数：")
scanf("%d,%d,%d",&a,&b,&c)
temp=(a>b)?a:b /* 将a和b中的大者存人temp中 */
max=(temp>c)?temp:c/* 将a和b中的大者与c比较，取最大者*/
printf("3个整数的最大数是%d\n”,max)
}
5.5 有一函数：
写一程序，输入x值，输出y值。
解：
#include <stdio.h>
main()
{int x,y
printf("输入x：")
scanf("%d",&x)
if(x<1) /* x<1 */
{ y=x
printf("x=%3d, y=x=%d\n",x,y)
}
else if (x<10) /* 1≤x-10 */
{ y=2*x-1
printf("x=%3d, y=2*x-1=%d\n",x,y)
}
else /* x≥10 */
{ y=3*x-11
printf("x=%3d, y=3*x-11=%d\n",x,y)
}
}
5.6 给一个百分制成绩，要求输出等级’A’、’B’、’C’、’D’、’E’。90分以上为’A’，80~90分为’B’，70~79分为’C’，60分以下为’D’。
解：
程序如下：
#include <stdio.h>
main()
{ float score
char grade
printf("请输入学生成绩：")
scanf("%f",&score)
while(score>100||(score<0)
{ printf("\n输入有误,请重新输入：")
scanf("%f",&score)
}
switch((int)(score/10))
{ case 10:
case 9: grade=’A’break
case 8: grade=’B’break
case 7: grade=’C’break
case 6: grade=’D’break
case 5:
case 4:
case 3:
case 2:
case 1:
case 0: grade=’E’
}
printf("成绩是%5.1f，相应的等级是%c。\n",score,grade)
}
说明：对输入的数据进行检查，如小于0或大于100，要求重新输入。(int)(score/10)的作用是将 (score/10) 的值进行强制类型转换，得到一个整型值。
5.7 给定一个不多于5位的正整数，要求：① 求它是几位数；② 分别打印出每一位数字；③ 按逆序打印出各位数字。例如原数为321，应输出123。
解：
#include <stdio.h>
main()
{ long int num
int indiv,ten,hundred,thousand,ten_thousand,place
/*分别代表个位、十位、百位、千位、万位和位数*/
printf("请输入一个整数（0~99999）：")
scanf("%ld",&num)
if (num>9999) place=5
else if(num>999) place=4
else if(num>99) place=3
else if(num>9) place=2
else place=1
printf("place =%d\n", place)
ten_thousand=num/10000
thousand=num/1000%10
hundred=num/100%10
ten=num%100/10
indiv=num%10
switch(place)
{ case 5: printf("%d,%d,%d,%d,%d",ten_thousand,thousand,hundred,ten,indiv)
printf("\n反序数字为；")
printf("%d%d%d%d%d\n",indiv,ten,hundred,thousand,ten_thousand)
break
case 4: printf("%d,%d,%d,%d",thousand,hundred,ten,indiv)
printf("\n反序数字为：")
printf("%d%d%d%d\n",indiv,ten,hundred,thousand)
break
case 3: printf("%d,%d,%d",hundred,ten,indiv)
printf("\n反序数字为：")
printf("%d%d%d\n",indiv,ten,hundred)
break
case 2: printf("%d,%d",ten,indiv)
printf("\n反序数字为：")
printf("%d%d\n",indiv,ten)
break
case 1: printf("%d",indiv)
printf("\n反序数字为：")
printf("%d\n",indiv)
break
}
}
5.8 企业发放的奖金根据利润提成。利润I低于或等于10万元时，奖金可提成10% ；利润高于10万元，低于20万元（100000<I≤200000）时，其中10万元按10%提成，高于10万元的部分，可提成7.5% ；200000<I≤400000时，其中20万元仍按上述办法提成（下同），高于20万元的部分按5%提成；400000<I≤600000时，高于40万元的部分按3%提成；600000〈I≤1000000时，高于60万的部分按1.5%提成；I>1000000时，超过100万元的部分按1%提成。从键盘输入当月利润I，求应发放奖金总数。要求：(1)用if语句编程序；(2)用switch语句编程序。
解：计算利润时，要特别注意不同利润的不同提成比例。例如，利润为15万元，其中有10万元按10%的比例提成，另外5万元则按7.5%提成。
(1) 用if语句编程序。
#include <stdio.h>
main()
{ long i
float bonus,bon1,bon2,bon4,bon6,bon10
bon1=100000*0.1 /*利润为10万元时的奖金*/
bon2=bon1+100000*0.075/*利润为20万元时的奖金*/
bon4=bon2+200000*0.05 /*利润为40万元时的奖金*/
bon6=bon4+200000*0.03 /*利润为60万元时的奖金*/
bon10=bon6+400000*0.015 /*利润为100万元时的奖金*/
printf("请输入利润i：")
scanf("%ld",&i)
if(i<=100000)
bonus=i*0.1 /*利润在10万元以内按0.1提成奖金*/
else if(i<=200000)
bonus=bon1+(i-100000)*0.075 /*利润在10万至20万元时的奖金*/
else if(i<=400000)
bonus=bon2+(i-200000)*0.05/*利润在20万至40万元时的奖金*/
else if(i<=600000)
bonus=bon4+(i-400000)*0.03/*利润在40万至60万元时的奖金*/
else if(i<=1000000)
bonus=bon6+(i-600000)*0.015 /*利润在60万至100万元时的奖金*/
else
bonus=bon10+(i-1000000)*0.01 /*利润在100万元以上时的奖金*/
printf(”奖金是%10.2f\n",bonus)
}
此题的关键在于正确写出每一区间的奖金计算公式。例如利润在10万元至20万时，奖金应由两部分组成：①利润为10万元时应得的奖金。即100000ⅹ0.1；②10万元以上部分应得的奖金。即(num-100000)ⅹ0.075。同理，20万～40万这个区间的奖金也应由两部分组成：①利润为20万元时应得的奖金，即100000ⅹ0.1ⅹ10万ⅹ0.075；②20万元以上部分应得的奖金，即(num-200000)ⅹ0.05。程序中先把10万、20万、40万、60万、100万各关键点的奖金计算出来，即bon1、bon2、bon4、bon6、hon10；然后再加上各区间附加部分的奖金。
(2) 用switch语句编程序。
输入利润i，确定相应的提成等级branch
根据branch确定奖金值
0奖金=i*0.1
1奖金=bon1+(i-105)*0.075
2奖金=bon2+(i-2*105)*0.05
3
4奖金=bon4+(i-4*105)*0.03
5
6奖金=bon6+(i-6*105)*0.015
7
8
9
10奖金=bon10+(i-106)*0.01
输出奖金
#include <stdio.h>
main()
{ long i
float bonus, bon1, bon2, bon4, bon6, bon10
int c
bon1=100000*0.1
bon2=bon1+100000*0.075
bon4=bon2+200000*0.05
bon6=bon4+200000*0.03
bon10=bon6+400000*0.015
printf("请输入利润i：")
scanf("%ld",&i)
c=i/100000
if(c>10) c=10
switch(c)
{ case 0: bonus=1*0.1break
case 1: bonus=bon1+(i-100000)*0.075break
case 2 :
case 3: bonus=bon2+(i-200000)*0.05break
case 4:
case 5: bonus=bon4+(i-400000)*0.03break
case 6:
case 7:
case 8:
case 9: bonus=bon6+(i-600000)*0.015break
case 10: bonus=bon10+(i-1000000)*0.01
}
printf("奖金是%10.2f",bonus)
}
5.9 输入4个整数，要求按由大到小的顺序输出。
解：此题采用依次比较的方法排出其大小顺序。在学习了循环和数组以后，可以有更多的排序方法。
#include <stdio.h>
main()
{ int t,a,b,c,d
printf("请输入4个整数：")
scanf("%d,%d,%d,%d",&a,&b,&c,&d)
printf("\n a=%d,b=%d,c=%d,d=%d\n",a,b,c,d)
if(a>b) {t=aa=bb=t}
if(a>c) {t=aa=cc=t}
if(a>d) {t=aa=dd=t}
if(b>c) {t=ab=cc=t}
if(b>d) {t=bb=dd=t}
if(c>d) {t=cc=dd=t}
printf("排序结果如下：\n")
printf("%d, %d, %d, %d\n",a,b,c,d)
}
5.10 有4个圆塔，圆心分别为(2，2)、(-2，2)、(2，-2)、(-2，-2)，圆半径为1。这4个塔的高度分别为10m。塔以外无建筑物。今输入任一点的坐标，求该点的建筑高度(塔外的高度为零)。
程序如下：
#include <stdio.h>
main()
{ int h=10
float x1=2,y1=2,x2=-2,y2=2,x3=-2,y3=-2,x4=-2,y4=-2,x,y,d1,d2,d3,d4
printf("请输入一个点(x,y)：")
scanf("%f,%f",&x,&y)
d1=(x-x1)*(x-x1)+(y-y1)*(y-y1) /*求该点到各中心点的距离*/
d2=(x-x2)*(x-x2)+(y+y2)*(y+y2)
d3=(x+x3)*(x+x3)+(y-y3)*(y-y3)
d4=(x+x4)*(x-x4)*(y+y4)*(y+y4)
if(d1>1&&d2>1&&d3>1&&d4>1) h=0 /*判断该点是否在塔外*/
printf("该点高度为%d\n",h)
}
第六章
第六章 循环控制
6.1输入两个正整数m和n，求其最大公约数和最小公倍数。
main()
{long m,n,i=1,j,s
scanf("%ld,%ld",&m,&n)
for(i<=m&&i<=ni++)
{if(m%i==0&&n%i==0) s=i}
if(m>=n) j=m
else j=n
for(!(j%m==0&&j%n==0)j++)
printf("s=%ld,j=%ld\n",s,j)
}
6.2输入一行字符，分别统计出其中英文字母、空格、数字和其他字符的个数。
#include"stdio.h"
main()
{char cint i=0,j=0,k=0,l=0
while((c=getchar())!=’\n’)
{if(c>=65&&c<=90||c>=97&&c<=122) i++
else if(c>=48&&c<=57) j++
else if(c==32) k++
else l++}
printf("i=%d,j=%d,k=%d,l=%d\n",i,j,k,l)
}
6.3求Sn=a+aa+aaa+…+aa…aaa（有n个a）之值，其中a是一个数字。例如：2+22+222+2222+22222（n=5），n由键盘输入。
#include"math.h"
main()
{int n,sum=0,i=1,s=2
scanf("%d",&n)
while(i<=n)
{sum=sum+ss=s+2*pow(10,i)
i++}
printf("sum=%d\n",sum)
}
6.4 求 ，（即求1!+2!+3!+4!+5!+…+20!）
main()
{int n,i=1long sum=0,s=1
scanf("%d",&n)
while(i<=n) {s=s*isum=sum+si++}
printf("sum=%ld\n",sum)
}
6.5 求
main()
{double i=1,j=1,k=1,s1=0,s2=0,s3=0,sum
for(i<=100i++) s1=s1+i
for(j<=50j++) s2=s2+j*j
for(k<=10k++) s3=s3+1/k
sum=s1+s2+s3
printf("sum=%f\n",sum)
}
6.6打印出所有"水仙花数"，所谓"水仙花数"是指一个三位数，其各位数字立方和等于该本身。例如：153是一个水仙花数，因为153=1^3+5^3+3^3。
#include"math.h"
main()
{int x=100,a,b,c
while(x>=100&&x<1000) {a=0.01*xb=10*(0.01*x-a)c=x-100*a-10*b
if(x==(pow(a,3)+pow(b,3)+pow(c,3))) printf("%5d",x)x++}
}
6.7一个数如果恰好等于它的因子之和，这个数就称为"完数"。例如，6的因子为1、2、3，而6=1+2+3，因此6是"完数"。编程序找出1000之内的所有完数，并按下面格式输出其因子：
6 its factors are 1、2、3
main()
{int m,i,j,s
for(m=6m<10000m++)
{s=1
for(i=2i<mi++)
if(m%i==0) s=s+i
if(m-s==0)
{printf("%5d its fastors are 1 ",m)for(j=2j<mj++) if(m%j==0)
printf("%d ",j)printf("\n")}
}
}
或
main()
{int m,i,j,s
for(m=6m<1000m++)
{s=m-1
for(i=2i<mi++)
if(m%i==0) s=s-i
if(s==0)
{printf("%5d its fastors are 1 ",m)for(j=2j<mj++) if(m%j==0)
printf("%d ",j)printf("\n")}
}
}
6.8有一分数序列：
求出这个数列的前20项之和。
main()
{int i=1,ndouble t,x=1,y=2,s,sum=0
scanf("%ld",&n)
while(i<=n) {s=y/xsum=sum+st=yy=y+xx=ti++}
printf("%f\n",sum)
}
6.9一球从100米高度自由下落，每次落地后返回原高度的一半，再落下。求它在第10次落地时共经过多少米？第10次反d多高？
main()
{int i,ndouble h=100,s=100
scanf("%d",&n)
for(i=1i<=ni++)
{h*=0.5if(i==1) continues=2*h+s}
printf("h=%f,s=%f\n",h,s)
}
6.10猴子吃桃问题。猴子第一天摘下若干个桃子，当即吃了一半，还不过瘾，又多吃了一个。第二天早上又将剩下的桃子吃掉一半，又多吃一个。以后每天早上都吃了前一天剩下的一半零一个。到第10天早上想再吃时，见只剩下一个桃子了。求第一天共摘多少桃子。
main()
{int i=1,sum=0
for(i<=10sum=2*sum+1,i++)
printf("sum=%d\n",sum)
}
6.11用迭代法求 。求平方根的迭代公式为：
要求前后两次求出的得差的绝对值少于0.00001。
#include"math.h"
main()
{float x0,x1,a
scanf("%f",&a)
x1=a/2
do
{x0=x1x1=(x0+a/x0)/2}
while(fabs(x0-x1)>=0.00001)
printf("%.3f\n",x1)
}
6.12 用牛顿迭代法求方程在1.5附近的根。
main()
{double x,yx=1.5
do{y=2*x*x*x-4*x*x+3*x-6
x=x-y/(6*x*x-8*x+3)}
while(y!=0)
printf("x=%.3f\n",x)
}
6.13用二分法求方程在(-10，10)之间的根
main()
{double x1,x2,y1,y2x1=-10x2=10
do{y1=2*x1*x1*x1-4*x1*x1+3*x1-6
x1=x1-y1/(6*x1*x1-8*x1+3)}
while(y1!=0)
do
{y2=2*x2*x2*x2-4*x2*x2+3*x2-6
x2=x2-y2/(6*x2*x2-8*x2+3)}
while(y2!=0)
printf("x1=%.3f,x2=%.3f\n",x1,x2)
}
6.14打印以下图案
*
* * *
* * * * *
* * * * * * *
* * * * *
* * *
*
MATLAB程序设计与应用（第二版）实验参考答案%实验一 MATLAB运算基础
%第一题
%（1）
z1=2*sin(85*pi/180)/(1+exp(2))
%（2）
x=[2,1+2i-0.45,5]
z2=0.5*log(x+sqrt(1+x.^2))
%(3)
a=-3.0:0.1:3.0
z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)
%(4)
t=0:0.5:2.5
z4=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3)
%第二题
A=[12 34 -434 7 873 65 7]
B=[1 3 -12 0 33 -2 7]
A+6*B
A-B+eye(size(A))
A*B
A.*B
A^3
A.^3
A/B
B\A
[A,B]
[A([1,3],:)B^2]
%第三题
A=[1 2 3 4 56 7 8 9 1011 12 13 14 1516 17 18 19 2021 22 23 24 25]
B=[3 0 1617 -6 90 23 -49 7 04 13 11]
C=A*B
F=size(C)
D=C(F(1)-2:F(1),F(2)-1:F(2))
whos
%第四题
%(1):
A=100:999
B=rem(A,21)
C=length(find(B==0))
%(2):
A='lsdhKSDLKklsdkl'
k=find(A>='A'&A<='Z')
A(k)=[]
%实验二 MATLAB矩阵分析与处理
%第一题
E=eye(3)
R=rand(3,2)
O=zeros(2,3)
S=diag([2,3])
A=[E,RO,S]
A^2
B=[E,(R+R*S)O,S^2]
%第二题
H=hilb(5)
P=pascal(5)
Hh=det(H)
Hp=det(P)
Th=cond(H)
Tp=cond(P)
%第三题:
A=fix(10*rand(5))
H=det(A)
Trace=trace(A)
Rank=rank(A)
Norm=norm(A)
%第四题:
A=[-29,6,1820,5,12-8,8,5]
[V,D]=eig(A)
%数学意义略
%第五题方法一:
%(1):
A=[1/2,1/3,1/41/3,1/4,1/51/4,1/5,1/6]
b=[0.95,0.67,0.52]'
x=inv(A)*b
%(2):
B=[0.95,0.67,0.53]'
x=inv(A)*B
%(3):
cond(A)
%第五题方法二:
A=hilb(4)
A(:,1)=[]
A(4,:)=[]
B=[0.95,0.67,0.52]'
X=inv(A)*B
B1=[0.95,0.67,0.53]'
X1=inv(A)*B1
N=cond(B)
N1=cond(B1)
Na=cond(A) %矩阵A为病态矩阵
%第六题
A=[1,4,916,25,3649,64,81]
B=sqrtm(A)
C=sqrt(A) %sqrtm函数是以矩阵为单位进行计算，sqrt函数是以矩阵中的元素进行计算
%实验三 选择程序结构设计
%第一题程序一
x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0]
y=[]%建立存放所有y值的矩阵
for x0=x
if x0<0&x0~=-3
y=[y,x0*x0+x0-6]
elseif x0>=0&x0<5&x0~=2&x0~=3
y=[y,x0*x0-5*x0+6]
else
y=[y,x0*x0-x0-1]
end
end
x %输出所有x
y %输出所有y
%第一题程序二
x=[-5,-3,1,2,2.5,3,5]
y=[]
for a=1:7
if x(a)<0&x(a)~=-3
y=[y,(x(a))^2+x(a)-6]
elseif x(a)>=0&x(a)<5&x(a)~=2&x(a)~=3
y=[y,(x(a))^2-5*x(a)+6]
else
y=[y,x(a)*x(a)-x(a)-1]
end
end
%第二题程序一:
x=input('请输入一个百分制成绩：')
if x>100|x<0
disp('您输入的成绩不是百分制成绩，请重新输入。')
else
if x<=100&x>=90
disp('A')
elseif x<=89&x>=80
disp('B')
elseif x<=79&x>=70
disp('C')
elseif x<=69&x>60
disp('D')
else
disp('E')
end
end
%第二题程序二：
s=input('请输入一个成绩（0分到100分之间）：') %s用于存放成绩
while 1 %判断输入成绩的合理性
if s<0|s>100
disp('输入的成绩需在0到100之间，请重新输入：')
s=input('请输入一个成绩（0分到100分之间）：')
else
break
end
end
switch fix(s/10) %对成绩做出等级判断
case {9,10}
disp('A')
case 8
disp('B')
case 7
disp('C')
case 6
disp('D')
otherwise
disp('E')
end
%第三题
n=input('请输入员工工号：')
h=input('该员工工作时数是：')
if h>120
x=(h-120)*84*(1+0.15)+120*84
elseif h<60
x=h*84-700
else
x=h*84
end
disp([num2str(n),'号员工','的应发工资为',num2str(x)])
%第四题（还可以用switch语句实现）
a=fix(10+(99-10)*rand(1,2))%产生两个随机整数
x=a(1)
y=a(2)
t=input('请输入运算符号：','s')
if t=='+'
z=x+y
elseif t=='-'
z=x-y
elseif t=='*'
z=x*y
elseif t=='/'
z=x/y
end
disp([num2str(x),t,num2str(y),'=',num2str(z)]) %输出运算结果
%第五题
a=rand(5,6) %产生5x6的随机矩阵
n=input('请输入您要输出矩阵的第几行：')
if n>5
disp('超出了矩阵的行数，矩阵的最后一行为：')
a(5,:)
else
disp(['矩阵的第',num2str(n),'行为：'])
a(n,:)
end
%实验四 循环结构程序设计
%第一题程序一
s=0
n=input('n=?')
for i=1:n
s=s+1/i/i
end
PI=sqrt(6*s)
pi
%第一题程序二
n=input('n=?')
a=1:n
b=1./a.^2
PI=sqrt(6*sum(b))
pi
%第二题
y=0
n=1
while(y<3)
y=y+1/(2*n-1)
n=n+1
end
y=y-1/(2*(n-1)-1)
n=n-2
%第三题
a=input('a=?')
b=input('b=?')
Xn=1
Xn1=a/(b+Xn)
n=0
while abs(Xn1-Xn)>1e-5
Xn=Xn1
Xn1=a/(b+Xn)
n=n+1
if n==500
break
end
end
n
Xn1
r1=(-b+sqrt(b*b+4*a))/2
r2=(-b-sqrt(b*b+4*a))/2
%第四题
for i=1:100
if i==1
f(i)=1
elseif i==2
f(i)=0
elseif i==3
f(i)=1
else
f(i)=f(i-1)-2*f(i-2)+f(i-3)
end
end
max(f)
min(f)
sum(f)
length(find(f>0))
length(find(f==0))
length(find(f<0))
%第五题：
s=0n=0
for i=2:49
b=i*(i+1)-1
m=fix(sqrt(b))
for j=2:m
if rem(b,j)==0
break
end
end
if j==m
n=n+1
s=s+b
end
end
n
s
%实验五 函数文件
%第一题
function y=mat1(x)%建立函数文件mat1.m
y=[exp(x),log(x),sin(x),cos(x)]
%在命令窗口调用上述函数文件：
y=mat1(1+i)
%第二题程序一
function [a,b,N,M]=shiyanwu2(m,n,t)
A=[m*cos(t*pi/180),-m,-sin(t*pi/180),0m*sin(t*pi/180),0,cos(t*pi/180),00,n,-sin(t*pi/180),00,0,-cos(t*pi/180),1]
B=[0,9.8*m,0,9.8*n]
C=inv(A)*B'
a=C(1)
b=C(2)
N=C(3)
M=C(4)
%在命令窗口调用该函数文件:
m1=input('m1=')
m2=input('m2=')
theta=input('theta=')
[a1,a2,N1,N2]=shiyanwu2(m1,m2,theta)
%第二题程序二
function X=mat2(m1,m2,t)
g=9.8
A=[m1*cos(t*pi/180),-m1,-sin(t*pi/180),0m1*sin(t*pi/180),0,cos(t*pi/180),00,m2,-sin(t*pi/180),00,0,-cos(t*pi/180),1]
B=[0m1*g0m2*g]
X=inv(A)*B
%在命令窗口调用该函数文件：
X=mat2(1,1,60)
%第三题
function flag=mat3(x)
flag=1
for i=2:sqrt(x)
if rem(x,i)==0
flag=0
break
end
end
%在命令窗口调用该函数文件：
for i=10:99
j=10*rem(i,10)+fix(i/10)
if mat3(i)&mat3(j)
disp(i)
end
end
%第四题
function y=fx(x)
y=1./((x-2).^2+0.1)+1./((x-3).^4+0.01)
%在命令窗口调用该函数文件：
y=fx(2)
a=[1,23,4]
y=fx(a)
%第五题
%(1)
function f1=mat5(n)
f1=n+10*log(n*n+5)
%在命令窗口中调用该函数文件：
y=mat5(40)/(mat5(30)+mat5(20))
%(2)方法一
function f2=mat6(n)
f2=0
for i=1:n
f2=f2+i*(i+1)
end
%在命令窗口中调用该函数文件如：
y=mat6(40)/(mat6(30)+mat6(20))
%(2)方法二
function f2=mat7(n)
i=1:n
m=i.*(i+1)
f2=sum(m)
end
%在命令窗口中调用该函数文件如：
y=mat7(40)/(mat7(30)+mat7(20))
%实验六 高层绘图 *** 作
%第一题：
x=linspace(0,2*pi,101)
y=(0.5+3*sin(x)./(1+x.^2)).*cos(x)
plot(x,y)
%第二题：
%（1）
x=linspace(-2*pi,2*pi,100)
y1=x.^2
y2=cos(2*x)
y3=y1.*y2
plot(x,y1,'b-',x,y2,'r:',x,y3,'y--')
text(4,16,'\leftarrow y1=x^2')
text(6*pi/4,-1,'\downarrow y2=cos(2*x)')
text(-1.5*pi,-2.25*pi*pi,'\uparrow y3=y1*y2')
%（2）
x=linspace(-2*pi,2*pi,100)
y1=x.^2
y2=cos(2*x)
y3=y1.*y2
subplot(1,3,1)%分区
plot(x,y1)
title('y1=x^2')%设置标题
subplot(1,3,2)
plot(x,y2)
title('y2=cos(2*x)')
subplot(1,3,3)
plot(x,y3)
title('y3=x^2*cos(2*x)')
%（3）
x=linspace(-2*pi,2*pi,20)
y1=x.^2
subplot(2,2,1)%分区
bar(x,y1)
title('y1=x^2的条形图')%设置标题
subplot(2,2,2)
stairs(x,y1)
title('y1=x^2的阶梯图')
subplot(2,2,3)
stem(x,y1)
title('y1=x^2的杆图')
subplot(2,2,4)
fill(x,y1,'r')%如果少了'r'则会出错
title('y1=x^2的填充图')
%其他的函数照样做。
%第三题
x=-5:0.01:5
y=[]%起始设y为空向量
for x0=x
if x0<=0 %不能写成x0=<0
y=[y,(x0+sqrt(pi))/exp(2)] %将x对应的函数值放到y中
else
y=[y,0.5*log(x0+sqrt(1+x0^2))]
end
end
plot(x,y)
%第四题：
a=input('a=')
b=input('b=')
n=input('n=')
t=-2*pi:0.01:2*pi
r=a*sin(b+n*t)
polar(t,r)
%第五题
x=linspace(-5,5,21)
y=linspace(0,10,31)
[x,y]=meshgrid(x,y)%在[-5,5]*[0,10]的范围内生成网格坐标
z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2)/4)
subplot(2,1,1)
surf(x,y,z)
subplot(2,1,2)
contour3(x,y,z,50)%其中50为高度的等级数，越大越密
%第六题
ezsurf('cos(s)*cos(t)','cos(s)*sin(t)','sin(s)',[0,0.5*pi,0,1.5*pi])%利用ezsurf隐函数
shading interp %进行插值着色处理
%实验七 低层绘图 *** 作
%第一题
h=figure('MenuBar','figure','color','r','WindowButtonDownFcn','disp(''Left Button Pressed'')')
%第二题
x=-2:0.01:2
y=x.^2.*exp(2*x)
h=line(x,y)
set(h,'color','r','linestyle',':','linewidth',2)
text(1,exp(2),'y=x^2*exp(2*x)')
%第三题
t=0:0.00001:0.001
[t,x]=meshgrid(t)
v=10*exp(-0.01*x).*sin(2000*pi*t-0.2*x+pi)
axes('view',[-37.5,30])
h=surface(t,x,v)
title('v=10*exp(-0.01*x).*sin(2000*pi*t-0.2*x+pi)')
xlabel(Ct'),ylabel('x'),zlabel('v')
%第四题
x=0:0.01:2*pi
y1=sin(x)
y2=cos(x)
y3=tan(x)
y4=cot(x)
subplot(2,2,1)
plot(x,y1)
subplot(2,2,2)
plot(x,y2)
subplot(2,2,3)
plot(x,y3)
subplot(2,2,4)
plot(x,y4)
%第五题
cylinder(5)
light('Position',[0,1,1])
material shiny
%实验八 数据处理与多项式运算
%第一题
%(1)
A=rand(1,30000)
b=mean(A)
std(A,0,2)
%(2)
max(A)
min(A)
%(3)
n=0
for i=1:30000
if A(i)>0.5
n=n+1
end
end
p=n/30000
%第二题
%(1)
A=45+51*rand(100,5)
[Y,U]=max(A)
[a,b]=min(A)
%(2)
m=mean(A)
s=std(A)
%(3)
sum(A,2)
[Y,U]=max(ans)
[a,b]=min(ans)
%(4)
[zcj,xsxh]=sort(ans)
%第三题
h=6:2:18
x=6.5:2:17.5
t1=[18,20,22,25,30,28,24]
t2=[15,19,24,28,34,32,30]
T1=spline(h,t1,x)
T2=spline(h,t2,x)
%第四题
x=1:0.1:101
y1=log10(x)
p=polyfit(x,y1,5)
y2=polyval(p,x)
plot(x,y1,':',x,y2,'-')
%第五题
%(1)
p1=[1,2,4,0,5]
p2=[1,2]
p3=[1,2,3]
p=p1+[0,conv(p2,p3)] %为使两向量大小相同，所以补0
%(2)
A=roots(p)
%(3)
A=[-1,1.2,-1.40.75,2,3.50,5,2.5]
polyval(p,A)
%(4)
polyvalm(p,A)
%实验十五 综合实验
%第一题
%定义函数文件
function zp=fs(z,n)
zp=0
z=0
for k=1:n
z=1+imag(z)-1.4*real(z)^2+i*real(z)*0.3
zp=[zp,z]
end
%在主窗口调用该函数：
n=30000
z=0
zp=fs(z,n)
plot(zp,'p','markersize',2,'color','r')
%第二题
quad('1/sqrt(2*pi)*exp(-x.*x/2)',0,1)
syms x
int(1/sqrt(2*pi)*exp(-x.*x/2),0,1)
%第三题
%定义函数文件
function dx=apollo(t,x)
mu=1/82.45
mustar=1-mu
r1=sqrt((x(1)+mu)^2+x(3)^2)
r2=sqrt((x(1)-mustar)^2+x(3)^2)
dx=[x(2)2*x(4)+x(1)-mustar*(x(1)+mu)/r1^3-mu*(x(1)-mustar)/r2^3
x(4)-2*x(2)+x(3)-mustar*x(3)/r1^3-mu*x(3)/r2^3]
%在主窗口调用该函数：
x0=[1.200-1.04935751]
options=odeset('reltol',1e-8)
%该命令的另一种写法是options=odesetoptions.reltol=1e-8
[t,y]=ode45(@apollo,[0,20],x0,options)
plot(y(:,1),y(:,3))%绘制x和y图像，也就是卫星的轨迹
title('Appollo卫星运动轨迹')
xlabel('X')
ylabel('Y')
%第四题
%这是一个力矩平衡问题，可列方程：500×1.5×cosα=400×2.0×cos(π/3-α)
%解得α＝26.8021°。
%程序一：迭代法解方程
x1=input('x1=?')
x=pi/3.0-acos(15.0*cos(x1)/16)
while abs(x-x1)>=10e-8
x1=x
x=pi/3.0-acos(15.0*cos(x1)/16)
if x>pi/3.0
disp('error')
end
end
y=x1*180/pi
%程序二：利用绘图功能
%分别绘制两个小孩所产生力矩随α变化的曲线，两曲线的交点即是跷跷板平衡时的α
alpha=linspace(0,pi/2,80000)
m1=500*1.5*cos(alpha) %第一个小孩产生的力矩
m2=400*2*cos(pi/3-alpha) %第二个小孩产生的力矩
k=find(abs(m1-m2)<1e-2) %找平衡点
alpha0=alpha(k)
m0=500*1.5*cos(alpha0)
plot(alpha,m1,alpha,m2,alpha0,m0,'kp')
alpha=alpha0*180/pi%化为角度
%第五题
clearclose
fplot('18-t^(2/3)',[0,20])grid onhold on
fplot('5+t+2*t^(2/3)',[0,20],'r')hold off
%发现t约为4
[t,f,h]=fsolve('18-x^(2/3)-5-x-2*x^(2/3)',4)
%求得t=4.6465
t=linspace(0,t,100)y=18-t.^(2/3)-5-t-2*t.^(2/3)
trapz(t,y)-20
%最大利润6.3232(百万元)
解析:学生成绩的三个等级由是否小于等于60分和是否小于等于79分这两个条件控制，因此可以将其分为两部分来设计框图，即判断框的选择问题，还应注意将全部学生的相应成绩等级全部显示. 答案:依题意设计以下程序框图： 点评:该题也可以先用判断条件“是否小于等于79分”先行判断，再根据“是否小于等于60分”，从而将成绩分成三个等级欢迎分享，转载请注明来源：内存溢出
微信扫一扫
支付宝扫一扫
评论列表（0条）