【转载】matlab 学习点滴

晓夜寒风

在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数；

1. a=1:500;
   
2. b=a(find(rem(a,3)==2));
   
3. c=b(find(rem(b,5)==3));
   
4. d=c(find(rem(c,7)==2));
   
5. disp('在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数:')
   
6. disp(d)

复制代码

在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数:
/ A2 t& K4 M4 Y+ t$ E  23   128   233   338   443

晓夜寒风

matlab上标标注

1. subplot(1,1,1);
   + X. r0 s6 ^7 _% [
2. x=0:0.1:2*pi;
   
3. y1=sin(x);
   % R2 ^0 r9 G8 M% D% j
4. y2=exp(-x);
   8 t" u- l* `. y
5. plot(x,y1,'--*',x,y2,':o');
   5 [6 N* ^# k, u9 y
6. xlabel('x=0 to 2\pi');
   
7. ylabel('values of sin(x) and e^{-x}');
   : d" [5 r7 m9 r* B' ^
8. title('function plots fo sin(x) and e^{-x}');
   6 K7 K8 \" ^1 j+ O7 V
9. legend('sin(x)','e^{-x}');

复制代码

晓夜寒风

请问在matlab里面怎么把一个三段的分段函数表示出来呢

1. f(x)=x-1 (x>1);f(x)=0 (abs(x)<1);f(x)=x+1 (x<-1)
   $ r8 h- K' O* H  c& T* s
2. 1)x=[-10:1:10];
   
3. for i=1:21
   
4. f(i)=(x(i)-1)*(x(i)>1)+(x(i)+1)*(x(i)<-1);
   
5. end
   ' g$ Y2 o8 s2 K1 Q  s# x2 E
6. plot(x,f,'r*')
   
7. 2)x=[-10:10];
   + F) D, l2 Y. D' E( U
8. f=zeros(size(x));
   
9. a=x>1;
   
10. b=x<-1;
    3 U' d# E% E! v- j- x
11. f(a)=x(a)-1;
    
12. f(b)=x(b)+1;
    8 n$ k8 A7 \2 f9 \5 n+ G
13. 3)f(x)=(x-1).*(x>1)+0.*(abs(x)<1)+(x+1).*(x<-1)

复制代码

晓夜寒风

Q1: 对于x = [0.1200 0.2400 0.3600 0.5800 0.6600 0.7400 0.8100 0.8700 0.9100 1.0000]' 随机产生一个数a = rand
怎么样把大于a的第一个x（i）选出来？

1. x = [0.1200 0.2400 0.3600 0.5800 0.6Q600 0.7400 0.8100 0.8700 0.9100 1.0000]'
   
2. a=rand;
   
3. k=find(x>a);
   + m3 P; B  p0 P: {" k$ Q4 [, x
4. x(k(1))

复制代码

Q2: 如何从列矩阵(n*1)中分段截取数据组成新的矩阵?比如用列矩阵a的前k个数据组成一个新的列矩阵b?中间k个呢?

1. b=a(1:k);
   
2. b=a(l:l+k); (l+k<=n)
   . U8 P$ \6 v4 z! X

复制代码

Q3: [1，2，3，NaN，5，NaN，7，8，9，10]如何将NaN以外的8个数加起来，然后求这8个数的平均数。

1. mean_a=mean(a(~isnan(a)))

复制代码

Q4 matlab里面怎么将符号常数变为数值

1. syms x;x1=solve('x+1=2','x');
   
2. if double(x1)>1   
   ! B0 x3 U3 Q2 d7 H9 T* ~0 ^, B
3.   x=1
   
4. else   
   ) o/ r" p4 r2 _: D& ]" t3 I9 s& g! e
5.   x=0
   
6. end

复制代码

晓夜寒风

Q5: MATLAB中用pie命令画饼图

1. function PlotPie
   $ a2 D4 d& r' T% m( w+ u: L- Q
2. x=[1 2 3 4 5 6 7 8 9 10];
   5 y: m3 ~( |; }2 G; i
3. y=[40 50 60 70 80 90 100 110 120 130 ];
   7 W' i$ ]# ]$ W6 ^
4. str1={'A','B','C','D','E','F','G','H','I','J'};
   
5. por=y./sum(y);
   
6. for i=1:length(por)   str{i}=[str1{i},' ',num2str(100*por(i)),'%'];
   
7. end
   + i5 q3 e3 u6 f3 j! d2 D8 R  T
8. pie(x,y,str)

复制代码

晓夜寒风

Q6 非线性拟和y＝a*x^b+c的程序，求a,b,c.
# c+ {. ^5 Q! w2 b  I/ ^function FitPaperData
format long
- ^3 R, r2 e2 c& u# Y  |9 fE1=rand(89,1);
: x2 d4 x9 C: q1 Y  ~. V" x: KE2=rand(89,1)+3;
F=inline('p(1).*x.^p(2)+p(3)','p','x');
[p,r,j]=nlinfit(E1,E2,F,[1,2,3]);
P

晓夜寒风

Q7 已知向量A=[1 7 6 3 2 9 8];B=[1 6 3 8] ;如果我想把向量A中与向量B中相同的元素清除掉，例如最后得到C=[7 2 9]
应该怎么做
) S, w8 l7 R5 I7 c" ^, K如果各个向量中均无重复元素，提供一个简便的方法
A=[1 7 6 3 2 9 8];
B=[1 6 3 8] ;
setdiff(A,B)

晓夜寒风

Q8 内联函数如何部分赋值传参

1. tic
   4 t2 t7 L# _4 w! S
2. a=rand(1,10);
   
3. for i=1:length(a)
   
4. ff=['x.^2+',num2str(a(i))];
   0 J4 @6 x4 `; g/ T0 A  F
5. f=inline(ff,'x');
   & g/ l. b( L- y8 T
6. d(i)=quadl(f,-1,1);
   ( m$ X( C8 J1 s6 s# f$ b7 c
7. end
   
8. d
   4 m3 b% o- G  n8 H5 s" l( |! D5 f
9. toc

复制代码

晓夜寒风

Q13 我想写一个函数，输入参数是一个字符串比如是＂１６０＂，函数主要目的是完成数据存储，我现在的问题是存储文件时，文件命名的问题，比如我想把输入的参数作为文件名的一部分，比如希望输出的文件名为"temp１６０",这个该怎样实现啊，谢谢！！！i = 160;
  d5 v7 s% k: m" Y8 `name = ['temp' , num2str(i), '.txt']; %假设你想把当前工作区的数据存储为temp１６０.txt；
eval('save(name)');
. n3 V0 y5 i' V9 ?
$ D. M% ]7 P& I; c' C具体的请看eval save fopen fread fwrite fprintf等函数的帮助

晓夜寒风

Q10 [1，2，3，NaN，5，NaN，7，8，9，10]
6 ~) m; H2 ]! y& d( p) r如何将NaN以外的8个数加起来，然后求这8个数的平均数。
3 V! r# M, _* x# Y-----------------------------------------------------------------

>> a=[1 2 3 nan 5 nan 7 8 9 10];
5 n& \: G) M$ r5 y7 ?>> b=a;
/ H) G" Q9 U! h) e8 Q" S>> b(find(isnan(a)))=[]
6 ^/ ~, W+ K" \9 N& i
b =
- V0 ]! @! I& v3 e2 _! _; e
1 2 3 5 7 8 9 10
4 T; s/ ?& X1 e: L. J" l9 l# r  ^" s2 Y
* x0 t4 B7 |8 V$ Y; H5 h1 Y7 `. F
a=[1 2 3 nan 5 nan 7 8 9 10];
简洁些可以这样:
mean_a=mean(a(~isnan(a)))