【转载】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));
   . c5 }" Z5 m+ a) o1 y
5. disp('在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数:')
   9 m( E* M' o; |7 C: u& @
6. disp(d)

复制代码

在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数:
  23   128   233   338   443

晓夜寒风

matlab上标标注

1. subplot(1,1,1);
   
2. x=0:0.1:2*pi;
   
3. y1=sin(x);
   
4. y2=exp(-x);
     N( k* W- k8 C2 h
5. plot(x,y1,'--*',x,y2,':o');
   - D. U$ |* e5 O) L  f
6. xlabel('x=0 to 2\pi');
   
7. ylabel('values of sin(x) and e^{-x}');
   
8. title('function plots fo sin(x) and e^{-x}');
   
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)
   5 K( T% P: ~6 X6 p$ F7 E/ V4 W
2. 1)x=[-10:1:10];
   ; ?8 X+ _8 D" B. n
3. for i=1:21
   
4. f(i)=(x(i)-1)*(x(i)>1)+(x(i)+1)*(x(i)<-1);
   
5. end
   4 x, |0 N- T, W4 _# V. b0 j
6. plot(x,f,'r*')
   ) i4 @; y* e# N# e  C
7. 2)x=[-10:10];
   
8. f=zeros(size(x));
   
9. a=x>1;
   1 p8 F4 `5 h7 |  C8 U; S
10. b=x<-1;
    
11. f(a)=x(a)-1;
    
12. f(b)=x(b)+1;
    
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
( o0 P2 b# J! D5 D$ r怎么样把大于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);
   8 [+ l8 ?! t1 K% ]9 v( A5 i7 H7 t
4. x(k(1))

复制代码

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

1. b=a(1:k);
   6 m% n$ _7 v2 a/ G  y# }6 F3 n+ n
2. b=a(l:l+k); (l+k<=n)
   2 q7 D/ P2 H% ^. f" Y* H; @

复制代码

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   
   
3.   x=1
   
4. else   
   4 {' t& o9 J! A' m1 b
5.   x=0
   
6. end

复制代码

晓夜寒风

Q5: MATLAB中用pie命令画饼图

1. function PlotPie
   $ `. i4 Q7 G8 m1 v
2. x=[1 2 3 4 5 6 7 8 9 10];
   
3. y=[40 50 60 70 80 90 100 110 120 130 ];
   0 k' V4 A; L! n
4. str1={'A','B','C','D','E','F','G','H','I','J'};
   + P. L$ \" Z) T5 F
5. por=y./sum(y);
   
6. for i=1:length(por)   str{i}=[str1{i},' ',num2str(100*por(i)),'%'];
   
7. end
   % ?- D+ @  a  |% [
8. pie(x,y,str)

复制代码

晓夜寒风

Q6 非线性拟和y＝a*x^b+c的程序，求a,b,c.
function FitPaperData
6 f! n) k& i3 b; b5 qformat long
E1=rand(89,1);
" N/ J% s2 g6 VE2=rand(89,1)+3;
4 h/ w# L; U, }& _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]
应该怎么做
( E& [/ e  ~' {1 {# v如果各个向量中均无重复元素，提供一个简便的方法
A=[1 7 6 3 2 9 8];
B=[1 6 3 8] ;
setdiff(A,B)

晓夜寒风

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

5 c: M$ d) E% q6 {9 n

1. tic
   0 x" T6 p- n" s  h* c( `! m" S
2. a=rand(1,10);
   ) t& ]- p* t8 K/ a% ^$ w  ~
3. for i=1:length(a)
   / l- {' ~3 L5 a; E' t; X
4. ff=['x.^2+',num2str(a(i))];
   
5. f=inline(ff,'x');
   
6. d(i)=quadl(f,-1,1);
   
7. end
   ; D" {6 B) {6 ^" ?
8. d
   ) G6 a5 i! M# H2 u6 p  \" W) G) \
9. toc

复制代码

晓夜寒风

Q13 我想写一个函数，输入参数是一个字符串比如是＂１６０＂，函数主要目的是完成数据存储，我现在的问题是存储文件时，文件命名的问题，比如我想把输入的参数作为文件名的一部分，比如希望输出的文件名为"temp１６０",这个该怎样实现啊，谢谢！！！i = 160;
name = ['temp' , num2str(i), '.txt']; %假设你想把当前工作区的数据存储为temp１６０.txt；
eval('save(name)');
; j2 N# r( o2 P! f+ R
  C+ _( e' |( i具体的请看eval save fopen fread fwrite fprintf等函数的帮助

晓夜寒风

Q10 [1，2，3，NaN，5，NaN，7，8，9，10]
如何将NaN以外的8个数加起来，然后求这8个数的平均数。
7 N/ X$ O7 p7 ^* f6 |- ]6 F: a-----------------------------------------------------------------

' g4 j9 o5 y0 d6 g+ L5 I! x>> a=[1 2 3 nan 5 nan 7 8 9 10];
, s4 t3 V3 |5 w! m$ ~. }>> b=a;
* u, |( w' B+ g>> b(find(isnan(a)))=[]
1 }( S$ o& w6 S$ P" w) z, F
- V: T% j' p3 o2 Tb =

1 2 3 5 7 8 9 10

: i8 l% S5 g0 V, ~; n$ |
1 j9 X. V8 Z3 k' _a=[1 2 3 nan 5 nan 7 8 9 10];
简洁些可以这样:
mean_a=mean(a(~isnan(a)))